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[57] ABSTRACT 


A sample accommodator for the measurement of dielec- 
tric properties includes a conductive strip having two 
planar surfaces, and two ground planes, one parallel to 
each of the planar surfaces of the conductive strip. At 
each end, the strip passes through but is electrically 
isolated from an endplate which is electrically and phys- 
ically connected to each ground plane and which main- 
tains the relative positions of the ground planes. Dielec- 
tric properties are measured by placing a sample be- 
tween each face of the strip and the ground plane paral- 
lel to it. 


13 Claims, 6 Drawing Sheets 
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SAMPLE ACCOMMODATOR AND METHOD FOR 
THE MEASUREMENT OF DIELECTRIC 
PROPERTIES 


This application is a continuation of application Ser. 
No. 907,763, filed Sept. 15, 1986, abandoned, which 
application is, in turn, a continuation of application Ser. 
No. 637,376, filed Aug. 2, 1984 and now abandoned. 


BACKGROUND OF THE INVENTION 


The present invention relates in general to methods 
and apparatus for the measurement of dielectric proper- 
ties and in particular to methods and apparatus for the 
measurement of dielectric properties of core samples 
encountered in petroleum exploration. 

The recent advent of a dielectric well logging tech- 
nique in petroleum exploration has necessitated labora- 
tory measurement of the dielectric permittivity of core 
samples in the same way that resistivity and induction 
logging techniques necessitated laboratory measure- 
ment of the conductivity of core samples. Known tech- 
niques of laboratory dielectric permittivity measure- 
ments in the VHF (30-300 MHz) and the UHF 
(300-1000 MHz) ranges (roughly covering the 20-1100 
MHz range of frequencies of interest in dielectric well 
logging) require solid samples to be precision-machined 
to fit sample accommodators having precise dimen- 
sions. However, may types of core samples do not 
readily lend themselves to such machining. 


SUMMARY OF THE INVENTION 


Accordingly, apparatus according to the present in- 
vention involves a sample accommodator (i.e. a holder 
for the sample) for the measurement of dielectric prop- 
erties. The accommodator comprises two ground 
planes, each of which is positioned opposite to and 
parallel with a different one of two surfaces of a con- 
ductive strip. 

A method of measuring dielectric properties accord- 
ing to the present invention includes placing each of 
two sample portions between a different planar surface 
of a conductive strip and one of two ground planes. The 
scattering parameters of the sample portions are then 
measured either in the frequency domain or in the time 
domain. These parameters are then used to compute the 
real and imaginary parts of the complex dielectric con- 
stant, which represent the dielectric constant and the 
resistivity (or conductivity) of the sample, over a broad 
range of frequencies. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a perspective view of a prior art sample 
accommodator shown partially in phantom; 

FIG. 2 is an exploded view of a sample accommoda- 
tor according to the present invention; 

FIG. 3 is an exploded view of a transition for connec- 
tion of a coaxial cable to a center strip according to the 
present invention; 

FIG. 4 is a cross-sectional view along line C of FIG. 
2 through a sample accommodator containing a sample 
according to the present invention; 

FIG. 5 is a cross-sectional view along line C of FIG. 
2 through an empty sample accommodator according to 
the present invention; 

FIG. 6 is a two-port network diagram of the sample 
accommodator of FIG. 2; 
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FIGS. 7A-7D show a schematic, transmission line 
representation of the sample accommodator of FIG. 2; 

FIG. 8 is a view in top plan of a sample accommoda- 
tor having compartments for liquid sample containment 
according to the present invention; 

FIGS. 9A and 9B show a schematic transmission line 
representation the sample accommodator of FIG. 8; 

FIGS. 10A and 10B are a set of curves showing real 
and imaginary parts of the dielectric constant as mea- 
sured by the apparatus according to the present inven- 
tion; and 

FIG. 11 is a view in top plan of an alternative embodi- 
ment of the sample accommodator of FIG. 8 for moni- 
toring the dielectric constant of a constantly flowing 
fluid. . 


DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 


Of the techniques currently available for the measure- 
ment of machinable dielectric samples, the one that has 
found the widest application in the oil industry is the 
coaxial sample accommodator technique. 

As illustrated in FIG. 1, a coaxial sample accommo- 
dator according to the prior art requires a sample 14 to 
be machined in the form of an annular cylinder which 
has to conform to very precise dimensions and concen- 
tricity. While such requirements are not difficult to 
satisfy in the case of standard dielectrics, such as plexi- 
glass and the sort of fluorine-containing resiuns sold 
under the trademark Teflon, or, with some effort, in the 
case of well-consolidated rocks of low porosity and 
homogeneous constitution, it is very hard to satisfy 
these requirements for the large majority of core sample 
varieties commonly encountered in the oil industry. 
Furthermore, it has been shown that when deviations 
from an ideal annular cylindrical shape occurs, it is 
necessary to apply large corrections to measured data. 

Coaxial sample accommodators have an inner con- 
ductor 10 with a diameter a and an outer conductor 12 
with a diameter b, as shown in FIG. 1. Sample 14, hav- 
ing an axial dimension X, is coaxial with conductors 10 
and 12. Conductors 10 and 12 are connected at either 
end to a first connector 16 and a second connector 18. 

The coaxial sample accommodators currently in use 
derive their a and b dimensions from the standard 50- 
ohm connectors available in the market. Commonly, 
connectors are of the APC-7 and GR-900 types for 
which the b dimensions are 7 mm and 14 mm, respec- 
tively. These small dimensions make machining very 
difficult, and accentuate any errors in dimension. 

However, even when the b dimension is increased to 
as much as 20.65 mm, close to the largest practical size ~ 
at which undesirable wave propagation modes may be 
avoided, the concentricity requirement (the coinci- 
dence of the axes of the inner and the outer cylindrical 
surfaces) is very hard to satisfy for most samples of 
interest. This requirement for concentricity is easier to 
achieve for short sample lengths (in the axial direction), 
but shorter lengths introduce larger errors in the mea- 
surements of dielectric properties, and make the mea- 
surements uncertain at relatively lower frequencies. 

FIG. 2 illustrates a sample accommodator 200 ac- 
cording to the present invention. A conductive strip 20, 
having flat surfaces 202 and 20b (not shown), is sus- 
pended between a first ground plane 22, and a second 
ground plane 24. Although shown in an exploded view 
in FIG. 2, when sample accommodator 200 is closed 
around a sample, ground plane 22 is parallel to and 
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3 
opposite to surface 20a, as indicated in FIG. 4. Ground 
plane 24 is parallel to and opposite to surface 205 (not 
shown). 

Strip 20 has a first fluted end 20c (not shown) which 
forms a part of a first transition 270. Strip 20 also has a 
second fluted end 20d which forms a part of a second 
transition 290. 

Two sample portions, in the form of blocks A and B, 
are placed in sample accommodator 200. Each sample 
block is machined to approximate a rectangular prism. 
Each has a length dimension X parallel to a long dimen- 
sion of strip 20. Each has a width dimension Z perpen- 
dicular to the long dimension of strip 20 and approxi- 
mately coextensive with a dimension of ground plane 24 
which dimension is perpendicular to the long dimension 
of strip 20. Finally, each block has a depth dimension Y 
normal to flat surfaces 20a and 208. 

In order to contain samples A and B within sample 
accommodator 200, each ground plane is attached to a 
first end plate 28 at a first end and to a second end plate 
26 at a second end. This attachment is accomplished by 
tightening a screw (not shown) through a hole 28a in 
end plate 28 and into threaded hole 23a and in ground 
plane 22, respectively. Screws (not shown) are also 
tightened through the corresponding pairs of holes 28¢ 
and 23; 26a and 22a; 26b and 22h; 26c and 24c (not 
shown); 26d and 24d (not shown) 28c (not shown) and 
25c (not shown); and 28d and 25d (not shown). 

Threads on the outer casing of each of two tuning 
screws 220c and 220d are respectively engaged with 
threads in two holes 22c and 22d above opposite ends of 
strip 20. Likewise, below opposite ends of ground plane 
24, threads on tuning screws 240c (not shown) and 240d 
(not shown) are engaged with threads in holes 24¢ (not 
shown) and 24d (not shown). 

End plate 26 has a central, round aperture aligned 
with strip 20 and into which a part of first transition 270 
is fitted. Similarly, end plate 28 has a central, round 
aperture aligned with strip 20 and into which a part of 
second transition 290 is fitted. 

As shown in FIG. 3, wherein structures also found in 
FIG. 2 are identified by the reference numerals used to 
identify them therein, each transition, exemplified by 
transition 290, includes a center conductor contact 291 
which has a male threaded end 2914 for screwing into a 
female threaded end 292a of a connector center conduc- 
tor 292. 

Each fluted end of strip 20 has two projections form- 
ing a notch, such as notch 20¢ in fluted end 20d. A male 
end of center conductor 2926 is soldered to fluted end 
20d within notch 20e. The two projections of fluted end 
20d are placed between the ends of two opposing semi- 
cylindrical, notched fillers 295a and 295. The ends of 
fillers 295a and 2956 fit together around the projections 
of strip 20 so that the combined notches of fillers 295¢ 
and 2956 form a cylindrical aperture 295c aligned with 
notch 206. 

The combination of fillers 295a and 295 and fluted 
end 20d is inserted within the aperture in end plate 28. 
A retaining ring 294 is placed around fillers 295¢ and 
2955. A central cylindrical aperture in a cylindrical 
dielectric 293a is passed over end 292a so that connec- 
tor 293 may be fitted into the aperture in end plate 28. 
Ring 294 is held in place by gluing it to an inner surface 
of connector 293 using a suitable glue such as the one 
sold under the trademark Locktite. Conductor contact 
end 291a is tightened into female end 292a of center 
conductor 292 to complete assembly of transition 290. 
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The same steps are performed to assemble the same 
sorts of components into transition 270. 

Ground planes 22 and 24 and end plates 26 and 28 
may be readily machined out of aluminum blocks com- 
monly available to those skilled in the art. Connector 
parts suitable for use as elements 291, 292, 293 and 294 
are available as an MPC 14 (GR-900) connector, part 
number 2480A, from Maury Microwave Corporation, 
Cucamonga, California. The gradual “flute” expansion 
at ends 20c and 20d brings center strip 20 from the 
dimension of the inner conductor of the GR-900 con- 
nector to its full width w. Center strip 20 may be fabri- 
cated out of 0.0024 inch-thick brass sheet readily avail- 
able to those skilled in the art. For example, a strip 
constructed for a sample accommodator according to 
the present invention has an impedance of 50 ohms 
when empty and has a separation between ground 
planes of 0.6221 in. Fluted portions 20c and 20d of ra- 
dius 0.375 in. extend for 0.276 in. at either end of strip 20 
has a total length of 0.399 inches, a width of 0.847 in. 
and having a 0.197 inch deep notch at either end L/S 
band tuning screws suitable for use as elements 220c, 
220d, 240c and 240d may be obtained under part number 
JMC P/N 6965-0 from the Johanson Manufacturing 
Company, Boonton, NJ. These screws are adjusted to 
obtain a desired impedance characteristic over a desired 
frequency band when the sample accommodator is 
empty. 

Filers 295a and 295d are readily fabricated out of any 
of a number of fluorine-containing resins, such as those 
sold under the trademark Teflon, for example. In the 
case of a 50 ohms stripline, the dielectric constant of the 
material is approximately 2. 

The separation between ground planes 22 and 24 may 
be chosen to equal the inside diameter of the outer con- 
ductor of an unassembled commercial GR-900 rigid line 
connector, so as to obtain a good electrical match be- 
tween transitions 270 and 290 and strip 20. Strip 20 is 
chosen to be paper thin (i.e., t< <w), in order to elimi- 
nate one dimension from the calculations. The strip 
width w may be chosen to obtain a 50-ohm characteris- 
tic impedance for the empty sample accommodator. 
After calculating the theoretical value of w, a trial-and- 
error technique may be used to trim this width in order 
to fine-adjust the impedance. 

As shown in FIG. 2, samples A and B are cut in the 
form of bars of rectangular cross-section. Two such 
bars are used for each sampie. They fit in the sample 
accommodator in the way shown in FIG. 2. The X- 
dimension of each bar enters into the calculations. Thus, 
it is desirable that the X-dimension be as nearly the same 
for the two bars as possible. The actual value of this 
dimension is not critical, as long as it is known accu- 
rately. For each bar the y-dimension, corresponding to 
the separation between each ground plane and an op- 
posing surface of strip 20, has to be such that the sample 
fits snugly in the accommodator. The z-dimension of 
the sample accommodator, is somewhat arbitrary, but it 
is important that this dimension be large enough so that 
fringing fields 40, as shown in FIG. 4, are contained 
within the sample blocks. In FIG. 4, structures also 
found in FIG. 2 are identified by the reference numerals 
used to identify them therein. 

As shown in the partial cross-section of FIG. 5, 
wherein structures also found in FIG. 2 are identified by 
the reference numerals used to identify them therein, an 
empty sample accommodator according to the present 
invention is in the form of a transmission line, where the 
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central conductor is a metallic strip 20, sandwiched 
between two conducting planes 22 and 24. A dielectric 
material may be placed between each surface of strip 20 
and the plane opposite to it. Assuming that the width w 
of strip 20 is large enough so that the fringing fields at 
the edges of strip 20 do not interact, the characteristic 
impedance Z, a dielectric-filled line is: 


94.15 
wl, Cf 
T— (@/b * 0.0885, 


= @) 


Zo —== ohms 


._ where the capacitance Crdue to fringing fields is given 
by 


Q) 


‘}} 


ee? Se | 2 pee Ka = 
{[; (77) fol 1-Wbh* | 


1 1 7 
[on Je 1 — Wop 


in picofarad/cm and where: 

€-=the permittivity of the dielectric; 
w=the width of strip 20; 

t=the thickness of strip 20; and 

b=the distance between planes 22 and 24. 

Dielectric constants are obtained according to the 
method of the present invention from measurements of 
scattering parameters of a sample by connecting transi- 
tions 270 and 290 to ports of an appropriate apparatus. 
In the frequency domain measurements may be made by 
using a network analyzer, such as the Model 8505A 
Automatic Network Analyzer System 300 available 
from Hewlett-Packard, Inc., Palo Alto, CA. In the 
time-domain a waveform processing system, such as the 
Model 7854/7812 Time-Domain Reflectometer avail- 
able from Tektronix, Inc., Beaverton, Oregon, may be 
used:to make measurements with various pulses genera- 
tor/sampler combinations. Either technique measures 
the S parameters (scattering parameters) S11 and S21 (or, 
equivalently, S22 and S2). These are complex quantities 
giving the amplitudes and phases of the signals reflected 
from and transmitted through the sampie, respectively. 
From these parameters, the complex dielectric permit- 
tivity (as well as the magnetic permeability) of the sam- 
ple can be calculated, given the width of the sample in 
the direction of propagation which corresponds to the 
X-dimension in FIG. 2. 

The network formed by the sample accommodator 
and attached measurement apparatus may be dia- 
grammed as shown in FIG. 6. 

The two-port network illustrated in FIG. 6 may be 
described by the following equations: 


5, =S11 a1 +512 a2 G) 


b2=S21 a1 +522 a2 (4) 

where: 

ai=a voltage measured across source impedance Zs 
looking into a two-port network 60 at 61; 

a2=a voltage measured across a impedance Z; looking 
into network 60 at 62; 

bi =a source voltage measured looking out of network 
60 at 64; and 
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b2=a voltage measured looking out of network 60 at 63 
and where S11, S22, S12, and S21 are the S parameters 
(scattering parameters) which are defined as follows: 


(5) 


an input reflection coefficient with the output 
port terminated in a matched load (Zz = Zo sets 


az = 0); 


(6) 


322 = = 
42 Ig1=0 


an output reflection coefficient with the input 
port terminated in a matched load (Z; = Zo and 


V; = O set a) = 0 where V; is a source voltage); 


(1) 


a forward transmission gain with the output port 
terminated in a matched load; and 


(8) 


a\=0 


a reverse transmission gain with the input port 
terminated in a matched load. 


FIG. 7A shows a diagrammatic, edge-on view of two 
rectangular slabs of rock inserted into a strip transmis- 
sion line. The equivalent transmission line model of this 
view is shown in FIG. 7B. FIG. 7C depicts the trans- 
mission line terminated in a matched load and FIG. 7D 
is an equivalent circuit for the circuit in FIG. 7C. An air 
region has a magnetic permeability jz, and electric per- 
Mittivity €, while a rock region has a permeability 2 
and a permittivity ¢,2. This permittivity is a complex 
quantity defined by 

: 


(9) 


oO 


€c2 = €0 («2 —I Ge 

where: 

€,2=the relative dielectric constant of a rock sample, 

o2=the conductivity of the rock sample, and 

«@==the angular frequency of an electromagnetic wave 
propagating through the line. 

In the equivalent transmission line model of FIG. 7B, 
Zo and Z2 are the impedances of the air portion and 
rock-filled portions of the line, respectively, and k2 the 
propagation constant. These quantities are defined as 
follows: 


(10) 


ee 83,988 ae 
wb Cf 
T— (/b + T0885 
(14) 
N12€0 83,988 
= nT 2 
22 oa w/b Ceo ohms; and 
T— G/b) * 00885 ep 
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-continued 


(12) 
Q=a \ H2€c2 

where: 

€.2=the permittivity of the sample; 

Cr=the fringe capacitance given by 


2 1 
(an | +1|- 


am tae ae 
T= Wy ~3 inf 1 woe 


(13) 


“2 
in picofarads/cm; 
é€=the permittivity of the medium in the stripline sam- 
ple accommodator, either €, for the Zo portion or 
€9€72 for the Z2 portion; and 
where all other constants and variables are as defined 
above. 

For a two-port network we need only consider the 
input reflection coefficient, S11, and the forward trans- 
mission gain, S21 (or, alternatively, S12 and S22). Since 
they are defined with the output terminated in a 
matched load, we consider the case illustrated in FIG. 
7C where the input impedance Zin (W) between W= —1 
and W=0 is given by 


0.0885 ¢€ 


Cr = eg 


: 1+ Ro9e2k2¥ (14) 
in (FW) = 1 Rape 
where 
Zo — Z2 (15) 


Ro = By FZ 


which is the reflection coefficient of the line at W=0. 
W=the distance from the edge at which an incident 
wave leaves the sample in the direction of propaga- 
tion (as indicated by the arrow in FIG. 7B); 
l=the length of the path through the sample in the 
direction of propagation; and 
where all other constants and variables are as defined 
above. 
Using the definition of S1; in Equation (5), the model 
in FIG. 7C can be simplified to that shown in FIG. 7D 
so that the expression for S1; becomes 


Zin(—) — Zo 
Zinl—) + Zo 


16 
Su = (16) 


where all constants and variables are as defined above. 
Inserting Equation (14) into equation (16) with 
W=-1 gives 


Roz + Rage #2! (17) 
Si = ~~ or 

1 + Ro2Re- 
where 

22 — Zo (18) 
Ro2 = Zt zo 3 and 


where all other constants and variables are as defined 
above. 
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8 
If the incident wave has an amplitude of unity and its 
phase at W = —1 is zero, then from FIG. 7C, the voltage 
V1 between W=—1 and W=— o is 
Vi(W) =e HAW) +81) oko +h (19) 
and the voltage V2 between W=0 and W= —1 is 
VW) =Vte_ ja 4 V— ok (20) 


where all constants and variables are as defined above. 
Upon satisfying the boundary condition at W=—1, 
Vi(—)D=V2(—)) we have 


14+S11=V to?! (1-+Ro9° Hh 21) 


where 
R29=V-/V+, (b 22) 
where 
pe ee Ie es (23) 
TT + Rage var 


where all constants and variables are as defined above. 
From the definition of S21, 


¥2(0) 
21 = 1 


(24) 


= V+ + V- = Vt + Rp). 


The “1” in the denominator represents the unit ampli- 
tude and zero phase of the incident wave and V2(O) is 
the voltage measured at W=0. All other constants and 
variables are as defined above. Now we substituted 
Equations (22) and (23) into Equation (24) and obtain 


spy = LE Sue AL + Rao) es) 
a 1 + Rye Jk2! 
combining Equation (17) with Equation (25) yields 
(1 — R2pje~se2! (26) 
Sy = Toe enor < 
1—- Rogge 2k! 
We can simplify this expression by letting 
q=e—/*2!, and (27) 
Zz — Zo (28) 
R= R= ~Ru = Wy: 


Consequently, Equations (17) and (26) may be rewritten 
as 


se RQ @) (29) 
1 — Rg? 

and 

see (1 — Rg (30) 
ren R2q2 . 


Adding Equations (29) and (30) we obtain 
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Su + Sa =H (31) 
Solving for q gives 5 
_ Gut+sy-R (32) 
~ 1— Sit + SR 
Equations (29) and (30) also give 10 
Sut+¢S+1=4 te (33) 
Su + Sy -1 = B52 4) 15 
Si — Sati= + eGo (35) 
Si) —- S21 -1= Sa = 36) 20 


which are used to write 


(G7) 25 
E + S21 + aE — Sa+  |- 1+ 2R + R2 
Si +52 —1 JL Su —sa-t J” 1—-2R + R2 
or, 
2 2g, ,2 30 
$11R*—(S117—S21*+ DR+S11=0. (38) 
Solving for R we have 
‘ (83, — 3, + 1) -A G9) 35 
= 2811 
where 
A=[(S11+Sa1+ 1)(S11+S21—) 
(S1—S21+ 1)S11-Sa— DF (40) 40 


and where all other constants and variables are as de- 
fined above. 

The quantity q is obtained by substituting Equation 
(=) into Equation (34) such that 


= (Sh) — Sp + —A ey 


2821 
Using Equation (27) we have 


ip = i (42) 


and from Equation (28) we have 
(43) 


n= otek. 


If we recall that ko =0V p2€02, we can use Equation 
(40) to write 


(44) 
Wine? = tn 


Substituting Equations (10) and (11) into Equation (43) 
gives 
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€c2 = €0 1-—R 
Ee -\= P+R° 


The complex permittivity of the rock, €,2, is found from 
Equations (44) and (45) to be 


o- (VE ole] 


where R and q are complex numbers. Comparison of 
this result with Equation 9 gives: 
€2 = 1 —- Re [ amo ( )I and 
Ninoco 
o2= — | < + Im [ano (+4)} 


Thus, q and R are determined from the S parameters 
using Equations (39) and (41), and the dielectric con- 
stant and conductivity are calculated from Equations 
(47) and (48). ; 

As illustrated in FIG. 8, for measurement of liquid 
samples, hollow compartments 80 and 81 may be made 
of materials sold under the trademarks made of standard 
dielectric materials (e.g., Teflon and Rexolite) and hav- 
ing a rectangular cross-section may be used in place of 
the solid sample bars shown in FIG. 2. Each of these 
compartments accommodates an equal portion of the 
sample and is equipped with inlet and outlet holes (not 
shown) for introducing liquids into the compartments. 
By knowing the dielectric permittivity and the wall 
thickness (in the direction of propagation 82) of the 
compartment, its effect may be removed from the mea- 
surement, and thus, the dielectric permittivity of the 
liquid derived. The filler and drainage holes (not 
shown) in the compartments also allow the insertion of 
a temperature probe (not shown). 

The measurement technique for solid samples may be 
readily extended to liquid samples. In this case liquid 
compartments are placed parallel to each planar surface 
of strip 20. A transmitted wave must now pass through 
the container’s walls, one wall before entering the sam- 
ple and one after leaving it as indicated in FIG. 8. End 
walls 83 and 84 of compartment 80 are away from the 
stripline so that they do not interfere with wave propa- 
gation. 

In FIG. 9A, 1 and €; are the magnetic permeability 
and dielectric permittivity, respectively, of the walls of 
the liquid container. These values can be measured 
using the dry sample theory where a sample made of the 
same material as the liquid container is measured in the 
stripline sample accommodator. The permeability and 


(45) 


- 


(46) 


(47) 


1-R 
1+R 


(48) 


60 permittivity of the liquid sample are denoted by p22 and 


€¢2; respectively. 

FIG. 9B is the equivalent transmission line model of 
FIG. 9A. The values k; and Z; are the propagation 
constant and characteristic impedance of the wall por- 
tion of the sample accommodator. They can be obtained 
from dielectric measurements of the wall material. 

As the impedances and propagation constants may be 
expressed as follows: 
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11 
Vine (49) 
5 a Se 
Veet (w/b) ( Seo ) 
T — /b) 0.0885 €-1 
(50) 
k=@ H1€cl 
(51) 
he N2¢0 83,988 
Me (w/b) Creo 
2 any + Woes en 
(2) 
=o Vinee 
where: 
b=the distance between ground plates; 
w=the width of strip 20; 
t=thickness of strip 20; 
€,1=the permittivity of the liquid container; 
pi=the permeability of the liquid container; 
€,.2=the permittivity of the liquid sample; 
}42=the permeability of the liquid sample; 
=the angular frequency of the wave; and 
Cris the fringe capacitance given by 
(63) 
0.0885 ¢ 2 1 
c=—s ( T— Wb) n(; Wey * 1)- 


1 es: Seen 
(—= = i)a(s — (1/6) ')} 


in picofarads/cm where € is €,1 for the Z; portion and 
€72 for the Z2 portion. 

The line is divided into five regions and W==d is 
taken as a reference point where the incident wave 
(which is assumed to have an amplitude of 1) has a 
phase of zero as in FIG. 9. W is the distance from the 
boundary of at which an incident wave enters a liquid 
sample in the direction of propagation (as shown by 
arrow 82 in FIG. 8). If d is the wall thickness of liquid 
container 80 and | is the length of the liquid sample in 
the direction of propagation, then the voltages and 
currents in the five regions may be written as follows: 


Region 1: W S-—d 


Vo(W) = e~FOW+4) 4. 5) ek W +4) (54) 
IW) = a [eA +4) _ 51 elk +0)) (55) 
Region 2: ~d = WO 

V\(W) = Ae—JRIW 4. Belk W (56) 
1M) = de — Bolt on 
Region 3:05 W=1 

Wy(W) = Ce~HR2W 4. Dejk2W (58) 
IW) = Jo [CeO — Delta @) 
Region 4: /2 WS!l+d 

V3(W) = BeWJKIW 4 FeikiW (60) 
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-continued 


1 [ze—ikIW _ Reiki 


Zi 


(61) 


(W) = 


Region 5:/+d3 W 


Va(W) = Spe KW 1-4) (62) 


IW) = al [Sare#O—I—-a om 


where A through F are constant coefficients and S71 
and S11 are measured scattering parameters. 

At each interface, the boundary conditions consisting 
of continuity in the phase and amplitude of the voltages 
and currents must be satisfied. Consequently, at 
W=-—d we have 


VA{—d)=Vi(—d) (64) 
If—d)=I(—d) (65) 
From Equations (54) and (57) 
1 + Sy; = Aokld 4. Be—ikld (66) 
1— Sit Agikid _ Be—jkid (67) 
“ZX  — +4 
Solving these equations for A and B gives us 
(68) 
; Z| 
A = he~Jhld [a + Si) + - sy | 
(69) 
kid ci 
B= je*l4) (1 + Sy) — (1 — Si) Zo 
Similarly, at the point W=1+d, we have 
(70) 
E = jokll+4) 5) € + 3) 
(71) 
F = jeH1+9) $9) ( - 3) 
For W=0 
(72) 
22 v4) 
cma[ (+ )e+-#)e | 
(73) 
22 22 
D= LC: ot yy +(1 +z Je | 
For W = | we have 
(74) 


, 22 , 22 
= lelkal — —jkul — AL 
C = bel [G + Zi ac + (: Zi, ee | 
. v4) P 22 ‘ 
= tect ee ER —jkil Joe kil 
D=the [G Zi ez + ( + Zi ) ] 


Now the product CD obtained from Equations (72) and 
(73) must equal that obtained from Equations (74) and 
(75). 


(75) 
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13 14 
76 84) 
2 en Zo _,| Sc2 | BO : 
Leo pee 2, N pm 0 
: thus 
[(:-#)e 0-4) +04) ]- 9 
Zi Z 41 7 €c2 Zo s 
“nm Zr ON 
22 = 10 
[G + Zi Jez sas Combining Equation (82) with Equation (85) gives 
(86) 
(: = Be tt |G - 2) pe + iin _20 <9 
1 Zi 15 ae ae ig 


We can simplify this equation by letting 20 
X=Z7/Z1 
‘ G=Ee=/kil 
; 25 
H=FeAi! 
Thus Equation (76) becomes 
(U+%)4+(1—X) BY —X- . 30 
)4+(14+X)B]=(1+4)G+(1-XY ANA -~2G+C 
14+X)H] (77) 
After some rearrangement of this expression we have: 
(G—-H)—(4—BY 1X? =(G+ H)—-(A+ BP (78) 35 
Consequently, 
79 
yo 22 [Gt -Gin as 
41 (4 — BY —(G— HY 


since the right hand sides of Equations 72 and 74 must 
be equal, we can write 
45 


_ a G+(1— XH (80) 


= ew ikl 
iia ED 9 EES 

The quantities A, B, G, H, and X are knowa, so q is 
known. The propagation constant of the strip transmis- 50 
sion line filled with liquid is: 


. (81) 
n=a N H2€c2 -in@ 


55 
so 
- (82) 
Ving = Ua 
60 


Recalling that in air the characteristic impedance for a 
stripline transmission line is given by 


(83) 


83,988 65 


w/b + of 
1 — (t/b) 0.0885 


Zo = yo 


and combining this with Equation (51) leads to 


From the basic definition of complex permittivity given 
in Equation (9), the dielectric constant and conductivity 
for any liquid can be obtained from Equation (86) as 


follows: 
[: In @) ] 


Zo 
Im [: In (g) 2] 


FIG. 10 shows data on Berea sandstones, saturated 
with distilled water and with brines of various salinities. 
These data are in good agreement with published re- 
sults. The critical sample dimension is the x-dimension 
shown in FIG. 2. As indicated in Equation (44), the 
calculated permittivity is inversely proportional to this 
dimension. However, this dimension can be machined 
and measured precisely enough so that it is not a signifi- 
cant source of error. 

Based on the above considerations, and also on expe- 


(87) 


1 
€2 = +- Re 


BOEO 


Veo 4 


pO i 


Zo 
Z2 


(88) 


= 


40 rience with measurements of standard solids and liquids, 


it is believed that the accuracy in the measurement of 
the complex permittivity for a prototype sample accom- 
modator is typically +5% in the range 20-1300 MHz, 

degrading to +10% for very lossy samples (e.g. 
strong brines and some alcohols). It is believed here is 
room for further improvement in the design (e.g. the 
transition), leading to even higher accuracies over 
broader bands. 

While the present invention has been described in 
terms of a preferred embodiment, further modifications 
and improvements will occur to those skilled in the art. 
For example, sample compartments, the apparatus ac- 
cording to the present invention, may be modified to 
measure the dielectric constant of a flowing fluid. 

In FIG. 11, sections of pipe segments 110 and 111 (not 
shown) are rectangular cross-section similar to contain- 
ers 80 and 81. Each of segments 110 and 111 is con- 
nected to an in-line mixer 113 by way of small diameter 
pipes 112. In-line mixer 113 has an input connected to a 
feed line 114 which may be connected to a source of a 
flowing fluid. A return line is connected to an end of 
each of segments 110 and 111 which end is directly 
opposite to the end connected to lines 112. By connect- 
ing transitions 270 and 290, as described for the embodi- 
ment of the present invention having sample containers 
80 and 81, the complex dielectric constant of the flow- 
ing fluid may be continually monitored. Such apparatus 
may be used, for example, to monitor salinity of a flow- 
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15 
ing brine (such as produced water) or to measure the 
water/oil fraction in a small diameter pipeline. 

I desire it to be understood, therefore, that the present 
invention is not limited to the particular form shown 
and that I intend in the appended claims to cover all 
such equivalent variations which come within the scope 
of the invention as claimed. 

What is claimed is: 

1. Stripline apparatus for the measurement of scatter- 
ing parameters of a core sample comprising: 

(a) a sample holder comprising: 

a single conductive strip having two planar sur- 
faces; and 

two ground planes, each of said ground planes 
being positioned opposite to and parallel to a 
different one of said planar surfaces of said strip, 

wherein said holder has a characteristic impedance 

unchanging along its length when said holder is 

empty of samples; 

(b) means for containing the core sample in two por- 
tions between each of said ground planes and said 
surface opposite to it, said sample portions having a 
z dimension large enough so that the main electric 
field and the fringing electric field from said con- 
ductive strip are substantially contained within the 
portions; 

(c) means for measuring amplitudes and phase shifts 
of the signals reflected from and transmitted 
through the core sample. 

2. Apparatus as recited in claim 1 wherein said means 
for containing comprises two endplates, each of said 
endplates being connected to each ground plane and 
being electrically isolated from said strip. 

3. Apparatus as recited in claim 2 wherein said means 
for containing further comprises two coaxial cable con- 
nectors, each coaxial cable connector being connected 
to different end of said strip. 

4. Apparatus as recited in claim 3 wherein each end of 
said strip is fluted. 

5. Apparatus as recited in claim 2 wherein said means 
for containing further comprises two compartments, 
each of said compartments being positioned between a 
different surface of said strip and one of said ground 
planes. 
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6. Apparatus as recited in claim 5 wherein said means 
for containing further comprises two coaxial cable con- 
nectors, each coaxial cable connector being connected 
to different end of said strip. 

7. Apparatus as recited in claim 6 wherein each end of 
said strip is fluted. 

8. Apparatus as recited in claim 7 wherein said means 
for containing further comprises two pipe segments 
compartments, each of said pipe segments being posi- 
tioned between a different surface of said strip and one 
of aid ground planes. 

9. Apparatus as recited in claim 8 wherein said means 
for containing further comprises two coaxial cable con- 
nectors, each coaxial cable connectors being connected 
to different end of said strip. 

10. Apparatus as recited in claim 9 wherein each end 
of said strip is fluted. 

11. A method for measuring scattering parameters of 
core samples comprising the steps of: 

(a) placing each of two core sample portions between 

a different planar surface of a conductive strip and 
one of two ground planes, wherein said conductive 
strip and two ground planes form a sample holder 
having: a characteristic impedance unchanging - 
along its length when said sample holder is empty 
of samples, said core samples having dimensions 
large enough so that the main electric field and the 
fringing fields from said conducive strip are con- 
tained within said core sample portions; and 

(b) calculating the amplitude and phase shift parame- 

ters of said core sample portions in the frequency 
domain based on reflected and transmitted signal 
data from the samples. 

12. The method as recited in claim 11 wherein said 
placing step comprises the step of introducing each core 
sample portion in the form of a block of rectangular 
cross-section between a different planar surface of the 
conductive strip and one of the two ground planes. 

13. The method as recited in claim 11 wherein said 
placing step comprises the step of introducing each core 
sample portion in liquid form into one of two compart- 
ments, each compartment being located between a dif- 
ferent planar surface of the conductive strip and one of 


the two ground planes. 
* * * * * 
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METHOD AND APPARATUS FOR BROADBAND 
MEASUREMENT OF DIELECTRIC PROPERTIES 


FIELD OF THE INVENTION 


The present invention relates generally to the analysis 
of materials having dissimilar dielectric properties. 
More specifically, this invention provides an antenna 
which can combine the functions of various resistivity 
and dielectric constant measuring devices into a single 
tool, capable of operating over a wide range of frequen- 
cies. It is particularly useful in the field of well logging. 


BACKGROUND OF THE INVENTION 


It is well known to log or record certain electrical 
characteristics of earth formations surrounding a well 
borehole as a function of depth in order to determine 
the location and extent of oil-bearing strata and to quan- 
tify the amount of oil present in such strata. A log of 
formation resistivity versus depth may indicate the pres- 
ence of hydrocarbons, since hydrocarbon-bearing for- 
mations typically exhibit a higher resistivity than forma- 
tions containing mostly salt water. There are only three 
material parameters which affect an electromagnetic 
wave. They are conductivity (resistivity), magnetic 
permeability and dielectric constant. Conductivity pro- 
vides an indication of the energy absorbing characteris- 
tics of the medium, while magnetic permeability and 
dielectric constant give a measure of the energy storing 
capacity of a material. The magnetic permeability of 
most earth materials is the same and is equal to the 
magnetic permeability of free space. It is therefore of 
very little use in electrical logging techniques. 

It is well known that conductivity or resistivity 
(which is the reciprocal of the conductivity) has wide 
variation in value for earth materials and strongly af- 
fects electromagnetic waves. A propagating electro- 
magnetic wave has two fundamental characteristics, 
amplitude and phase. By comparing the amplitude and 
phase of an electromagnetic wave as it passes receivers, 
propagation characteristics due to formation may be 
studied. Measurement of these two characteristics (or 
equivalently, the wave travel time and attenuation) may 
be used to determine the dielectric constant and the 
resistivity of the media through which the wave is prop- 
agated. 

Dielectric oil well logging is aimed at determining 
the water saturation and the water salinity in the forma- 
tion from measurements of the dielectric constant and 
the resistivity of the formation. A number of criteria 
enter into the selection of the frequency of the electro- 
magnetic wave used to probe the formation. Depending 
on these criteria, service companies have fielded a num- 
ber of different logging tools at different frequencies. 
However, no one tool in the prior art is capable of 
probing a formation over a broad band of frequencies. It 
is therefore advantageous to extend the frequency range 
to include the function of the dielectric and the resistiv- 
ity (induction) logging tools. 

The objective of dielectric logging is to measure the 
dielectric constant and the conductivity of earth forma- 
tions at a specific frequency, and deduce therefrom the 
water saturation and the salinity of this water. It has 
been found that logging at various frequency domains 
has various specific advantages. For example, logging 
at relatively high frequencies, near 1000 MHz, one 
probes a depth of only a few inches into the formation, 
so that here one measures essentially the properties of 
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the invaded zone (i.e., the zone where the formation 
fluid has been displaced by the fluid filtrate from the 
mud in the borehole). On the other hand, logging at 
relatively low frequencies, near 20 MHz, allows one to 
probe much deeper into the formation (~ a few feet), 
yielding perhaps information on the virgin formation. 
Logging at intermediate frequencies causes one to 
probe intermediate depths. The low frequency tools are 
of the centralized mandrel type, since at these frequen- 
cies the loss of energy in the mud annulus between the 
tool and the formation can be tolerated. At high fre- 
quencies, however, this loss would be prohibitive, and 
hence the high frequency tools must necessarily be of 
the pad type, where the pad directly contacts the forma- 
tion or borehole wall. 

The development of tools at these many different 
frequencies shows that it is advantageous to create a 
single tool capable of logging: 

(i) at a number of discrete frequencies 

(ii) over a continuous range of frequencies employing 

the swept frequency technique, and 

(iii) in the time-domain, resulting - upon inversion - in 

the same information as in (ii) above. 

The main advantages of such a tool would be: 

(i) it would yield the profile of water saturation and 

water salinity with distance from the wellbore, and 

(ii) since a large amount of information would be 

gathered, this makes it possible to deduce, from the 
dielectric log alone, the porosity, in addition to the 
water saturation, and the water salinity. 

U.S. Pat. No. 4,899,894, issued to Katahara et al., 
describes an acoustic logging tool, and does not discuss 
electromagnetic logging. It uses a plurality of acoustic 
frequencies only to accommodate a wide variety of 
formations having different preferred logging frequen- 
cies, and does not teach a method of broadband logging 
over each formation. 

U.S Pat. No. 3,982,176, issued to Meador, discloses a 
conventional induction log (20 KHz) and a conven- 
tional dielectric log (16 MHz) in single tool, using dif- 
ferent antennas (coils) for the different frequencies. This 
is a mandrel type tool, that hangs freely in a wellbore 
with an annulus of mud around it. There is no mention 
of broadband logging. 

U.S. Pat. No. 4,451,789, also issued to Meador, dis- 
closes a method of depth-probing the formation by 
changing the transmitter-receiver spacing. Frequency 
variation is not discussed. Although Meador uses three 
frequencies, they are very closely spaced, and are not 
used for the purpose of broadband logging. Meador 
simply uses slightly different frequencies as a way of 
labeling the transmitter-receiver pairs It is therefore 
essentially a single frequency tool. It is also a mandrel 
type tool, which hangs freely in a wellbore. 

USS. Pat. No. 4,774,471, issued to Sims et al., discloses 
a- single-frequency (between 10 MHz and 200 MHz) 
mandrel type tool Sims does not discuss how their an- 
tenna can accommodate this frequency range, or 
whether different antennas are used for different fre- 
quencies. This patent only uses a wider range of fre- 


‘quencies within the operation of an earlier tool. The 


patent does teach that broadband dielectric logging is 
useful, however. 

The largest hurdle to developing a broadband dielec- 
tric logging too] has been the lack of a suitable broad- 
band antenna that can couple electromagnetic energy to 
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and from a formation, and that is compact enough to fit 
within the confines of a logging tool. 

An additional advantage is that a dielectric measuring 
apparatus can be applied to the field of medical technol- 
ogy. For example, U.S. Pat. No. 4,240,445 issued to 
Iskander et al. teaches a method of coupling electro- 
magnetic energy into a material such as tissue, to mea- 
sure water content. Measuring lung water content is an 
especially useful application. However, Iskander’s de- 
vice is so large that only a few antennas can be place on 
the chest, and the antenna cannot be described as a point 
source. Also, the electric field vanishes at some distance 
from the antenna, as the electric fields in the two paral- 
Jel slots are oppositely directed. Furthermore, a resistor 
is included in the antenna, which dissipates much of the 
electromagnetic energy in the antenna itself. Additional 
prior work includes M. F. Iskander and C. H. Durney 
(1980): “Electromagnetic Techniques for Medical Di- 
agnosis: A Review”, Proceedings of IEEE, vol. 68, no. 
1. and M. F. Iskander et al (1982): “Two-dimensional 
Technique to Calculate the EM Power Deposition Pat- 
tern in the Human Body”, Journal of Microwave 
Power, vol. 17, no. 3. 

The prior work is limited in the attempts at broad- 
band logging (or measuring) in that no suitable single 
antenna element has been designed which can couple 
electromagnetic energy into a material, whether it be a 
geologic formation or tissue, over a broad range of 
frequencies, that is also sufficiently compact and is ca- 
pable of handling high power levels. There is therefore 
a need for a device and a method for use in such broad- 
band logging/measuring/heating applications. 


SUMMARY OF THE INVENTION 


The present invention is surprisingly successful in 
providing a method and apparatus for combining the 
functions of various resistivity and dielectric constant 
measuring devices and electromagnetic energy cou- 
pling device into a single tool, capable of operating over 
a wide range of frequencies. It is especially useful in 
well logging applications. The apparatus is capable of 
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taking broadband measurements of materials having 


dissimilar dielectric properties. 

A measuring tool such as a logging tool, having a tool 
face, also has a novel transmitting and a novel receiving 
antenna. Electromagnetic energy is transmitted to a 
transmitting antenna. A stripline adapter permits the 
energy to flow to a stripline having a metallic central 
strip. A strip face of the central strip is bent at approxi- 
mately right angles, and has a height that is compatible 
with desired frequency coverage. 

A ground plane extends from the stripline adapter to 
the right angle bend, so that a distal end of the central 
strip extends away from it, and a void is created be- 
tween the center strip and the ground plane. 

A dielectric is positioned to nearly fill the void. The 
dielectric is comprised of a material having a very high 
dielectric constant and a very low energy loss. The 
transmitting antenna is positioned so that the ground 
plane is fixedly connected to the measuring tool, and the 
strip face lies flush with the tool face, so that electro- 
magnetic energy can be transmitted into the material to 
be analyzed. ; 

An enclosure surrounding the stripline is comprised 
of four metallic walls which are positioned in electrical 
contact with the ground plane and the stripline adapter, 
so that the strip face is nearly centered in the opening 
created by the walls and the ground plane. 
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A loss-less, non-conducting material fills in any re- 
maining open space in the enclosure, so that the non- 
conducting material forms an additional wall that is 
really flat with the strip face. 

A receiving antenna is comprised in essentially the 
same manner as the transmitting antenna, and is posi- 
tioned in the tool so that it can receive the electromag- 
netic energy which has traveled through the material 
being probed. A means for monitoring the received 
energy detects changes as the apparatus is moved across 
the material. 

In another embodiment of this invention, broadband 
measurements are taken to determine the nature of a 
fluid in a material such as a geologic formation. 

It is one object of this invention that electromagnetic 
energy is transmitted and received over a wide fre- 
quency range, specifically from 2 KHz to 1 GHz. 

The tool may further comprise a pad, which con- 
forms to the inside of a wellbore and holds the antennas. 
At least one transmitting and one receiving antenna are 
necessary. A plurality of each is often desirable. Some 
antennas may be staggered in a spiral fashion, to prevent 
undesirable coupling between adjacent antennas. 

Fluids analyzed include water and hydrocarbons. 
Electromagnetic energy is monitored to provide an 
indication of the porosity, salinity, and water saturation 
of the material to be analyzed. Depth probing a geo- 
logic formation is a particularly useful application, to 
measure the nature of fluids as a function of the distance 
from the antennas. 

The above and other embodiments, objects, advan- 
tages, and features of the invention will become more 
readily apparent from the following detailed descrip- 
tion of the invention, which is provided in connection 
with the accompanying drawings. 


DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a schematic, sectional view of the inventive 
logging (or measuring) tool positioned in a wellbore. 

FIG. 2 shows a top, front, and side view of the novel 
transmitting antenna. 

FIG. 2A is the same view as FIG. 2, further illustrat- 
ing the enclosing metallic walls. 

FIG. 3 shows an antenna mounted on a tool face. 

FIG. 4A-4E shows schematic side views of the in- 
ventive apparatus as applied to well logging technol- 
ogy. 

FIG. 5A-5E show three graphs of transmission and 
return loss as a function of frequency. 

FIG. 6 is a graph of transmission and return loss as a 
function of frequency, for low frequencies. 

FIG. 7A-7D show four graphs of time-domain trans- 
mission measurements at various distances from a metal 
refiector plate in a brine. 


DETAILED DESCRIPTION OF THE 
~ INVENTION 


‘In accordance with the present invention, a new im- 
proved method and apparatus for determining the na- 
ture of various materials and the fluids contained 
therein, using a broadband measuring apparatus, has 
been developed. 

Referring to the drawings, a first embodiment of the 
inventive broadband logging (or measuring) tool 101 is 
shown in FIG. 1, positioned in a wellbore 103. The 
wellbore wall 105 is sometimes lined with a layer of 
mud cake 107. Pressure arm 109 positions tool face 111 
near wellbore wall 105, so that transmitting antennas 
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such as T; and T2 and receiving antennas such as Ry and 
R2 are positioned close to the mud cake 107 or the well- 
bore wall 105. The transmitting antennas, such as T,, 
and the receiving antennas, such as R, are mounted on 
the tool face 111 above and below the tool pad 119, and 
are positioned opposite the wallbore 105, and are sepa- 
rated from It by the mud cake 107 and the mud in the 
wellbore 103. The tool face 111 is defined as the outer 
surface of the logging tool 101. 

A material, in this case a geologic formation 113 
surrounds the wellbore. An invaded zone 115 represents 
the region of the formation 113 that has been permeated 
by the mud filtrate The outer region 117 is that part of 
the formation 113 that has not been permeated by the 
mud filtrate In a formation evaluation program, it is 
desirable to know the properties of both regions, as a 
function of distance from the antenna or tool face 111. 
This may be called depth-probing. 

In FIG. 1, the logging tool 101 has a pad 119 or pad- 
type device so that some or all of the antennas can be 
positioned very close to the formation so that electro- 
magnetic energy can by transmitted through the forma- 
tion. The antennas work in pairs. That is, antennas 
T1-T4 and T,, are transmitting 1 antennas, and antennas 
R)-R4and R, are receiving antennas. Antenna T} trans- 
mits electromagnetic energy though the formation, and 
antennas Rj and R>2 through R,, receive it. Antenna T2 
transmits electromagnetic energy to antennas R) and R2 
through R,,, and so on, through antennas T,. The anten- 
nas.can all be positioned in a vertical line, as shown, or 
they can be staggered, in the azimuthal direction. It is 
possible for any given antenna to alternately play the 
role of a transmitting or a receiving antenna. 

For any pair of antennas, T;R;, by measuring the 
phases and amplitudes of the received power at Rj and 
of the reflected power at Tj, the dielectric constant and 
the conductivity of the region probed by the pair T;R; 
can be determined. From the dielectric constant and the 
conductivity, water saturation and salinity can be calcu- 
lated using theoretical expressions well known in the 
art. It is desirable that the antennas can transmit and 
receive energy over the entire frequency range of inter- 
est. For dielectric well logging, a frequency range of 
approximately 10-1000 MHz is desirable. To include 
resistivity and induction logging, this range can be ex- 
tended down to a few kilohertz. 

The depth of investigation for a pair such as T3R3 is 
a function of the frequency of the probing radiation. As 
this frequency is gradually reduced from a high value, 
the depth increases until at some frequency f; it equals 
the depth dj, roughly obtained by drawing a semicircle 
with T3R3 as the diameter-as seen on FIG. 1 (dashed 
lines). Decreasing the frequency further does not in- 
crease the depth of investigation much further Thus, 
each pair T;R; probes a zone of depth dj; appropriate to 
that pair and to the frequency f; for that pair. In this 
way, a set of broadband antennas placed with different 
spacings permits information about a set of zones of 
varying depth to be obtained (the longer the spacing, 
the lower is f;). The above discussion also illustrates that 
for a given transmitter-receiver pair, some depth-prob- 
ing is possible, by changing the frequency. Thus, by 
having different transmitter-receiver spacings and the 
broadband capability, redundant depth-probing is possi- 
ble. 

Another use of the first embodiment incorporates the 
concept of standalone logging. In conventional dielec- 
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tivity of a zone is measured at a single frequency, one 
needs an independent knowledge of the porosity of the 
zone in order to convert this information to water satu- 
ration and salinity. This porosity is usually derived from 
a separate logging tool. However, if the same zone 
could be probed in two very distinct frequencies, then 
the dual frequency data can be inverted to simulta- 
neously calculate the porosity, water saturation and 
salinity, without the need of a separate logging tool. 

Referring to FIG. 1 for the pair T4R4, if the fre- 
quency is raised from f;1 (which is the frequency cor- 
responding to the maximum depth of investigation dj+1 
for the pair), the zone of investigation will shrink from 
dj+1 until at some frequency F;+, this zone (dotted line 
120) will have a depth equal d;. Thus, having broadband 
antennas with different spacings permits us to probe the 
same depth d; with two different frequencies. 

An example of the inventive transmitting antenna 150 
is shown in FIG. 2. A coaxial connecting means, such as 
coaxial connector 151 is electrically connected to a 
stripline adapter 153 which is capable of transmitting 
electromagnetic energy from the coaxial connector 151 
to a stripline section with metallic central strip 155. An 
especially useful stripline adapter is a model No. 
3070-1404-10 designed by Omni-Spectra, or other types 
of microwave stripline adapters. Other types of trans- 
mission means may be utilized to transmit electromag- 
netic energy to the antenna. For example, a strip trans- 
mission line may be electrically connected to the strip- 
line section having the metallic central strip 155. As a 
commercial coaxial-to-stripline transition means has 
been utilized, the dimensions included herein reflect this 
means. One knowledgeable in the art would realize that 
the dimensions may be altered to change frequency 
coverage and to fine-tune performance. 

Metallic center strip 155 has a front end 157, a flat 
strip body 159, a flat strip face 161, and a distal end 163. 
The front end 157 is electrically connected to the center 
conductor 169 of the stripline adapter 153. Solder is a 
particularly useful connecting means. Flat strip body 
159 may also be tapered to come to a point at front end 
157 to provide a smooth electrical transition between 
the center conductor 169 and the center strip 155. The 
strip face 161 is bent at approximately right angles to 
strip body 159, and has a height that is measured from 
the right angle bend to distal end 163. The height is 
compatible with the desired frequency coverage. The 
longer the height, the more lower frequency coverage 
is allowed. A 3" height permits a frequency range of 
approximately 2 KHz—1 GHz. The metallic center 
strip 155 can be made of any metal. Copper, brass, or 
aluminum are especially useful. 

A ground plane 165 extends from stripline adapter 
153 to the right angle bend in the center strip 155, so 
that the distal end 163 extends away from the ground 
plane 165 and so that a void exists between the center 
strip 155 and the ground plane 165. Ground plane 165 is 
comprised of a metal. Commercial grade stainless steel 
is particularly useful. It is desirable to keep the ground 
plane and center strip as short at possible, to permit the 
apparatus to remain as compact as possible and to allow 
the use of as many antennas as possible. 

The void between the ground plane 165 and the cen- 
ter strip 155 is largely filled with a dielectric 167. The 
dielectric 167 should have a very high dielectric con- 
stant and a very low loss. By loss, we mean the dissipa- 
tion of energy. The dielectric 167 can be a ceramic 
dielectric, and comprised of material such as Barium 
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Titanate or Lead Zirconate Titanate. A crystalline di- 
electric may also be used, although more expensive. 
The thickness of the dielectric 167 is determined by the 
stripline adapter 153 used. The dielectric 167 acts to 
make the capacitance of the center strip 155 very large. 

The construction of the antenna is completed by 
enclosing the center strip 155 by metallic walls 181, 182, 
183, and 184, which contact the ground plane 165 and 
the adapter 153 electrically, as shown in FIG. 2A. The 
walls add rigidity and prevent leakage of the electro- 
magnetic radiation. The strip face 161 is approximately 
centered in the rectangular opening created by the 
edges of the walls and the edge of the ground plane 165. 
Thus, the: distance between an edge of the strip face 161 
and the adjacent edge of a wall is substantially the thick- 
ness of the dielectric 167. The entire void space in the 
antenna enclosed by the walls, including the set back 
168 at the dielectric edge, is filled with a loss-less, non- 
conducting material such as a mixture of epoxy and 
alumina which sets hard, seals the antenna, and makes it 
more rugged. 

The ground plane 165 and the walls 181, 182, 183, and 
184 are fixedly connected to a logging (or measuring) 
device or tool as seen in FIG. 3. The strip face 161 is 
positioned to lie flush with the tool face 171 (or a tool 
pad), so that the transmitting antenna 150 can transmit 
electromagnetic energy into a material such as a geo- 
logic formation or mammal tissue. Void space 173 is 
filled with a loss-less, non-conducting material such as 
an epoxy-aluminum compound. The ground plane 165 
and the walls 181, 182, and 183 connect to the tool face 
(or tool pad) and the center strip 161 becomes the sensor 
“button”. 

A receiving electromagnetic antenna is comprised in 
essentially the same manner as the transmitting antenna, 
and is positioned in the logging tool in the same manner 
as the transmitting antenna, so that the receiving an- 
tenna can receive the electromagnetic energy which has 
traveled through the material that is analyzed. 

The logging tool also contains a means for monitor- 
ing the received electromagnetic energy for changes, as 
the tool is moved across the material and moves across 
the interface of materials having dissimilar dielectric 
properties so that the interface can be located, and the 
nature of the material can be evaluated. An example of 
frequently studied dielectrically dissimilar materials are 
hydrocarbons and water, both of which are normally 
contained in the pore spaces of rocks. 

The logging tool is capable of transmitting and re- 
ceiving electromagnetic energy that has a range of fre- 
quencies. The range of frequencies is from about 2 KHz 
to about 1 GHz. In this embodiment, the electromag- 
netic energy is of sufficient strength to travel through 
the material (or formation) to the receiving antenna as a 
propagating wave. The amplitude and the phase of the 
received signal contain information (dielectric constant 
and conductivity) about the formation. A duplicate 
receiving antenna is usually employed, thereby obtain- 
ing measurements that eliminate certain unknown fac- 
tors. Additionally, by employing a duplicate transmit- 
ting location above or below the pair, the effects of 
alternate upward and downward transmission are aver- 
aged to achieve “borehole compensation”, known in 
the art. In a second embodiment, the logging (or mea- 
suring) tool (a device) determines the nature of a fluid in 
a material. The material may be, for example, a geologic 
formation or may be a mammal tisstie. In this embodi- 
ment, the device is constructed as in the first embodi- 
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ment, except that the means for monitoring the electro- 
magnetic energy for changes as the apparatus is moved 
across the material determines the nature of the fluid. 
For a geologic formation, the fluids analyzed can be 
hydrocarbons and water. 

In either the first or second embodiment, the electro- 
magnetic energy can be monitored to provide an indica- 
tion of the salinity and the water saturation of the mate- 
rial. If the same zone could be probed in two very dis- 
tinct frequencies, then the dual frequency data can be 
inverted to simultaneously calculate the porosity, also. 

The apparatus transmits electromagnetic pulses of 
designated profile shapes into the material to be ana- 
lyzed, such as a geologic formation. The pulses are 
received by antennas spaced at different distances from 
the transmitter. The different receiver spacings corre- 
spond to different depths of investigation. When the 
received pulses are analyzed, one obtains a broadband 
dielectric spectrum for different depths in the forma- 
tion. Together, all this information yields the depth- 
dependent water saturation, salinity, porosity and ma- 
trix dielectric constant (matrix identification). The 
depth-dependent water saturation could then be ana- 
lyzed to give an estimate of the permeability. The other 
conventional parameters such as formation factor can 
also be derived from the data, if there is present in the 
well, a zone with 100% water saturation. 

The present invention may operate either in time or 
frequency domain, or both. The above paragraph de- 
scribes the application in the time domain. Because of 
the compact size of the invention, many such sensors 
can be installed on a tool, to comprise an array of sen- 
sors capable of “depth-probing” the formation (i.e., 
measuring the saturation and the salinity as a function of 
distance from the wellbore). The depth-probing and the 
wide frequency coverage aspects can possibly be com- 
bined to give the tool a stand-alone character: It can, in 
and of itself, determine the porosity, the salinity, and the 
water saturation without the need of an additional po- 
rosity tool (which is required in the present art). 

FIGS. 4A-4E illustrate various applications of the 
inventive apparatus as applied to well logging technol- 
ogy. The antennas (both receiving and transmitting) can 
be mounted on a pad 201 or pad-type device which is 
designed to substantially conform to the inside of the 
wellbore 203, as seen in FIGS. 4B, 4D, and 4E. The pad 
acts to reduce or prevent propagation of the electro- 
magnetic energy through the lossy mud 205. At low 
frequencies, however, the loss in the mud is tolerable 
and the antennas can be mounted on the tool face 207, 
on a centralized tool, as shown in FIGS. 4A and 4C. A 
combination of two arrangements is also possible, as 
seen in FIG. 4E, where antennas are positioned both on 
tool face 207 and on a pad 201 or pad-type device to 
include low frequency waves. The vertical spacings and 
azimuthal orientations are selected from theoretical 
considerations known in the art. Therefore, the device 
may comprise a single transmitting antenna and one or 
more receiving antennas, or it may comprise a plurality 
of transmitting antennas. The antennas may be oriented 
in several arrangements. They may be staggered about 
the tool face and/or pad or may be staggered in a spiral 
fashion. Staggering of the antennas minimizes unwanted 
coupling between adjacent antennas, while still permit- 
ting a short spacing distance between antennas. 

It is advantageous to use several antennas at pre- 
selected spacings. This provides a variable depth of 
investigation, as one transmitting antenna could be used 
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with a plurality of receiving antennas. A variable fre- 
quency of investigation is also possible, as larger spac- 
ings are required with larger wavelengths. Also, redun- 
dant spacings provide enhanced quality control. 

The apparatus can operate in the frequency domain, 
using a single frequency, multiple frequencies (such as 
simultaneous, selectable, or time-multiplexed for exam- 
ple), or swept frequency techniques Or, the apparatus 
can operate in the time domain, using pulses of a wide 
variety of shapes, widths, rise and fall times, etc. When 
the pulses are transformed to the frequency domain, 
either electronically using a spectrum analyzer, or nu- 
merically using mathematical transforms, the same in- 
formation is obtained as would be given by a frequency 
domain tool. 


THEORY OF TIME-DOMAIN DIELECTRIC 
LOGGING 


A time-domain tool eliminates much of the cumber- 
some electronics of a swept frequency or a multiple 
frequency tool, and shifts the burden of the electronics 
to software. Consider a too! with a broadband transmit- 
ting antenna located at Z=0 and two identical broad- 
band receiving antennas located at Z; and Z2 (Zo, Z1 
and Z2 being three points along the axis of the well- 
bore). Suppose that an electromagnetic pulse f (t, 0) of 
a finite duration is launched into the formation at Z=0 
beginning at t=0 (f (t, 0) =O for t=0 and t>to, say). 

The Fourier frequency components of this pulse are 
given by: 


1 q) 
Fw, 0) = a f Kit, 0) edt 
oO 


where w is the circular frequency of the radiation. Simi- 
larly, the pulses f (t, Z1) and f(t, Z2) received at Z; and 
Z2 can be decomposed into their spectral components: 


1 a ' 2) 
Fo, Z)) = a= [Re 21) ema 
o 


1 a ; (3) 
F(w, Z2) = a i Kt, Z2) edt 
0 


Now, the relationship between the transmitted and re- 
ceived Fourier components is 


F(, Z})=F(o, 0) e*21 M(w) G(w, Z1) (4) 
F(w, Z2)=F(@, 0) e*22 M(w) G(w, Z2) (5) 
where 
k=a+if is the complex wave number and 
a= phase shift of the wave in radians/meter 


B= attenuation of the wave in nepers/meter 

M(w) = factor resulting from mismatches at the 
transmitting antenna/formation and receiving an- 
tenna/formation interfaces 

G(w, Z)= factor resulting from geometric spread of 
radiation. 

Dividing (5) by (4) eliminates the term M(@), and we 

have 
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(6) 
Fla, Z2) 


Qs G(w, 22) Rel 
™\ Fo, 2) 


B BZ Go,Z\) Im 


whence the relative dielectric constant € and the con- 
ductivity o are obtained using the relations: 


ee (7) 


€optow? 


c= 2aB. mho/m 8) 
Bow 


where a, 8, € and o are all function of frequency. 

The geometric factor in Equation (6) has to be deter- 
mined experimentally by pulsing the tool in air and in 
other lossy media of known, lossy dielectric properties, 
such as brine. 

In the alternative to Equation (6), the dielectric prop- 
erties may also be derived by referencing to only one 
receiver, and using equations (4) or (5). Here, the quan- 
tity M(w) is determined by measuring the power re- 
flected back into the transmitting antenna due to the 
mismatch, as described in U.S. Pat. No. 4,831,331 issued 
to De and Keetch. 

A prototype logging tool pad was constructed, with 
the inventive antennas. The pad consists of one trans- 
mitting and one receiving antenna, the distance between 
them being variable. 

An acceptable dielectric logging antenna must meet 
the following criteria: 

(i) It must be able to couple sufficient energy into and 
from the formation at its operating frequency to 
allow probing of the formation; 

(ii) This probing energy must penetrate into the for- 
mation, rather than clinging to the surface of the 
tool pad (i.e., it must travel as a freely propagating 
wave rather than a surface wave guided along the 
pad). 

In the present instance, the above two conditions must 
hold over the entire range of the frequency of opera- 
tion. 

The first of the above criteria is tested by measuring 
the return loss for the transmitting antenna, and the 
transmission loss from the transmitting to the receiving 
antenna - both as a function of frequency. These mea- 
surements are shown in FIG. 5 where the antenna pad 
is placed in air and against brine of conductivity 0.5 
mho/m (to represent a formation). The return loss 
curve in brine shows that sufficient energy is entering 
the brine over the frequency range of the measuring 
device (Hewlett-Packard HP8505A Network Analyzer; 
500 KHz - 1300 MHz) to permit probing. The transmis- 
sion loss shows that sufficient energy is being received 
at the receiving antenna to permit measurements. 

Measurements were made by using another measur- 
ing device (HP3577A Network Analyzer, SHz ~ 200 
MHz) to test the low frequency limitation of the an- 
tenna. The results are shown in FIG. 6, showing that 
the low frequency limitation is about 5 KHz. The im- 
proved return loss performance in the 200 MHz region 
(at FIG. 5) results from a drying (curing) of the epoxy 
alumina filling between measurements. 

FIG. 7 shows time-domain transmission measure- 
ments at various distances (d) to a metal reflector plate 
in the brine The change in amplitude of the received 
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pulse as a function of the distance of the metallic reflec- 
tor shows that the energy has penetrated into the brine 
out to the location of the plate. 

While a preferred embodiment of the invention has 
been described and illustrated, it should be apparent 
that many modifications can be made thereto without 
departing from the spirit or scope of the invention. 
Accordingly, the invention is not limited by the forego- 
ing description, but is only limited by the scope of the 
claims appended hereto. 

What is claimed is: 

1. Apparatus for making broadband measurements of 
materials having dissimilar dielectric properties com- 
prising an enclosed measuring tool having a tool face, 
said measuring tool further comprising an electromag- 
netic transmitting antenna, said transmitting antenna 
further comprising: 

(a) a coaxial cable connecting means and means to 

transmit electromagnetic energy therethrough; 

(b) a stripline adapter capable of transmitting electro- 
magnetic energy from said coaxial cable connect- 
ing means to a stripline having a metallic central 
strip, said center strip having a front end, a flat strip 
body, a flat strip face, and a distal end, said front 
end electrically connected to a center conductor of 
said stripline adapter, said strip face bent at approx- 
imately right angles to said strip body and having a 
height measured from said right angle bend to said 
distal end that is compatible with a desired fre- 
quency coverage; 

(c) a ground plane which extends from said stripline 
adapter to said right angle bend, so that said distal 
end extends away from said ground plane and so 
that a void exists between said center strip and said 
ground plane; 

(d) a dielectric largely filling said void, said dielectric 
comprised of a material having a very high dielec- 
tric constant and a very low energy loss, so that 
said transmitting antenna is positioned so that said 
ground plane is fixedly connected to said measur- 
ing tool and said strip face is positioned to lie flush 
with said tool face so that said transmitting antenna 
can transmit electromagnetic energy into said ma- 
terial; 

(e) an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in electri- 
cal contact with said ground plane and said strip- 
line adapter, so that said strip face is nearly cen- 
tered in the opening created by said walls and said 
ground plane; 

(f) a loss-less, non-conducting material which fills in 
any remaining open space in said enclosure so that 
said non-conducting material forms an additional 
wall that is nearly flat with said strip face; 

(g) a receiving electromagnetic antenna comprised in 
essentially the same manner as said transmitting 
antenna, said receiving antenna positioned in said 
measuring tool in the same manner as said transmit- 
ting antenna, so that said receiving antenna can 
receive said electromagnetic energy which has 
traveled through said material; and 

(h) means for monitoring said received electromag- 
netic energy for changes as said apparatus is moved 
across said materials whereby said dielectric prop- 
erties can be measured. 

2. Apparatus as recited in claim 1 wherein said trans- 

mitting antenna transmits and said receiving antenna 
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receives electromagnetic energy over a frequency 
range of 2 KHz to 1 GHz. 

3. Apparatus as recited in claim 1 further comprising 
a means for positioning said tool face near said material. 

4. Apparatus as recited in claim 1 wherein a strip 
transmission line is electrically connected to said strip- 
line, so that electromagnetic energy can be transmitted 
thereto. 

5. Apparatus as recited in claim 1 wherein said anten- 
nas are positioned on a pad-type device. 

6. Apparatus as recited in claim 1, further comprise 
a plurality of receiving antennas. 

7. Apparatus as recited in claim 6 further comprising 
a plurality of transmitting antennas. 

8. Apparatus as recited in claim 7 wherein said mate- 
rials having dissimilar dielectric properties are hydro- 
carbons and water. 

9. Apparatus as recited in claim 1, wherein said 
broadband measurements are taken to determine said 
dielectric properties as a function of distance from said 
antenna. 

10. An apparatus for making broadband measure- 
ments to determine the nature of a fluid in a material, 
said apparatus having a tool face and further comprising 
a first electromagnetic transmitting antenna, said first 
transmitting antenna further comprising: 

(a) a coaxial cable connecting means and means to 

transmit electromagnetic energy therethrough; 

(b) a stripline adapter capable of transmitting electro- 
magnetic energy from said coaxial cable connect- 
ing means to a stripline having a metallic central 
strip, said center strip having a front end, a flat strip 
body, a flat strip face, and a distal end, said front 
end electrically connected to a center conductor of 
said stripline adapter, said strip face bent at approx- 
imately right angles to said strip body and having a 
height measured from said right angle bend to said 
distal end that is compatible with a desired fre- 
quency coverage; 

(c) a ground plane which extends from said stripline 
adapter to said right angle bend, so that said distal 
end extends away from said ground plane and so 
that a void exists between said center strip and said 
ground plane; 

(d) a dielectric filling most of said void, said dielectric 
composed of a material having a very high dielec- 
tric constant and a very low energy loss, so that 
said first transmitting antenna is positioned so that 
said ground plane is fixedly connected to said log- 
ging tool and said strip face is positioned to lie flush 
with said tool face so that said first transmitting 
antenna can transmit electromagnetic energy into 
said material; 

(e) an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in electri- 
cal contact with said ground plane and said strip- 
line adapter, so that said strip face is nearly cen- 
tered in the opening created by said walls and said 
ground plane; 

(f) a loss-less, non-conducting material which fills in 
any remaining open space in said enclosure so that 
said non-conducting material forms an additional 
wall that is nearly flat with said strip face; 

(g) a receiving electromagnetic antenna comprised in 
essentially the same manner as said transmitting 
antenna, said receiving antenna positioned in said 
apparatus in the same manner as said transmitting 
antenna, so that said receiving antenna can receive 
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said electromagnetic energy which has traveled 
through said material; and 

(h) means for monitoring said electromagnetic energy 
for changes as said apparatus is moved across said 
material, so that the nature of said fluid can be 
determined. 

11. Apparatus as recited in claim 10 wherein said 
transmitting antenna transmits and said receiving an- 
tenna receiver electromagnetic energy over a frequency 
range of 2 KHz to 1 GHz. 

12. Apparatus as recited in claim 1 or 10 wherein said 
transmitting antenna can alternately function as a re- 
ceiving antenna and said receiving antenna can alter- 
nately function as a transmitting antenna. 

13. Apparatus as recited in claim 11 further compris- 
ing a pad, said pad substantially conforming to the in- 
side of a wellbore and holding said transmitting and 
receiving antennas. 

14. Apparatus as recited in claim 13 further compris- 
ing a plurality of receiving antennas. 

15. Apparatus as recited in claim 14 further compris- 
ing a plurality of transmitting antennas. 

16. Apparatus as recited in claim 15 wherein some of 
said antennas are positioned on said tool face and some 
antennas are positioned on said pad. 

17. Apparatus as recited in claim 15 wherein some of 
said antennas are staggered. 

18. Apparatus as recited in claim 15 wherein said 
antennas are staggered in a spiral fashion. 

19. Apparatus as recited in claim 10 wherein said 
fluids are hydrocarbons and water. 

20. Apparatus as recited in claim 1 or 10 wherein said 
ground plane is no greater than 10 mm in length. 

21. Apparatus as recited in claim 1 or 10 wherein said 
strip face has a height that is no greater than 5 mm. 

22. Apparatus as recited in claim 1 or 10 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the porosity of said material. 

23. Apparatus as recited in claim 1 or 10 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the water saturation of said material. 

24. Apparatus as recited in claim 1 or 10 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the salinity of said fluid. 

25. Apparatus as recited in claim 1 or 10 wherein a 
strip transmission line is electrically connected to said 
stripline, so that electromagnetic energy can be trans- 
mitted thereto. 

26. Apparatus as recited in claim 10 wherein said 
nature of said fluid is determined as a function of the 
distance from said antenna. 

27. Method for making broadband measurements of 
materials having dissimilar dielectric properties com- 
prising the steps of: forming a measuring tool having a 
tool face, an electromagnetic transmitting antenna and a 
receiving antenna, said transmitting antenna further 
comprising: 

(a) a coaxial cable connecting means and means to 

transmit electromagnetic energy therethrough; 

(b) a stripline adapter capable of transmitting electro- 
magnetic energy from said coaxial cable connect- 
ing means to a stripline having a metallic central 
strip, said center strip having a front end, a flat strip 
body, a flat strip face, and a distal end, said front 
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distal end that is compatible with a desired fre- 
quency coverage; 

(c) a ground plane which extends from said stripline 
adapter to said right angle bend, so that said distal 
end extends away from said ground plane and so 
that a void exists between said center strip and said 
ground plane; 

(d) a dielectric largely filling said void, said dielectric 
comprised of a material having a very high dielec- 
tric constant and a very low energy loss, so that 
said transmitting antenna is positioned so that said 
ground plane is fixedly connected to said measur- 
ing tool and said strip face is positioned to lie flush 
with said tool face so that said transmitting antenna 
can transmit electromagnetic energy into said ma- 
terial; 

(e) an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in electri- 
cal contact with said ground plane and said strip- 
line adapter, so that said strip face is nearly cen- 
tered in the opening created by said walls and said 
ground plane; “a 

(f) a loss-less, non-conducting material which fills in 
any remaining open space in said enclosure so that 
said non-conducting material forms an additional 
wall that is nearly flat with said strip face; 

(g) said receiving antenna comprised in essentially the 
same manner as said transmitting antenna, and posi- 
tioned in said measuring tool in the same manner as 
said transmitting antenna, so that said receiving 
antenna receives said electromagnetic energy 
which has traveled through said material; 

interconnecting said measuring tool with a means for 
monitoring said electromagnetic energy for changes as 
said measuring tool is moved across said materials 
whereby said dielectric properties can be measured; 
interconnecting said measuring tool with a source of 
electromagnetic energy; and 

moving said measuring tool across said materials. 

28. Method as recited in claim 27 wherein said trans- 
mitting antenna transmits and said receiving antenna 
receives electromagnetic energy over a frequency 
range of 2 KHz to 1 GHz. 

29. Method as recited in claim 27 further comprising 
a means for positioning said tool face near said material. 

30. Method as recited in claim 27 wherein a strip 
transmission line is electrically connected to said strip- 
line, so that electromagnetic energy can be transmitted 
thereto. 

31. Method as recited in claim 27 wherein said anten- 
nas are positioned on a pad-type device. 

32. Apparatus as recited in claim 27 further compris- 
ing a plurality of receiving antennas. 

33. Apparatus as recited in claim 32 further compris- 
ing a plurality of transmitting antennas. 

34. Apparatus as recited in claim 27 wherein said 
materials having dissimilar dielectric properties are 
hydrocarbons and water. 

35. Method as recited in claim 27 wherein said broad- 
band measurements are taken to determine said dielec- 


_ tric properties as a function of distance from said an- 


end electrically connected to a center conductor of 65 


said stripline adapter, said strip face bent at approx- 
imately right angles to said strip body and having a 
height measured from said right angle bend to said 


tenna. 

36. Method for making broadband measurements to 
determine the nature of a fluid in a material, comprising 
the steps of: forming an apparatus having a tool face, an 
electromagnetic transmitting antenna, and a receiving 
antenna, said transmitting antenna further comprising: 
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(a) a coaxial cable connecting means and means to 
transmit electromagnetic energy therethrough; 

(b) a stripline adapter capable of transmitting electro- 
magnetic energy from said coaxial cable connect- 
ing means to a stripline having a metallic central 
strip, said center strip having a front end, a flat strip 
body, a flat strip face, and a distal end, said front 
end electrically connected to a center conductor of 
said stripline adapter, said strip face bent at approx- 
imately right angles to said strip body and having a 
height measured from said right angle bend to said 
distal end that is compatible with a desired fre- 
quency coverage; 

(c) a ground plane which extends from said stripline 
adapter to said right angle bend, so that said distal 
end extends away from said ground plane and so 
that a void exists between said center strip and said 
ground plane; 

(d) a dielectric largely filling said void, said dielectric 
comprised of a material having a very high dielec- 
tric constant and a very low energy loss, so that 
said transmitting antenna is positioned so that said 
ground plane is fixedly connected to said measur- 
ing tool and said strip face is positioned to lie flush 
with said tool face so that said transmitting antenna 
can transmit electromagnetic energy into said ma- 
terial; 

(e) an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in electri- 
cal contact with said ground plane and said strip- 
line adapter, so that said strip face is nearly cen- 
tered in the opening created by said walls and said 
ground plane; 

(f) a loss-less, non-conducting material which fills in 
any remaining open space in said enclosure so that 
said non-conducting material forms an additional 
wall that is nearly flat with said strip face; 

(g) said receiving antenna comprised in essentially the 
same manner as said transmitting antenna, and posi- 
tioned in said apparatus in the same manner as said 
transmitting antenna, so that said receiving antenna 
receives said electromagnetic energy which has 
traveled through said material; 

interconnecting said measuring tool with a means for 
monitoring said electromagnetic energy for changes as 
said apparatus is moved across said material whereby 
said nature of said fluid can be determined; 
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interconnecting said apparatus with a source of electro- 
magnetic energy; and 
moving said apparatus across said material. 

37. Method as recited in claim 36 wherein said trans- 
mitting antenna transmits and said receiving antenna 
receiver electromagnetic energy over a frequency 
range of 2 KHz to 1 GHz. 

38. Method as recited in claim 27 or 36 wherein said 
transmitting antenna can alternately function as a re- 
ceiving antenna and said receiving antenna can alter- 
nately function as a transmitting antenna. 

39. Method as recited in claim 37 further comprising 
a pad, said pad substantially conforming to the inside of 
a wellbore and holding said transmitting and receiving 
antennas. : 

40. Method as recited in claim 39 further comprising 
a plurality of receiving antennas. 

41. Method as recited in claim 40 further comprising 
a plurality of transmitting antennas. 

42. Method as recited in claim 41 wherein some of 
said antennas are positioned on said tool face and some 
antennas are positioned on said pad. 

43. Method as recited in claim 40 wherein said anten- 
nas are staggered in a spiral fashion. 

44. Method as recited in claim 37 wherein said fluids 
are hydrocarbons and water. 

45. Method as recited in claim 27 or 36 wherein said 
ground plane is no greater than 10 mm in length. 

46. Method as recited in claim 27 or 36 wherein said 
strip face has a height that is no greater than 5 mm. 

47. Method as recited in claim 27 or 36 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the porosity of said formation. 

48. Method as recited in claim 27 or 36 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the water saturation of said formation. 

49. Method as recited in claim 27 or 36 wherein said 
electromagnetic energy is monitored to provide an indi- 
cation of the salinity of said formation. 

50. Method as recited in claim 27 or 36 wherein a strip 
transmission line is electrically connected to said strip- 
line, so that electromagnetic energy can be transmitted 
thereto. 

51. Method as recited in claim 36 wherein said nature 
of said fluid is determined as a function of the distance 


from said antenna. 
* * * * * 
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MICROWAVE SYSTEM FOR MONITORING 
WATER CONTENT IN A PETROLEUM PIPELINE 


FIELD OF THE INVENTION 


The present invention relates to the field of two- 
phase measurement..More specifically, the present in- 
vention relates to measurement of the water content of 
petroleum in a petroleum pipeline. It can also be applied 
to multi-phase measurement of other material. 


BACKGROUND OF THE INVENTION 


Crude oil produced from subterranean wells is fre- 
quently contaminated with undesirable materials. The 
most common of these materials is water. The fluid 
produced from an oil well can contain anywhere from 
less than 1% water to 99% water or more. Produced 
fluids containing these high “water cuts” are treated to 
remove most or all of the water because crude oil sales 
contracts typically allow a maximum of | to 3% water. 

After treatment, the crude oil stream is monitored to 
determine the amount of water remaining in the crude. 
Various methods of monitoring the water content of a 
crude oil stream have been proposed. For example, U.S. 
Pat. No. 4,596,136 describes a system in which a sample 
of oil is removed. This system (and all others which 
require the removal of a sample) assume that the water 
content of the sample is representative of the water 
content of the entire flow stream. This is frequently not 
a valid assumption, particularly in large pipelines which 
contain significant concentrations of water. 

Other. methods have also been proposed. For exam- 
ple, spinners, venturies, or capacitance probes have 
been used to determine the water content of crude oil. 
Among other disadvantages, all of these systems require 
the insertion of one or more devices into the crude oil 
stream. These intrusive systems can significantly in- 
crease the pressure drop in a pipeline, requiring addi- 
tional pumping capacity. Many of these devices also 
assume that the conditions in one portion of the pipeline 
are representative of the entire pipeline. This again may 
not be the case. 

In summary, there is a need to provide a method and 
means for determining the water content of flowing 
crude oil streams which is highly accurate, nonintru- 
sive, and which does not require the removal of a sam- 
ple that is representative of the entire crude oil stream. 


BRIEF SUMMARY OF THE INVENTION 


The present invention provides a nonintrusive 
method and means for accurately determining the water 
content of crude oil flowing in a pipeline. The device is 
comprised of two antenna pots mounted flush in the 
walls of a petroleum pipeline. In each of the ports are 
various S- and X-band antennas. These antennas are 
used to derive the real and imaginary parts of the com- 
plex dielectric constant of the oil/water mix (e*) from 
which the volume fraction of water ($) can ultimately 
be determined. 

In the first port is an S-band (or other low frequency) 
transmitting antenna and in the second port is an S-band 
receiving antenna. Since an ideal mixing formula for the 
complex dielectric constant in the low frequency region 
is not known, the volume fraction of water cannot be 
calculated directly. Therefore, an “adjustment factor” 
(p) must be calculated. In addition, it is necessary to 
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know the conductivity of the water (ay) to calculate } 
when e* is known. 

The adjustment factor and water conductivity are 
determined by utilizing a set of low frequency (S-band) 
and (X-band) transmitting and receiving antennas con- 
tained entirely within one of the antenna ports. Across 
the short distance of the antenna port, the attenuation 
and the phase shift of a transmitted signal in either the 
X- or S-band can be determined unambiguously. From 
the measurments, the value of €;* and €,*, i.e., the value 
of e* in the S- and X-band for the mixture near the wall, 
are determined. Using €,*, the value of oy and the vol- 
ume content of water near the wall (@) can also be 
determined since n@ adjustment factor to the mixing 
formula is necessary for X-band radiation. By knowing 
oy and @, the value of the S-band adjustment factor to 
the mixing law (e) can be calculated for the near-wall 
region. Since the S-band adjustment factor to the mix- 
ing formula determined at any given time and the di- 
electric constant of water can be assumed to be constant 
throughout the volume of the pipe, these values of p and 
ow can also be applied to the S-band link across the 
pipe. While it is not possible to unambiguously measure 
the phase shift in the S-band link across the pipe, it is 
possible to calculate it if the distance between the S- 
band antennas in the wall, and the two ports is known. 

The wave attenuation in the S-band link across the 
pipe is then measured, making it possible to determine 
the overall volume fraction of water contained within 
the pipe ($) in the region between the antenna ports. 

The apparatus for carrying out the above calculations 
consists of a first antenna port mounted flush in a wall of 
the pipeline, the first antenna port contains an X-band 
link and a sidewall S-band link. The X-band link con- 
sists of an X-band transmitting antenna and an X-band 
receiving antenna located within one-half of the mini- 
mum expected effective wavelength of the X-band radi- 
ation. The sidewall S-band link consists of a sidewall 
S-band transmitting antenna and a sidewall S-band re- 
ceiving antenna located likewise within one-half of 
thewavelength of the S-band radiation. A main link 
S-band transmitting antenna is also included in the first 
antenna port. : 

A second antenna port is mounted flush in the wall of 
the pipeline substantially diametrically opposite the first 
antenna port. The second antenna port contains a main 
link S-band receiving antenna. The device also includes 
a means for generating waves, a means for measuring 
phase shift, and attenuation of a signal transmitted 
through the X-band link, a means for measuring phase 
shift and attenuation of a signal transmitted through said 
sidewall S-band link, and a means for measuring attenu- 
ation of a signal transmitted from the main link S-band 
transmitting antenna to the main link S-band receiving 
antenna. 

Duplicates of the above antennas may be provided to 
increase the accuracy of the device. For example, the 
sidewall X-band link may comprise an X-band transmit- 
ting antenna and two X-band receiving antennas. Simi- 
larly, the S-band links may comprise a single transmit- 
ting and two receiving antennas. Further, the second 
antenna port may be similarly arranged to the first an- 
tenna port. When the first antenna port transmits in the 
S-band link, the second antenna port receives and vice 
versa. It is also possible to provide two sets of the an- 
tenna ports along horizontal and vertical diameters of 
the pipe to detect and quantify gravitational settling. 
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Unlike the prior art, in utilizing the above approach it 
is not necessary to assume that the overall water content 
of the oil in the pipe is represented by an individual 
portion of the pipe since a bulk volume fraction water 
will be found. Further, it is unnecessary to insert a de- 
vice into the flowing stream of oil which would not 
only increase the pressure drop in the pipeline, but can 
also result in disturbances to the flow stream at the point 
of measurement that will result in measurement inaccu- 
racies. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 shows the ratio of power received to power 
transmitted for various frequency signals. 

FIG. 2 illustrates schematically the manner in which 
the described invention is operated. 

FIG. 3 shows the anticipated antenna configuration 
in the preferred embodiment. The antenna ports are 
shown in the pipeline in cross-section and in enlarged 
front views. 


DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 


In the discussion below it is not necessary to assume 
that the distribution of water over a cross section of a 
pipeline is uniform. It is, however, necessary to assume 
that the distribution of water conductivity is uniform 
over the cross section of the pipeline. In almost all cases 
this will be a valid assumption. It is also necessary to 
assume that the adjustment to the mixing formula for 
high frequency radiation is the same near the wall as it 
is throughout the pipe. This will generally be a valid 
assumption if the shape of the water droplets is substan- 
tially the same throughout the pipe because droplet 
shape largely determines p. 

By transmitting electromagnetic waves some distance 
through a fluid medium (flowing or not), it is well 
known that two parameters of the medium may be 
measured: 

a=phase shift of the wave (radians/meter) 

B=attenuation of the wave (neper/meter) 

From these two parameters it is also well known that 
the real and imaginary parts of the complex dielectric 
constant (€*) can be described as follows: 


e* = € —Vie" (1) 
feet 2 (2) 
peo? 
et = 2h (3) 
pe 
where: 


e’=the real part of the complex dielectric constant 
é’=the imaginary part of the complex dielectric 
constant 
p2=permeability of free space 
w=the circular frequency of the wave (=27f where 
f=the frequency of the wave 
i=the imaginary constant (=V —1) 
When the medium is a two-phase mixture of oil and 
water, it is useful to find a “mixing formula” that will 
relate the dielectric constant of oil (€,*), the dielectric 
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the mixture (€*) to the volume fraction of water (6). 
For higher frequency waves (X-band) the mixing law 
can be defined as follows: 


4 


ext =[(1—d)eo*? + bew*4? (4) 
(where the subscript “‘o” indicates oil, the subscript “w” 
indicates water, and the subscript “x” indicates the 
composite dielectric constant for X-band radiation). 

For lower frequency waves it is necessary to provide 
a “correction factor’ (p) to the above equation in order 
to account for the manner in which the oil is dissemi- 
nated in the oil (e.g., fine mist, oval shaped drops, 
streaks, etc.). Although this correction factor could be 
applied in any of several ways to the above equation, in 
the preferred embodiment the mixing law for S-band 
radiation is described as follows: 


és*=[(1-—d)eo*? + ey*P/2]2 (5) 


(where the subscript “‘s” indicates the composite dielec- 
tric constant for S-band radiation). 

It is also known that the dielectric constants for oil 
and water can be described by the following equations: 


€9* =a — ib (6) 
Es = Eq ioy (7) 
Ce te a 
1 + (iwr)!—¥ o 
where: 


a=real part of the complex dielectric constant for oil 

b=imaginary part of the complex dielectric constant 

for oil 

Ow=water conductivity 

€.) =dielectric cqnstant for water at infinitely high 

frequency 

€és=C—DT+ET? 

T=fluid temperature 

C, D, E=Constants (See later discussion) 

7 =a time constant related to the dielectric absorption 

feature for water 

v=a well known constant (see Shawn et al., J. Chem. 

Physics, vol. 67, p. 2257) 

It is desirable to rely on the highest possible fre- 
quency radiation in using the device for the following 
reasons. First, at low frequencies antenna dimensions 
would become impractically large in relation to the size 
of a pipeline. Second, at low frequencies unwanted 
interactions of the radiation with the metallic walls of 
the pipe set in and, third, at higher frequencies the mix- 
ing formula is relatively easily defined, [e.g., by equa- 
tion (5) above]. 

The upper bound of the frequency is limited by the 
power ratio which can be practically measured. With a 
high quality preamplifier, the minimum detectable 
power level at the receiving antenna would be 80 deci- 
bels below one milliwatt (dBm). A realistic value of the 
transmitted power is 10 watts or 40 dBm. Thus, the 
minimum power ratio is: 


40 dBm—(—80 dBm)=120 dBm 


The ratio of power received (P,) to power transmit- 
ted (P,) is: 
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Z=the distance between transmitting and receiving 

points 

T-=the power transmission coefficient at the anten- 

na/fluid interface 

G=the gain of the transmitting or receiving antenna 

(assumed to be identical) 

=the effective wavelength in the fluid 

FIG. 1 illustrates the minimum power ratio versus 
frequency for a typical petroleum pipeline by plotting 
the quantity 10 log (P,/P, as a function of frequency 
when: 

Power loss due to antenna/fluid mismatch= —3 dB 

Antenna gain=0 dB 

Pipeline diameter =24 inches 

€)=2.4—-i 0.1 

oOw=3.5 mho/m 

From this graph it can be found that in a typical 
application, the maximum frequency which can be em- 
ployed is 3 GHz. 

It should be noted, however, that at this frequency it 
is impossible to (1) measure a across the pipe because 
this value can only be measured unambiguously when 
the pipeline diameter is one-half the effective wave- 
length or less, and (2) measure because at this fre- 
quency it is necessary to apply an unknown correction 
factor p to the mixing equation. 

In order to overcome the above difficulties, the fol- 
lowing device/method is employed as schematically 
illustrated in FIG. 2 (not drawn to scale for illustration 
purposes). 

Two S-band links are installed in a pipeline 1. The 
first consists of an S-band transmitting antenna 2 and 
S-band receiving antennas 3 and 7 (the S-band sidewall 
link) and the second consists of an S-band transmitting 


antenna 2 and S-band receiving antenna 4 (the S-band 


main link). The antennas are mounted flush with the 
pipeline wall in nonmetallic, preferably ceramic, win- 
dows. The main link antennas are placed in the pipeline 
such that a wave passed through the main link will pass 
through a representative portion of the two-phase fluid 
contained within it. In addition, an X-band sidewall link 
consisting of an X-band transmitting antenna 5 and an 
X-band receiving antenna 6 is installed flush with the 
pipeline walls. 

Initially, S-band transmitting antenna 2 transmits a 
signal to S-band receiving antenna 3 and X-band trans- 
mitting antenna 5 transmits a signal to X-band receiving 
antenna 6. The phase shift a and attenuation B of both 
of these waves.can be readily and unambiguously deter- 
mined because of the short distance involved and, there- 
fore, the values of €;* and €,* can be calculated from 
equations (1), (2), and (3) above. Values of a, b, €.., C, 
D, E, 7, v, and T can also be readily determined as 
follows: 

C, D, and E are empirical constants which can be 
found in R. P. Wharton et al. doc. of Petroleum Engi- 
neers, Paper 9267 (1981). 

T can be measured continuously in the pipeline with 
a temperature probe; 

€~ can be readily located in standard tables (=4.6); 

a and b can be measured in the laboratory or derived 

from reference tables (since these values will vary 
only slightly with time and will also vary only 
slightly from one crude oil to the next). Note that 
these values are also frequency independent; 

7 can also be derived from standard reference tables; 

and 
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v can be determined from standard reference tables, 
and in any case can be assumed to be very small or 
approaching 0 (greatly simplifying the calcula- 
tions). It is assumed here to be 0.014. 

With knowledge of these variables, as well as the 
frequency of the X-band and S-band waves, it is possible 
to simultaneously solve equations (4) and (5) for ow, p, 
and ¢ by substituting equations (6) and (7) in for €) and 
€y respectively. Note that it is possible to solve these 
two equations simultaneously to obtain 3 values because 
és* and €x* are complex values. 

The S-band transmitting antenna 2 now transmits a 
signal to S-band receiving antenna 4. The phase shift a 
of this signal cannot be determined unambiguously but 
can be calculated. First, a trial value of the phase shift, 
Qo, is calculated by multiplying the phase shift in the 
sidewall link by the ratio: 


ae Distance from antenna 2 to antenna 6 
J Distance from antenna 2 to antenna 3 


We then add 180° of phase as many times as needed to 
the measured phase shift a to bring the result closest to 
Gp. This will give the correct value of the main links 
phase shift (except in the event that the average water 
saturation is drastically different from the near wall 
saturation). 

The attenuation B of this signal can be easily mea- 
sured and, therefore, the value of €;* can be calculated 
across the pipeline with equations (1), (2) and (3). Fur- 
thermore, the value of p for the sidewall links should be 
approximately the same as for the main link for the 
reasons discussed above. Also, as stated above, it is 
assumed that the conductivity of the water is uniform 
throughout the pipeline and is assumed to be constant. 
Therefore, the spatially-averaged water fraction > 
across the pipeline can be determined from equation (5). 

As stated above, FIG. 2 is a simplified drawing in- 
tended to illustrate the operation of the system. The 
actual antenna configuration of the preferred embodi- 
ment is shown in FIG. 3. In the preferred embodiment, 
multiple transmitting and receiving antennas are pro- 
vided in both the sidewall and main links. In addition, 
the first antenna port is duplicated in the second antenna 
port. This allows accuracy improvement through mak- 
ing sidewall measurements near both the ports and aver- 
aging of their values. In addition, if the values differ 
radically, it will be an indication of asymmetric distribu- 
tion in the liquid mixture composition. 

The device consists of a first antenna port 1 and a 
second antenna port 2 mounted flush in the sidewall of 
a pipeline 3 through which a gas-free oil/water mixture 
4 is flowing. The first S-band main link transmitter is a 
two-element array of antennas designated SA; and SA. 
Similarly, the first S-band main link receiver is a two- 
element array of antennas designated SB, and SB2. This 
configuration is used to created narrow, focused beams 
at transmission and to concentrate energy on the receiv- 
ing antennas, which are likewise focused in the direc- 
tion of the transmitting antennas. The power thus re- 
ceived at the second antenna port is compared with the 
power received at the first antenna port with the an- 
tenna designated SA3. By doing so, the effect of mis- 
match at the antenna/fluid interface is eliminated and 
we obtain the absolute power loss for the path length 
that is the difference between the main link path length 
and the path length from the transmitting antennas SA, 
and SA2 to receiving antenna SA3. 
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Redundant measurements can be made by reversing 

the role of the antenna ports, i.e., by using SB; and SB2 

as transmitting antennas and SA; and SA2 as receiving 
antennas, with SB3 used to eliminate mismatch at the 
antenna/fluid interface. 

The S-band sidewall links consists of antenna SA, as 
the transmitter and SA2 and SA3 as the receivers. The 
power from SA, is sequentially received by SA2 and 
SA3, which are spaced at a distance less than one-half of 
the effective wavelength apart. Upon taking the differ- 
ence of the two measurements, a and 8 can be deter- 
mined and hence e¢’ and e” for the sidewall region 
(which in a typical case will be about one inch from the 
wall). Duplicate measurements are made to obtain 
greater accuracy by using SB; as a transmitting antenna 
and SB2 and SB3 as receiving antennas. 

Similarly, the X-band sidewall links consists of an- 
tenna XA1, used as a transmitting antenna, and XA2 and 
XA3 used as receiving antennas. Upon taking the ratio 
of the two measurements, a and 8 can be determined 
and, hence, ¢e’ and e”’. Again, duplicate measurements 
are made for accuracy in the opposite sidewall by using 
XB, as a transmitting antenna and XB? and XB3 as re- 
ceiving antennas. 

All of the above-mentioned attenuation measure- 
ments will need to be corrected for an inherent attenua- 
tion due to well known “geometric spreading” of radia- 
tion (also known as space loss). This correction factor 
can be calculated from theoretical considerations well 
known to one skilled in the art or may be determined for 
all fluids in the system in the laboratory with fluids of 
known dielectric properties such as alcohol and brines. 

From the X-band measurements of e’ and e”, we can 
readily obtain the values of @ and o, for the near wall 
region using equations (1) to (4) and (7). Using these 
two values and the values of ¢’ and €”’ obtained from the 
S-band sidewall link, we can calculate the value of (p) 
from equation (5). Knowledge of these values allows 
calculation of the water fraction across the cross section 
of pipe. 

All of. the antennas discussed above are circularly 
polarized to minimize the effects of stray radiation. 

It should be noted that the above system could be 
fully automated and computerized so as to provide a 
continuous measurement of the water fraction . This 
would also allow all measurements to be taken in a time 
period short enough that the water content does not 
change during the measurement. It should also be noted 
that by calculating the conductivity of the water that 
the salinity of the water can also be continuously moni- 
tored. 

It is to be understood that while the invention has 
been described above with a great deal of particularity, 
that it is manifest that many details of construction and 
the arrangement of the components may be possible 
without departing from the spirit and scope of the in- 
vention. 

For example, the above-described device could be 
used to measure the composition of materials other than 
water and oil. Further, with a simple extension of the 
above-described calculations, the device could be used 
for three or four-phase measurement. Still further, 
while the invention has been described with reference 
to X-band and S-band radiation, other frequencies may 
be optimum for various pipeline sizes and fluid types. 
Therefore, the invention is not to be limited by the 
embodiment set forth herein but is to be limited instead 
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by the attached claims, including the full range of equiv- 
alency to which each is entitled. 

What is claimed is: 

1. Apparatus for analyzing a multi-phase liquid com- 
prising: 

pipeline means for containing the liquid; 

a low frequency transmitting antenna substantially in 
the S-band near a wall of said means for containing 
a liquid; 

a low frequency receiving antenna substantially in the 
S-band located within one-half wavelength of said 
low frequency transmitting antenna; 

means for determining a phase shift of a low fre- 
quency wave; and 

means for determining an amplitude change of the 
low frequency wave. 

2. Apparatus as recited in claim 1 further comprising: 

a complementary high frequency transmitting an- 
tenna substantially in the X-band located substan- 
tially near said low frequency transmitting antenna; 

a first high frequency receiving antenna substantially 
in the X-band located within one-half wavelength 
of said high frequency transmitting antenna; 

a second high frequency receiving antenna substan- 
tially in the X-band located such that a wave trans- 
mitted from said high frequency transmitting an- 
tenna to said second high frequency receiving an- 
tenna passes through a representative portion of 
the liquid; 

means for determining a phase shift and an amplitude 
change of a wave transmitted from said high fre- 
quency transmitting antenna and said first high 
frequency receiving antenna; and 

means for determining an amplitude change of a 
wave transmitted from said high frequency trans- 
mitting antenna to said second high frequency re- 
ceiving antenna. 

3. Apparatus as recited in claim 2 wherein said means 

for containing a liquid is a petroleum pipeline. 

4. Apparatus as recited in claim 3 wherein said anten- 
nas are located in ceramic windows flush with the walls 
of said pipeline. 

5. Apparatus as recited in claim 4 further comprising: 

a duplicate first high frequency receiving antenna 
substantially in the X-band located adjacent said 
first high frequency receiving antenna; 

a duplicate second high frequency receiving antenna 
substantially in the X-band located adjacent said 
second high frequency receiving antenna; and 

a duplicate low frequency receiving antenna substan- 
tially in the X-band located adjacent said low fre- 
quency receiving antenna. 

6. Apparatus for determining the water content of 

crude oil in a pipeline comprising: 

(A) a first antenna port mounted flush in a wall of the 
pipeline further comprising: 

(1) an X-band link comprising: 

(a) an X-band transmitting antenna; 

(b) an X-band receiving antenna located with 
one-half the wavelength of said X-band trans- 
mitting antenna; 

(2) a sidewall S-band link comprising: 

(a) a sidewall S-band transmitting antenna; 

(b) a sidewall S-band receiving antennalocated 
with one-half wavelength of said first S-band 
transmitting antenna; 

(3) a main link S-band transmitting antenna; 
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(B) a second antenna port mounted flush in the wall 
of the pipeline, said second antenna port containing 
a main link S-band receiving antenna; 

(C) means for measuring phase shift and attenuation 
of a signal transmitted through said X-band link; 
(D) means for measuring phase shift and attenuation 
of a signal transmitted through said sidewall S- 

band link; and 

(E) means for measuring attenuation of a signal trans- 
mitted from said main link S-band transmitting 
antenna to said main link S-band receiving antenna. 

7. Apparatus as recited in claim 6 further comprising: 

(a) a duplicate X-band receiving antenna; 

(b) a duplicate sidewall link S-band receiving an- 
tenna; and 

(c) A duplicate main link S-band receiving antenna. 

8. Apparatus as recited in claim 7 further comprising 

a third and a fourth antenna port located at substantially 

right angles to said first and second antenna port. 

9. A method for determining the relative proportions 

of two mixed, dielectrically-dissimilar fluids in a con- 

tainer comprising: 

(a) transmitting a low frequency wave substantially in 
the S-band through said fluids; 

(b) measuring a phase shift and attenuation of the 
wave to obtain the complex dielectric constant of 
the mixed fluids; and 

(c) calculating the relative proportions of the mixed 
fluids. ; 

10. A method for determining the water content of oil 

in a flowing pipeline comprising: 

(a) transmitting a low frequency signal substantially 
in the S-band between a low frequency transmit- 
ting antenna and a low frequency receiving an- 
tenna, said antennas located within one-half wave 
length; 
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(b) transmitting a first high frequency signal substan- 
tially in the S-band between a high frequency trans- 
mitting antenna and a high frequency receiving 
antenna, said antennas located within one-half 
wavelength of each other and in close proximity to 
said low frequency transmitting and receiving an- 
tennas; and 

(c) transmitting a second high frequency signal sub- 
stantially in the X-band between a second high 
frequency transmitting antenna and a second high 
frequency receiving antenna. 

11. A method as recited in claim 10 wherein said 


mixed fluids are water and oil. 


12. The method as recited in claim 11 further com- 


prising the steps of: 


measuring a phase shift and attenuation of said low 
frequency signal; 

measuring a phase shift and attenuation of said first 
high frequency signal; 

measuring an attenuation of said second high fre- 
quency signal; 

calculating a value of an X-band complex dielectric 
constant of fluid along a wall of the container from 
said phase shift and attenuation of said low fre- 
quency signal; 

calculating a value of an S-band complex dielectric 
constant of said fluid along said wall from said 
phase shift and attenuation of said first high fre- 
quency signal; 

calculating a volume fraction of water of said fluid 
along said wall; 

calculating a correction factor to a mixing law; 

calculating a conductivity of said water; 

calculating a value of an S-band complex dielectric 
constant of fluid within said container; and 

calculating a bulk volume fraction water in said con- 


tainer. 
* * * * * 
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[57] ABSTRACT 


A method and apparatus for determining the interface 
between two materials having different dielectric prop- 
erties, especially in determining the thickness of a hy- 
drocarbon layer on ground water, is disclosed. The 
apparatus consists of an air-filled, cased surveillance 
well in which a logging tool is raised and lowered. The 
logging tool contains an electromagnetic receiving an- 
tenna and two transmitting antennas located above and 
below the receiving antenna. A first wave and a second 
wave are transmitted from the transmitting antennas. 
The first and second waves are offset 180° to produce 
“sum” and “delta” signals using a 180° hybrid. Alterna- 
tively, a single wave can be transmitted from a central 
transmitting antenna to receiving antennas located 
above and below the transmitting antenna. The re- 
ceived signals are then phase shifted 180° using a 180° 
hybrid to produce sum and delta signals. As the tool is 
raised and lowered in the well bore, changes in the sum 
and delta signals will be observed when an interface is 
passed. A logging tool for use in petroleum production 
and exploration wells is also disclosed. 
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METHOD AND APPARATUS FOR INTERFACE 
LOCATION DETERMINATION 


FIELD OF THE INVENTION 


This invention is directed generally towards the de- 
termination of the location of hydrocarbons in a forma- 
tion. Specifically, this invention is directed toward the 
location of fluid interfaces in a hydrocarbon bearing 
formation. 


BACKGROUND OF THE INVENTION 


In recent years serious concern for the protection of 
groundwater from various pollutants has risen. It has 
been found that contaminants may enter groundwater 
from various surface locations rendering groundwater 
unsuitable for drinking or other uses. Among the poten- 
tial sources of groundwater pollution are leaking under- 
ground tanks and, specifically, underground hydrocar- 
bon storage tanks, In many cases, gasoline storage tanks 
associated with service stations have been found to leak 
gasoline to the water table. 

Light hydrocarbons (e.g., gasoline) lost from refinery 
tanks, pipelines and service stations gravitates down- 
ward to settle on top of the water table in the shape of 
a plume. In the center of the plume, the hydrocarbon 
thickness could range up to a couple of feet. Recovery 
of spilled gasoline is usually effected by pumping the 
fluid back up to the surface, for example, see U.S. Pat. 
Nos. 4,630,677; 4,497,370; and 4,273,650. 

It is necessary to obtain an estimate of the total quan- 
tity of gasoline remaining in the plume at any given time 
to estimate the volume of hydrocarbon spilled as well as 
to assess the progress of the cleanup. This can be done 
by determining the thickness of the plume in several 
locations over its areal extent. The current technique 
for doing this is to drill an open well through the plume, 
and let the water and the gasoline flow into the well and 
stand to their hydrostatic equilibrium levels. The depths 
of these fluids and thickness of hydrocarbon in the well 
can then be measured by simple methods readily known 
to one skilled in the art (see U.S. Pat. No. 4,213,340). 

This wellbore is not representative of the fluid/solid 
matrix in the formation, however, and the thickness of 
the free-standing column of gasoline in the open well is 
generally greater than the actual formation thickness of 
the spilled fluid [J. Van Dam, “The Migration of Hy- 
drocarbons in a Water-Bearing Stratum”, Schwille, F., 
“Petroleum Contamination of the Subsoil—A Hydro- 
logical Problem”, both in “The Joint Problems of Oil 
and Water Industries, Symposium held at Brighten, 
England (1967), all incorporated by reference herein]. 

No good means of quantifying this discrepancy exist, 
but it is believed that the two numbers could differ bya 
factor of as much as 2 to 4. Using the thickness of the 
hydrocarbon layer in the wellbore to estimate the 
amount of hydrocarbon spilled would cause a greatly 
exaggerated estimate of the amount of gasoline spilled. 

It is, therefore, desirable to devise a method and 
means of accurately determining the actual thickness of 
hydrocarbon in a subsurface formation. 


SUMMARY OF THE INVENTION 


We have developed a system to measure the hydro- 
carbon thickness in the formation itself. The measure- 
ment consists of lowering a UHF electromagnetic trans- 
mission device down an unperforated air-filled well 
drilled to the lowest seasonal water table or below. A 
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non-metallic Well casing of a material such as PVC 
which allows electromagnetic transmission is used. 

The device is raised from the bottom of the well at 
approximately 10 to 20 feet per minute. Changes in the 
UHF transmission are recorded as it passes fluid inter- 
faces (e.g., hydrocarbon/water) in the formation. The 
changes detected are due to the differences in the di- 
electric properties of the fluid at two locations. The 
device is based on the “null balance” concept, i.e., as 
long as the tool is traversing section of the formation of 
a homogeneous fluid type a “delta” channel signal re- 
mains balanced at zero. When the tool traverses a sec- 
tion of the formation where the dielectric characteris- 
tics change (i.e., an interface), this balance is perturbed 
and a finite signal is recorded. A second channel 
(“sum”) indicates the relative magnitude of the dielec- 
tric constant of the fluids in the formation. 

The first embodiment of the tool is generally cylindri- 
cal with the exception of a protruding pad. The curva- 
ture of the pad surface matches that of the inside wall of 
the plastic casing. On the pad are a receiving antenna, 
and two transmitting antennas located equal distances 
above and below the receiving antenna. The three an- 
tennas are generally identical in all respects, and they 
are matched to radiate efficiently through the plastic 
casing. 

When the transmitting antennas are energized 
through a sum port, they transmit radiation having the 
same amplitude and phase. When energized through a 
delta port, the amplitudes are again the same, but the 
phases differ by 180°. Thus, in the sum-mode, the ampli- 
tudes of radiation received at the receiving antenna 
from the two transmitters are added to one another, 
while in the delta-mode, they are subtracted from one 
another. 

The sensing of the dielectric constant occurs predom- 
inantly via the internally reflected paths within the 
casing wall, which are influenced by the dielectric con- 
stant of the formation opposite the antenna. The paths 
that travel into the formation and then reenter the cas- 
ing wall make a relatively small contribution. 

To operate the device, it is raised and lowered in the 
wellbore and the elevation of the tool in the wellbore is 
carefully monitored. The received signal in the sum- 
mode remains at a constant level until a liquid/liquid or 
liquid/vapor interface is passed. It shifts to a different 
magnitude when an interface is crossed corresponding 
to the change in dielectric constant from one zone to the 
other. The delta-mode signal remains at zero (i.e., the 
difference between the transmissions is zero) before 
passing an interface and becomes finite (i.e., the differ- 
ence is non-zero) only during the crossing of a interface 
boundary. The probe thus senses variations in the di- 
electric constant properties of the formation, rather 
than its absolute values. With this information it is possi- 
ble to determine the exact location of liquid/liquid (i.e., 
water/hydrocarbon) and liquid/vapor (i.e., hydrocar- 
bon/air) interfaces and, therefore, the hydrocarbon 
layer thickness in the formation. 

In a second preferred embodiment, the device is de- 
signed to operate on a continuous basis by supplying 2 
receiving antennas with a single transmission antenna 
located between them. Signals are constantly received 
at the receiving antennas and one of the received signals 
is phase shifted 180° to create a delta signal. 

In summary, a method of determining the location of 
an interface of materials having differing dielectric 
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properties is disclosed comprising the steps of transmit- 
ting electromagnetic radiation from a first transmitting 
antenna and a second transmitting antenna in a well- 
bore, said radiation having a common phase and ampli- 
tude, said first transmitting antenna located above said 
second transmitting antenna; receiving a sum signal at a 
receiving antenna located between said first transmit- 
ting antenna and said second transmitting antenna; 
transmitting electromagnetic radiation having a com- 
mon amplitude and phases differing by 180° from said 
first and said second transmitting antennas; receiving a 
delta signal at said receiving antenna, and; monitoring 
said sum and said delta signals. 

Alternatively, the method of determining the location 
of an interface of materials having differing dielectric 
properties comprises the steps of: transmitting electro- 
magnetic radiation from a transmitting antenna; receiv- 
ing a first signal at a first receiving antenna located 
above said transmitting antenna; receiving a second 
signal at a second receiving antenna located below said 
receiving antenna; adding said first and said second 
signals to create a sum signal; phase shifting said first 
and said second signals a total of 180° to produce phase 
shifted signals; adding said phase shifted signals to cre- 
ate a delta signal, and; monitoring said sum and said 
delta signals. 

An apparatus for determining the location of an inter- 
face of materials having differing dielectric properties is 
disclosed comprising: a first radiation transmitting 
means; a second radiation transmitting means, said sec- 
ond radiation transmitting means alternately capable of 
transmitting radiation of the same phase as said first 
radiation transmitting means and capable of transmit- 
ting radiation 180° out of phase from said first radiation 
transmitting means; radiation receiving means; and 
means for monitoring a received signal. 

An alternative apparatus for determining the location 
of an interface of materials having differing dielectric 
properties is disclosed comprising: a radiation transmit- 
ting antenna; a first receiving antenna above said trans- 
mitting antenna; a second radiation receiving antenna 
below said transmitting antenna; means for phase shift- 
ing and combining a signal from said first and said sec- 
ond receiving antennas to produce a delta signal; and 
means for monitoring the delta signal. 

The concept of null balanced well logging will also 
find utility in logging crude oil production or explora- 
tion wells. A third embodiment is described that would 
be utilized in the logging of crude oil and natural gas 
wells. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 is a cross section of a subterranean formation 
showing the devices used to remove hydrocarbons from 
groundwater. 

FIG. 2 is a cross section of a subterranean formation 
illustrating the method of determining the volume of 
hydrocarbons in a subterranean formation. 

FIG. 3 illustrates schematically a prior art well log- 
ging device used in petroleum production wells. 

FIG. 4 is a schematic diagram generally illustrating 
the layout of the present well logging tool. 

FIG. 5 is a generalized RF diagram for the present 
invention. 

FIG. 6 schematically illustrates the method of opera- 
tion of the present invention. 

FIGS. 7 to 13 are idealized (i.e., numerical simula- 
tions) well logs using the present invention for hydro- 
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4 
carbon layer thickness of 0”, 1”, 2”, 3”, 4”, 5”, and 10”, 
respectively. The formation conditions for the corre- 
sponding depths are shown in each Figure on the left 
side of the FIGURE. 

FIG. 14 shows a schematic of the present invention in 
a subsurface formation. 

FIGS. 15a-c are a detailed drawing of the first pre- 
ferred embodiment in side view (FIG. 15a), front view 
(FIG. 155) and cutaway rear view (FIG. 15c). 

FIG. 16 is a detailed RF diagram of the first preferred 
embodiment. 

FIGS. 17a-c are detailed drawing of the second pre- 
ferred embodiment in side view (FIG. 17a), front view 
(FIG. 17b) and cutaway rear view (FIG. 17c). 

FIG. 18 is a detailed RF diagram of the second pre- 
ferred embodiment. 

FIG. 19 is a well log obtained with the second pre- 
ferred embodiment in a laboratory logging well in a 
model formation with a sand matrix. 

FIG. 20 is a well log obtained with the second pre- 
ferred embodiment in a subsurface well in Boise, Id. 

FIG. 21 is a well log obtained with the second pre- 
ferred embodiment in a well in Boise, Idaho with no 
hydrocarbon layer. 

FIGS. 22a-c are detailed drawing of the third pre- 
ferred embodiment inside view (FIG. 23a), front view 
(FIG. 235), and cutaway rear view (FIG. 23c). 

FIG. 23 is an RF diagram for the third preferred 
embodiment. 

FIG. 24 is an idealized (i.e., numerical simulation) 
well log for the third preferred embodiment. 

FIG. 25 is a simulated well log illustrating the use of 
an “omega” channel. 


. DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 


Recovery of spilled hydrdocarbons is usually ef- 
fected by pumping the fluid back up to the surface, as 
shown in FIG. 1. A slotted well liner 7 is inserted into 
a subterranean formation. The formation contains a 
water saturated region 8, and a hydrocarbon saturated 
region 9. A recovery pump 10 pumps the hydrocarbon 
9 to a surface storage facility 11. A water table depres- 
sion pump 12 is also used in most applications. The 
water table depression pump serves to draw down the 
water table below its natural level. This creates a “sink” 
in which hydrocarbon will tend to accumulate. 

In order to determine the effectiveness of this recov- 
ery procedure, it is necessary to obtain an estimate of 
the total quantity of gasoline remaining in the plume at 
any given time. This can be done as shown in FIG. 2 by 
determining the thickness of the plume in several loca- 
tions at a number of surveillance welis 13. The current 
technique for determining the thickness of a hydrocar- 
bon layer in a surveiliance well is to allow the water and 
the gasoline flow into the surveillance well and stand at 
their hydrostatic equilibrium levels. The thickness of 
these fluid layers in the well can then be measured by 
simple methods. However, the thickness of a free-stand- 
ing column of gasoline is known to be substantially 
greater than the actual formation thickness of the spilled 
fluid, causing a greatly exaggerated estimate of the 
amount of hydrocarbon spilled. 

Some conventional logging tools employing electro- 
magnetic or acoustic waves to probe a petroleum pro- 
duction formation are shown in FIG. 3. The wave is 
transmitted into the formation 1 from the borehole 2 
from an antenna 3, and is received at another point 4 
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above or below the transmitter after the wave has tra- 
versed some distance in the formation. The amplitude 
and the phase of the received signal contain information 
(dielectric constant and conductivity in case of electro- 
magnetic waves, velocity and attenuation in case of 
acoustic waves) about the formation. A duplicate re- 
ceiving location 5 is usually employed, thereby obtain- 
ing measurements that eliminate certain unknown fac- 
tors. Additionally by employing a duplicate transmit- 
ting location 6 (above or below the receiver pair), the 
effects of alternate upward and downward transmis- 
sions are averaged to achieve “borehole compensa- 
tion”. Examples of such tools are the Schlumberger TM 
Borehole-compensated BBHC) Sonic Tool, the Elec- 
tromagnetic Propagation Tool (EPT), and the electro- 
magnetic tool called the Deep Propagation Tool 
(DPT). 

The EPT or DPT device is especially useful in locat- 
ing an interface between materials with radically differ- 
ent dielectric properties (e.g., water and gasoline). 
However, the above-described devices have not been 
utilized in locating materials in groundwater and., more 
importantly, for reasons which are illustrated below, 
this type of tool is not useful in locating the interface 
between materials having similar dielectric properties 
(e.g., air and gasoline). Therefore, the prior art devices 
would not be useful in determining the thickness of a 
hydrocarbon layer in groundwater because they mea- 
sure absolute properties of the formation and it would 
not be possible to locate the air/hydrocarbon interface 
as a result of the tools’ insensitivity to small changes in 
dielectric properties. 

In a groundwater formation, assuming that the forma- 
tion is lossless, the sand matrix, the water, and the gaso- 
line can be represented by real dielectric constants €m, 
€w and €, respectively. The following simple mixing 
laws can then be employed: 

(1) For the water zone (“Zone I’) it is assumed that 
the sand matrix is fully saturated with water, so that the 
dielectric constant for this zone (e7) is 

€r=[0.—$)emt + dey? 6) 
where ¢ is porosity (i.e., the void volume fraction of the 
matrix). 

(2) For the contaminant zone (“Zone I”) it is as- 
sumed that the sand matrix is fully saturated with gaso- 
line, so that the dielectric constant for this zone (ez) is 


€=((1—o)emt + eg}? Q) 
(3) For the air zone (‘Zone IIT”) it is assumed that the 
sand matrix is completely dry (air-filled), so the dielec- 
tric constant for this zone (€777) is 
€1={(1—$)em+o4? (3) 
Typical values of €m, éw, and €g are €m=4.6, €y=78, 


€g2.4. For b=30%, we find: ¢€7=17.2, €7=3.9, 
€77= 3.2. 


The problem with prior art logging tools becomes 
immediately obvious from the above. While the con- 
trast between Zones I and II is great, the contrast be- 
tween Zones II and III is marginal. Such a contrast 
would not be noticeable by the conventional electro- 
magnetic (dielectric) logging tools, even if such tools 
were capable of logging through a plastic casing. The 
Schlumberger Electromagnetic Propagation Tool, for 
example, is a padtype tool and is designed to operate 
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when in direct contact with the formation. The anten- 
nas used herein are matched to look through a plastic 
casing into a formation, while EPT antennas are 
matched to look directly into the formation (with the 
size of the slot, the size of the cavity, and the type of 
filling). Without the ability to discriminate between 
Zones II and III (the air and gasoline zone), it would be 
impossible to quantify the thickness of the plume. The 
Schlumberger DPT tool is a centralized tool and uses a 
relatively low frequency transmission device for the 
purpose of deep penetration, but results in poor spatial 
resolution. It is, therefore, not useful for the present 
purpose where vertical resolution to a few inches is 
required. 

FIG. 4 schematically illustrates the well logging tool 
described herein. The probe is a pad 17 on the logging 
tool 18 that is inserted in plastic well casing 19. The 
curvature of the pad surface matches that of the inside 
wall of the plastic casing. On the pad are a receiving 
antenna 22 and two transmitting antennas 21 and 23 
which are located equal distances above and below the 
receiving antenna. 

FIG. 5 is a simplified RF (radio frequency) block 
diagram for the device. A signal from source 24 and 
transfer switch 25 is transmitted to either the sum (Z) or 
the difference (A) channel of 180° hybrid 51. The hybrid 
is a passive four-port device which, if fed through the 5 
port, divides the power into two equal portions whose 
amplitudes are in phase; if fed through the A port, these 
amplitudes are 180° out of phase. Receiving antenna 22 
monitors the transmitted signals. 

Referring now to FIG. 6, it can be assumed that the 
three antennas are identical in all respects, and that they 
are matched to radiate efficiently into the plastic mate- 
rial (having a dielectric constant €,). It can be further 
assumed that the antennas can be treated as point 
sources, and that the propagation of radiation can be 
treated in terms of ray optics. FIG. 6 shows the signal 
paths from antenna 21 to 22, and from antenna 23 to 22 
into and through the casing well 19. The paths through 
the formation (21, B, C, 22 and 23, E, F, 22) suffer larger 
geometric spread and one additional reflection (at the 
plastic/formation interface) as compared to the paths in 
the plastic material (21, A, 22 and 23, D, 22). Therefore, 
the contribution of the former paths can be neglected. 

If x is the coordinate of the receiving antenna mea- 
sured from some reference depth in the well bore (as 
shown in FIG. 6), let €(x) be the dielectric constant of 
the formation at location x and let d be the distance 
between the receiving and transmitting antennas (as 
shown in FIG. 6). When the wall of the casing 19 is thin, 
the formation is in the near field of the antennas, and 
influences their radiation efficiency. This efficiency is 
difficult to quantify. However, we can obtain an ad hoc 
estimate for this efficiency by assuming that the impe- 
dance contrast between the plastic and the formation in 
the near field of the radiating antenna causes the re- 
flected radiation to ga back into the antennas and onto 
the source. This efficiency is thus equal to the ratio of 
the amplitude transmitted into the formation to that of 
the amplitude of the forward-traveling wave at the 
antenna aperture, which is assumed to be unity Thus, 
the efficiency [n(x)] is: 
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4(x) 


This expression assumes that ¢(x) is always greater 
than €) which will always be true for a plastic cased 
well. This efficiency also applies to the receiving an- 
tenna. 

The power received at the receiving antenna for the 
sum and the delta modes respectively can now be writ- 
ten as 


Vecx + d/2) _\ New + 4/2) -Ne& 
W= Cl ene 1 0 ——_—_— eT, + nx)ne — 
Ve x + +\ ep 


Vex + d/2) — \ Ne + 4/2) —N & 
ex + 


where C is a geometric constant which need not be 
determined since the absolute value of = and A wiil not 
be important (i.e., only changes in these values will be 
important). 

Equations (5) and (6) can now be used to generate the 
series of idealized logs discussed below. The logs below 
assume €7=17.2, €77=3.9, €777=3.2 and €p=2.0, d=2", 
the thickness of the plastic= 4”, and the dry sand poros- 
ity is 30%. The probe response was calculated at every 
0.1 inch interval, and then a fivepoint averaging was 
applied to the calculated data in order to partially simu- 
late the effect of a finite antenna aperture. 

In FIG. 7 there is no spilled hydrocarbon layer, i.e., 
this is the case of a “clean” water table. As the tool is 
moved upward through a clean water table and into the 
air-saturated portion, the sum signal goes through a 
twin peak, and then shifts to a constant lower level. The 
delta signal also goes through a twin peak, and then 
drops back to zero. 

The separation distance between antennas is the prin- 
cipal parameter defining the spatial resolution of the 
probe. In these idealized well logs, this distance is 4 
inches. FIGS. 7 to 13 show the evolving features of the 
logs as the thickness of the spilled layer gradually in- 
creases and exceeds the distance between the transmit- 
ting antennas. 

In the sum channel, we see the gradual emergence of 
two twin-peaks, corresponding to the two formation 
zones of fluids with different dielectric properties. 
These two twin-peaks are clearly discernible starting in 
FIG. 11, where the thickness of the spilled layer equals 
2 times. In FIG. 13, we observe three fully developed 
segments of constant signal levels, corresponding to the 
water, the hydrocarbon, and the dry zones. 

For the delta channel we observe, again beginning in 
FIG. 11, the emergence of a distinct peak feature for the 
hydrocarbon layer/dry sand interface. By FIG. 13, the 
three formation zones have been well demarcated. 

Lastly, we note that the probe response is a convolu- 
tion of the array pattern of the antenna arrays 21, 22, 
and 23 with the dielectric profile 2(x). Thus, it may be 
possible to employ a deconvolution procedure to re- 
cover this profile even when the spilled layer thickness 
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is somewhat smaller than the distance between the 
transmitting antennas. 

While the above discussion reeers to the logging tool 
in general terms, prototype tools have been constructed 
in order to refine the design of the logging tool and its 
related equipment. FIG. 14 illustrates the logging tool 
as it is being raised and lowered (preferably raised to 
prevent hang-up) in the well. Wellbore 2 is completed 
with a non-perforated casing 19 which is Schedule 40 
PVC. The bottom of the casing is isolated with cap 26 
(PVC) and the casing is held centrally within the well- 
bore with non-metallic centralizers 27. A cement plug 
102 is preformed around the casing to weight it below 
the water table. In the preferred embodiment the well- 
bore is 12” in diameter and there is a 3.75” annulus 
between the casing and the wellbore. 


(5) 
on& —d/2) — Ne@ — d/2) — N & 
Ba —d/2) + pares ce Ep 


It is important to evenly sand/gravel pack the annu- 
lus of the well in order to achieve good results with the 
invention. The grain size of the gravel (or sand) 28 is 
chosen to match as closely as possible that of the sur- 
rounding formation. In the preferred embodiment for a 
formation in Boise, Id., this is chosen to be between 20 
to 60 mesh. Cement grout seal 29 is installed above the 
sandpack for the purpose of preventing the intrusion of 
surface waters into the gravel pack. 

The logging tool 18 consists of antenna pad 17 
mounted on body 30. Antenna 21, 22 and 23 are con- 
tained within the antenna pad. Pressure arms 31 are 
provided in order to keep the tool stable and pressed 
against the casing in the well. The tool is held in the 
wellbore with cable 32 which is raised and lowered in 
the wellbore with a constant speed winch 33 and tripod 
34. Information from the tool is transmitted to instru- 
mentation and recording equipment 36. This recording 
equipment provides both analog real-time results and 
digitally recorded data. A weight 100 is used to keep the 
cable taut. Pre-installed metallic bands 101 can be used 
as depth markers/calibrators. It may be desirable to use 
separate cables for raising/lowering the tool and for 
transmitting data. 

Greater detail regarding the construction of a first 
embodiment of the logging tool is provided in FIG. 15. 
To provide both high vertical resolution and deep prob- 
ing depths, two banks of antennas are provided. The 
lower bank of antennas 37 is designed to give the best 
vertical resolution (approximately 2 inches) while prob- 
ing a relatively shallow depth (approximately 1 inch). 
The three antennas (38, 39, and 40) are of the cavity- 
backed slot type having an aperture dimension of 
13” <4”. The small dimensions of the aperture necessi- 
tates that the antenna cavity be partially filled with 
deionized water and then filled completely with a sili- 
cone rubber compound (SilasticTM 732 RTV ad- 
hesive/sealant made by Dow Corning). The return loss 
characteristics of the antennas result in an operating 
frequency of 1225 MHz. 
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In order to minimize direct coupling between the 
antennas via surface currents, the antennas are stag- 
gered so as to be arranged in a spiral fashion. The plane 
of polarization of the lower bank of antennas 37 is verti- 
cal. 

. The upper bank of antennas 42 is designed to give 
somewhat deeper probing depth (approximately 2 
inches) at the sacrifice of some resolution. The two 
transmitting antennas, 43 and 44, are also the cavity 
backed, water/Silastic TM filled slot type, but have an 
aperture of 34 inches by 3 inch, and are turned 90° with 
respect to the lower bank antennas. Thus, the plane of 
polarization of these antennas is horizontal. The orthog- 
onality of polarization of the two antenna banks pro- 
vides sufficient isolation between them if they were to 
be operated simultaneously. 

In the upper bank the receiving antenna 45 differs 
from the transmitting antennas in that it is a slot in a 
metal sheet, backed by a cavity filled with deionized 
water and SilasticTM. The “long” dimension of this 
cavity is horizontal, thus permitting the two transmit- 
ting antennas to be in close proximity to each other (53 
inches). 

“Choke” slots 46 and 47 are placed between the trans- 
mitting antennas and the receiving antenna to minimize 
direct coupling. These slots are filled with RT Du- 
roid TM dielectric having a dielectric constant of ap- 
proximately 10. 

The return loss characteristics of the upper bank of 
antennas give rise to a common operating frequency 
1100 MHz. Pressure arm 31 holds the tool in place in the 
well. 

Internally, the logging tool includes RF source 48, 
variable attenuator 49, switches 50, 54 and 55, hybrid 
51, amplifier 56 and receiver 57. The function of these 
devices is discussed in greater detail below. 

FIG. 16 is an RF system block diagram for the tool in 
the first embodiment. The device consists of a 1-2 GHz 
microwave source 48 which transmits a signal through 
a variable attenuator 49. Switch 50 alternately transfers 
the microwave signal to the sum port 52 and the delta 
port 53 of the 180° hybrid, 51. Switches 54 and 55 alter- 
nately transfer the = and A signals to the upper bank 42 
and lower bank 37 transmitting antennas. 

Receiving antennas 39 and 45 send the received signal 
through switch 55 to 20 dB amplifier 56, and on to 
receiver 57. Appropriate cables are phase matched. 

FIGS. 17a, b, and ¢ show a second (preferred) em- 
bodiment of the logging tool. In this embodiment only 
a single bank of antennas consisting of a single transmit- 
ting antenna 60 and two receiving antennas 61 and 62 
are provided. The antennas are identical and no longer 
have a staggered arrangement; their phase centers are 
placed along a vertical line. The decoupling between 
the transmitting and receiving antennas is achieved by 
placing choke slots 59 filled with radio frequency ab- 
_ sorbing material between them. The spacing between 
the transmitting antenna 60 and the receiving antennas 
61 and 62 is 1 inch and the path of propagation of radia- 
tion is reversed, i.e., there is a single transmitting an- 
tenna 60 and two receiving antennas 61 and 62. Three 
amplifiers/receivers 76 are utilized in this embodiment. 

A third channel called the “mismatch” channel mea- 
sures the power reflected from the transmitting antenna 
due to mismatch between the antenna and the forma- 
tion/ casing combination. 

A 15-Volt DC power source 74 and a 20-Volt power 
source 75 supply the electrical energy for the device. 
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An RF source 63, variable attenuator 64, hybrid 68, and 
rollers 80 are also provided. Information is transmitted 
to the surface via multiple conductor cable 77. In the 
second preferred embodiment the device is held in the 
well with pressure wheels 80 which provide a smoother 
tool traverse during logging. Directional coupler 65 is 
used to sample the power reflected back from the trans- 
mitting antenna, so as to provide a MISMATCH signal. 

FIG. 18 provide an RF diagram for the second pre- 
ferred embodiment. As shown, this device operates and 
records data simultaneously on all three channels on a 
continuous basis, removing the need for any switches In 
this embodiment source 63 transmits a signal to variable 
attenuator 64 and through directional coupler 65 to 
transmitting antenna 60. Receiver 66 and amplifier 67 
provide information regarding mismatch between the 
antenna and formation. 

Signals received at the receiving antennas 61 and 62 
are combined in a 180° hybrid 68 to generate the S 
and A signals, which are amplified in amplifiers 69 and 
70 before being detected at the receivers 71 and 72. 
Note that in this embodiment the center antenna is a 
transmitter and the upper and lower antennas are re- 
ceivers. The path of the radiation here is reversed from 
that described earlier, but achieves the same result. 

As would be expected, actual well logs are often less 
than ideal and some degree of experience must be used 
in their interpretation. FIGS. 19 and 20 show well logs 
obtained with the above described tools. 

In general, the primary emphasis in interpreting the 
well logs is on the delta channel. The sum and the mis- 
match channels provide additional information which is 
used to corroborate the delta channel results, especially 
in the event of ambiguities. 

FIG. 19 shows a sample well log for a laboratory 
logging well using the second preferred embodiment. 
An 8-inch hydrocarbon layer was established at the 17 
and 5-inch levels. The hydrocarbon/water interface 
(labeled as 78) in the log occurs at approximately 163” 
(the middle of the two high delta peaks) and the hydro- 
carbon/ air interface (labeled as 79) occurs at 25}” indi- 
cating a very good correlation 

FIG. 20 shows a sample well log obtained with the 
second embodiment from an aquifier in Boise, Id. with 
a significant layer of gasoline on its surface. The forma- 
tion is a coarse sand, cobbly formation. The log indi- 
cates delta peaks at hydrocarbon/water interface at the 
2690.78’ elevation and an air/hydrocarbon interface at 
2691.31’ elevation. The sum decreases from a high level 
to a lower level roughly corresponding to these peaks 
The mismatch channel also changes levels correspond- 
ing to and supporting the delta peaks The very large 
negative peaks on the mismatch channel are used to 
electronically mark elevations during the logging pro- 
cess 

By contrast to the well logging results, a standard 
monitoring well indicated the top of the hydrocarbon 
layer to be at 2690.50 feet and the bottom of the hydro- 
carbon layer to be at 2689.86 feet. Therefore, the hydro- 
carbon thickness with a standard monitoring well 
would be estimated to be 7.7 inches while the logging 
indicates the thickness to be 6.4 inches. This difference 
is substantial and would be significantly greater in for- 
mation with a less coarse formation. 

FIG. 21 is a well log obtained in the same formation 
as that described in FIG. 21 in an area where no hydro- 
carbon is present on the water table. The delta peak for 
the water/air interface is observed at the 2694.56’ eleva- 
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tion. The sum and mismatch shift to different values at 
this same elevation. No other distinctive delta peak is 
observed, and no support for considering any of the 
noise peaks as another delta interface is seen in the sum 
and mismatch channels. 

The hydrocarbon layer detected by the tool in the 
above-described logs is the saturated hydrocarbon layer 
and its associated capillary fringe. (See the Van Dam 
article cited earlier and incorporated by reference). This 
information can be used to estimate the volume of hy- 
drocarbon spilled to groundwater. 

In all of the above well logs, log interpretation can be 
further simplified with curve smoothing techniques 
familiar to one skilled in the art. 

The concept of null-balanced logging would also fin 
utility in logging a typical crude oil production 0 ex- 
ploratory weil. The antenna for an open hole oil or gas 
well logging tool will be generally configured as shown 
above in the second preferred embodiment. However, 
in a crude oil production well it will be necessary to 
operate the device in an open hole environment since it 
will not be practical to utilize plastic casing and/or to 
maintain the hole in an air-filled condition. Therefore, 
the tool will be completely encased, water tight, and all 
internals will need to be designed for the temperature 
encountered in oil bearing strata. The antennas would 
be designed to operate efficiently when in direct contact 
with a lossy formation. 

In addition, since the material surrounding the tool in 
an open hole will frequently be lossy (conductive), the 
sum, delta, and mismatch signals will be complex, i.e., 
they have both amplitude and phase components. 

To account for those differences, phase discrimina- 
tors are used in place of receivers. These devices give 
the phase and amplitude of an unknown signal n relation 
to a standard or reference signal. 

The sum and delta signals may now be written as: 


2*(x) = (7) 
ik*(x+ fd ikem Ld 
GM*(x)[M*x + de + Mx —dje j 
A*) = (8) 
“ike Le a ik fe Yad 
GM*x)[M*x + de — M*x —dje J 
where: 
G = a “geometric spread” factor that can be 
determined empirically or from theoretical 
modeling as is done, for example, in the EPT 
tool. 


*Indicates a complex quantity. 
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-continued 
wave vector for the region between the 
transmitting antenna and the top 
receiving antenna. 
wave vector for the region between the 
transmitting antenna and the bottom 


receiving antenna. 


Osa a function of the dielectric constant 
and the conductivity. 
M*x) = factor arising from mismatch between the 
transmitting antenna and formation 
M*x+da)= factor arising from mismatch between the 
top receiving antenna and formation 
M*x — d= factor arising from mismatch between the 


bottom receiving antenna and formation. 


Information about M*(x) is recorded by the mismatch 
channel when the tool is at location x. Information 
about M*(x—d) has been recorded when the tool was at 
the location (x—d). Similarly, M(x+d) will become 
M*(x) when the tool is moved upward, a distance of d. 
Note that M*(x—d) and M*(x+d) refer to the receiving 
antennas but these quantities can be determined for the 
transmitting antenna since all three antennas are identi- 
cal. Therefore, the quantities M*(x), M*(x—d) and 
M*(x+d) are known. 

Adding and subtracting equations 7 and 8, it is found: 


(9) 
iG +} 

D*(x) + A*(x) = 2IGM*X)M*X + de 

(10) 
Ge +) 

Z*(x) — A*(x) = 2GM*X)M*x — dje 


Solving for k* (= + 4) ana k* (: - +) 


(11) 


« @)t =*(x) + A*(x 
y (= +7 )- id nf IGM*(xMXx + d) 
(12) 
(2 $= a( aaa) 
2 id Gio —d) 


Since all of the quantities on the right side are known 
for the sum, delta, and mismatch channels k*(x+d/2) 
and k*(x—d/2) can be determined. From these the di- 
electric constants and the conductivities for the regions 
between the transmitting antenna and each of the re- 
ceiving antennas can be found using simple relations 
known to one skilled in the art. 

As a tool travels in the well, the coupling between an 
antenna and the formation changes continuously be- 
cause of changing formation conditions, thickness of the 
mudcake layer and/or rugosity of the borehole wall. 
The prior art tools seek to eliminate this unknown quan- 
tity by taking the ratio of the measurements made by 
two adjacent receiving antennas. This assumes that the 
two receiving antennas at a given time have the same 
coupling. This assumption may or may not be true, but 
in any event, we have no way of knowing if it is true or 
not for a given formation. In our tool we actually mea- 
sure this mismatch, and take this into account in a direct 
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13 
fashion. Besides alleviating the above problem, the mis- 
match channel can also tell us about the borehole rugos- 
ity and the mudcake condition (qualitatively). 

FIGS. 22a, b, and ¢ illustrate the tool that would be 
used in the open hole application discussed above. The 
tool comprises of a body 81, antenna pad 82, and pres- 
sure arm 83. The antenna configuration is similar to the 
second preferred embodiment and includes receiving 
antennas 84 and 85 and transmitting antenna 86. Again, 
choke slots 87 are provided. Internally, the device in- 
cludes power sources 88 and 89, and RF source 90. Also 
included in the device are the following (which are 
discussed in greater detail below): variable attenuator 
91, power sampler 92, 3-way power divider 93, amplifi- 
ers 94, and phase discriminators 95, 96, and 97. Hybrid 
98 and directional coupler 99 perform the same function 
as discussed in relation to the first and second preferred 
embodiments. 

FIG. 23 shows an RF block diagram for the third 
preferred embodiment. FIG. 23 is best described by 
following the electromagnetic signal along its path. The 
source 88 emits a continuous-wave (CW) electromag- 
netic signal which passes through a forward directional 
coupler 92. Here, a small portion of the power is sam- 
pled off to a three-way power-divider 93 which will 
serve to provide three reference signals. 

The main signal is then fed to a variable attenuator 91 
and on through a directional coupler 99 to the transmit- 
ting antenna 86. The attenuator 91 serves to continually 


20 


25 


adjust the transmitted power depending on the nature of 30 


the formation encountered, and can be automatically 
controlled by a feedback loop that also samples the 
received power at the receiving antennas 84 and 85. 

Of the power that is incident at the terminal of the 
transmitting antenna 86, a portion is reflected back be- 
cause of mismatch between this antenna and the forma- 
tion opposite it. This reflected power is sampled by the 
directional coupler 99, amplified in the amplifier 94a 
and fed to a phase discriminator 95 which also receives 
a reference signal from the power divider 93 upon being 
amplified in the amplifier 94. The d.c. outputs of the 
phase discriminator 95 contain information about the 
phase and amplitude of the mismatch signal. 

Of the energy that leaves the transmitting antenna 
and enters the formation, portions are received by the 
receiving antennas 84 and 85. These signals are com- 
bined in a 180° hybrid 98. The outputs of this hybrid are 
the sum and the difference of the two input signals. 

The difference signal is amplified in an amplifier 94c 
and fed to a second phase discriminator 96 which also 
receives a reference signal from the power divider 93 
upon being amplified in amplifier 94d. The d.c. outputs 
of the phase discriminator 96 contain information about 
the phase and amplitude of the delta channel signal. 

The sum signal is likewise amplified in the amplifier 
94e and fed to a third phase discriminator 97 which also 
receives a reference signal from the power divider 93 
upon being amplified in amplifier 94. The d.c. outputs 
of the phase discriminator 97 contain information about 
the phase and amplitude of the sum channel signal. In 
FIG. 23 the cables are adjusted to have the correct 
electrical lengths. 

FIG. 24 illustrates a numerically generated well log 
for the third embodiment. This sample well log assumes 


35 


40 


45 


55 


the tool is operating at 1000 MHz and the thickness of 65 


the oil sand is 10 inches. The distance between the re- 
ceiving antennas and transmitting antenna is assumed to 
be 2 inches. The porosity, water saturation (Sp), and oil 


14 
saturation (So) for the three zones is shown on the Fig- 
ure. The water resistivity is assumed to be 0.2 ohm- 
meters. As with the first two embodiments, the delta 
channel (here, especially the phase signal) shows a clear 
shift at each interface. Support for interpreting these 
shifts as an interface is found in the remaining channels. 

FIG. 25 illustrates the use of an additional parameter 
which can aid in the use of the third embodiment. In this 
Figure, the formation and tool are the same as shown in 
FIG. 24. The plotted parameter, however, is “omega” 
(X). This empirical parameter is defined by: X=(Sum 
Channel amplitude) [cos (delta channel phase)]. As 
shown, this parameter shows a sharp delineation at the 
oil/water contact. 

The frequency of operation of any of the above- 
described embodiments should preferably be such that 
the distance between the transmitting and receiving 
antennas is less than one-half the effective wavelength 
of the radiation. For a given tool spacing, this estab- 
lishes an upper limit to the operating frequency. An 
approximate lower limit of the operating frequency is 
obtained from the consideration of resolution distance. 
If we wish to resolve a distance down to the spacing 
from the transmit to receive antenna, then this distance 
must be greater than approximately one-third the effec- 
tive wavelength. In this way we bracket the range of 
operating frequencies. In typical application, this range 
will be approximately 1000-1500 MHz. 

In the groundwater application, the tool should be 
operated near the top of this range, to maximize resolu- 
tion. In an openhole, the tool should be operated near 
the bottom of the range, in order to maximize the depth 
of investigation. Obviously, the operating frequency 
can readily be varied and the optimum frequency for a 
given application can readily be determined by one 
skilled in the art. 

While the above embodiments have been described 
with a great deal of specificity it is to be understood that 
the scope of the invention is not to be limited thereto, 
but instead should be limited only by the appended 
claims, along with the full scope of equivalents thereto. 

What is claimed is: 

1. A method of determining the location of an inter- 
face of dielectrically different material comprising the 
steps of: 

(a) transmitting electromagnetic radiation from a first 
transmitting antenna and a second transmitting 
antenna in a wellbore, the radiation from the first 
transmitting antenna and the second transmitting 
antenna having a common phase and amplitude, 
said first transmitting antenna located above said 
second transmitting antenna; 

(b) receiving a sum signal resultant from the step (a) 
at a receiving antenna located between said first 
transmitting antenna and said second transmitting 
antenna; 

(c) transmitting electromagnetic radiation from the 
first transmitting and the second transmitting an- 
tennas, the transmitted radiation from the first and 
second antennas having a common amplitude and 
phases differing by 180°; 

(d) receiving a delta signal resultant from step (c) at 
said receiving antenna; 

(e) moving said antennas across the interface; and 

(f) monitoring said sum and said delta signals for 
changes whereby the interface can be detected. 
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2. A method of determining the location of an inter- 
face of dielectrically different materials comprising the 
steps of: 

(a) transmitting electromagnetic radiation from a 

transmitting antenna; 

(b) receiving said electromagnetic radiation at a first 
receiving antenna located above said transmitting 
antenna to produce a first signal; 

(c) receiving said electromagnetic radiation at a sec- 
ond receiving antenna located below said receiving 
antenna to produce a second signal; 

(d) adding said first and second signals to create a sum 
signal; 

(e) Phase shifting said first signal from said second 
signal a total 180° to produce phase shifted signals; 

(f) adding said phase shifted signals to create a delta 
signal; 

(g) moving said antennas across the interface; and 

(h) monitoring said sum and said delta signals for 
changes whereby the interface can be located. 

3. The method as recited in claim 1 or 2 wherein said 
dielectrically dissimilar materials are hydrocarbons and 
air. 

4. The method as recited in claim 1 or 2 wherein the 
dielectrically dissimilar materials are a hydrocarbon in a 
groundwater formation and air. 

5. The method as recited in claim 1 or 2 Wherein the 
transmitting and receiving steps are performed on a 
pad-type device. 

6. An electromagnetic logging tool for determining 
the location of a hydrocarbon-bearing formation com- 
prising: 

(a) an uncased weil; 

(b) an enclosed logging tool, said logging tool further 

comprising: 

(i) an electromagnetic transmitting antenna; 

(ii) a first receiving antenna above said transmitting 
antenna; 

(iii) a second receiving antenna below said transmit- 
ting antennas; 

(iv) means for holding said antennas against a wall of 
said well; 

(v) means for phase shifting signals received by said 
first and said second antennas by substantially 180° 
to produce offset signals and means for combining 
said offset signals to produce a delta signal; 

(iv) means for measuring reflected signals due to 
impedance mismatch between the antenna and 
formation to produce a mismatch signal; 

(vii) means for producing reference signals of stable 
phase and amplitude; 

(viii) phase discriminator means for comparing the 
reference signals to the delta signal and mismatch 
signal to determine their phase and amplitude; 

(c) means for providing said logging tool in said well, 
and; ’ 

(d) means for monitoring the phase and amplitude of 
delta signal and the mismatch signal for changes 
whereby properties of the formation can be mea- 
sured. 

7. The method as recited in claim 5 wherein the fre- 
quency of said electromagnetic radiation is substantially 
in the range from 1000 to 1500 MHz. 

8. The method as recited in claim 5 wherein the fre- 
quency of said electromagnetic radiation is substantially 
1500 MHz. 

9. The method as recited in claim 4 further comprises 
the step of locating a hydrocarbon/water interface. 
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10. The method as recited in claim 1 or 2 wherein said 
transmitting and receiving steps are performed adjacent 
to a nonmetallic casing, said plastic casing defining an 
air-filled well. 

11. The method as recited in claim 1 or 2 further 
comprising the step of measuring power reflected from 
the material back into at least one transmitting antenna 
using a directional coupler attached to a power moni- 
tor. 

12. Apparatus for determining the location of an in- 
terface of dielectrically dissimilar materials comprising: 

(a) a first radiation transmitting means; 

(b) a second radiation transmitting means, said second 
radiation transmitting means alternately transmit- 
ting radiation of the same phase as said first radia- 
tion transmitting means and transmitting radiation 
180° out of phase from said first radiation transmit- 
ting means; 

(c) radiation receiving means located between said 
first and said second transmitting means for pro- 
ducing a received signal; 

(d) means for moving the apparatus across the inter- 
face; and 

(e) means for monitoring said received signal for 
changes as the apparatus is moved across the inter- 
face whereby the interface can be located. 

13. Apparatus for determining the location of an in- 

terface of dielectrically dissimilar materials comprising: 

(a) a radiation transmitting antenna; 

(b) a first receiving antenna vertically displaced from 
said transmitting antenna; 

(c) a second radiation receiving antenna vertically 
displaced from said transmitting antenna opposite 
said first receiving antenna; 

(d) means for phase shifting and combining a signal 
from said first and said second receiving antennas 
to produce a delta signal; and 

(e) means for monitoring the delta signal for changes 
whereby the interface can be located. 

14. Apparatus as recited in claim 13 further compris- 
ing means for adding a signal from said first and said 
second receiving antenna prior to phase shifting to pro- 
duce a sum signal and means for monitoring the sum 
signal. 

15. Apparatus as recited in claim 12 or 13 further 
comprising a pad, said pad coinciding with the inside of 
a wellbore and holding said transmitting and receiving 
antennas. : 

16. A method for determining the location of an inter- 
face in a petroleum formation comprising the steps of: 

(a) inserting a logging tool into an uncased well at a 
first location; 

(b) transmitting electromagnetic radiation from a 
transmitting antenna on said logging tool; 

(c) receiving a first signal at a first receiving antenna 
located above said transmitting antenna; 

(d) receiving a second signal at a second receiving 
antenna located below said transmitting antenna; 

(e) phase shifting signals received by said first and 
said second signals substantially 180° to produce 
offset signals; 

((f) combining said offset signals to produce a delta 
signal; 

(g) moving said logging tool to a second location; and 

(h) monitoring a phase of said delta signal for changes 
as said tool is moved whereby said interface loca- 
tion can be determined. 
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17. Apparatus as recited in claim 15 wherein said 
transmitting and said receiving antennas are adapted to 
operate substantially in the frequency range from 1000 
to 1500 MHz. 

18. Apparatus as recited in claim 15 wherein said 
transmitting and receiving antennas are adapted to op- 
erate at a frequency of substantially 1500 MHz. 

19, Apparatus as recited in claim 12 or 13 further 
comprising means to monitor power reflected from said 
transmitting antennas. 

20. Apparatus as recited in claim 12 or 13 further 
comprising an air-filled, plastic cased well, said well 
adapted to receive said antennas. 

21. Apparatus as recited in claim 20 wherein said well 
is drilled substantially into a groundwater formation. 

22. Apparatus as recited in claim 12 or 13 wherein 
said means for moving is a winch connected to said 
apparatus via a cable. 

23. Apparatus as recited in claim 12 or 13 further 
comprising a second bank of antennas, said second bank 
of antennas in closer proximity to each other than said 
first bank of antennas to give increased vertical resolu- 
tion. 

24. Apparatus as recited in claim 12 or 13 wherein 
said antennas are cavity backed slot antennas imped- 
ance-matched to transmit through a nonmetallic casing. 

25. Apparatus as recited in claim 24 wherein said 
cavity back slotted antennas are filled materials selected 
from the group (a) with a silicone rubber compound and 
deionized water and (b) ceramic barium titanate chips 
and epoxy. 

26. Apparatus as recited in claim 12 or 13 wherein 
said antennas are staggered in a spiral fashion. 

27. Apparatus as recited in claim 12 wherein a long 
dimension of said receiving antenna is horizontal. 

28. Apparatus as recited in claim 12 or 13 further 
comprising choke slots between said antennas. 

29. Apparatus as recited in claim 28 wherein said 
choke slots are filled with a dielectric material. 

30. Apparatus as recited in claim 12 or 13 wherein 
said antennas have their phase centers placed along a 
vertical line. 

31. Apparatus as recited in claim 30 further compris- 
ing choke slots between said antennas. 

32. Apparatus for determining the thickness of a hy- 
drocarbon on a water table comprising: 

(a) a surveillance well, said surveillance well being an 
air-filled, nonmetallic cased well inserted into said 
water table; 

(b) a logging tool, said logging tool further compris- 
ing: 

(i) a first antenna bank, said first antenna bank further 
comprising: a transmitting antenna; a first receiving 
antenna; and a second receiving antenna; 

(ii) a source of electromagnetic radiation operably 
connected to said transmitting antenna; 

(iti) a 180° hybrid, said 180° hybrid arranged to offset 
signals received by said first receiving antenna and 
said second receiving antenna by 180° to produce 
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offset signals and to combine said offset signals to 
produce a delta signal; 

(c) means for recording said delta signal operably 
connected to said hybrid; and 

(d) means for moving said logging tool in said well. 

33. Apparatus as recited in claim 32 further compris- 
ing a gravel pack surrounding said cased well. 

34. Apparatus as recited in claim 32 further compris- 
ing a second bank of antennas on said logging tool, said 
second bank of antennas spaced to provide high resolu- 
tion. 

35. Apparatus as recited in claim 32 wherein said first 
receiving antenna and said second receiving antenna 
and said transmitting antenna are a cavity partially filled 
with water and a silicone rubber compound. 

36. Apparatus as recited in claim 32 wherein said first 
bank of antennas is a staggered bank antennas. 

37. Apparatus as recited in claim 32 wherein said 
nonmetallic cased well is a polyvinyl chloride cased 
well. 

38. Apparatus as recited in claim 32 further compris- 
ing choke slots in said first antenna bank between said 
antennas. 

39. Apparatus as recited in claim 38 wherein said 
choke slots are filled with a dielectric material. 

40. Apparatus as recited in claim 32 wherein said first 
antenna bank has its phase centers along a vertical line. 

41. Apparatus as recited in claim 32 wherein said 
antennas are staggered in a signal fashion. 

42. Apparatus as recited in claim 32 further compris- 
ing a winch connected to said logging tool to raise said 
logging tool. 

43. Apparatus as recited in claim 32 further compris- 
ing means to monitor power reflected from material 
surrounding said casing into said transmitting antenna. 

44. Apparatus as recited in claim 32 further compris- 
ing a pad, said pad containing said first antenna bank 
and adapted to substantially coincide to the inside of 
said surveillance well. 

45. A method of determining the thickness of a hy- 
drocarbon layer on a groundwater table comprising the 
steps of: 

(a) installing an air filled nonmetallic well casing ina 
wellbore, said wellbore traversing the hydrocar- 
bon layer on the groundwater; 

(b) transmitting electromagnetic radiation from the 
well with a transmitting antenna; 

(c) receiving the electromagnetic radiation in the well 
with a receiving antenna; 

(d) moving the transmitting antenna and the receiving 
antenna in the well casing traversing through the 
hydrocarbon layer on the groundwater; and 

(e) monitoring radiation received by the receiving 
antenna for phase changes whereby the thickness 
of the hydrocarbon layer can be determined. 

46. Apparatus as recited in claim 44 wherein said 
antennas are adapted to operate in the frequency range 
from 1100 to 1500 MHz. 

47. Apparatus as recited in claim 44 wherein said 
antennas are adapted to operate at a frequency of sub- 
stantially 1500 MHz. : 

* 


* * *& & 


UNITED STATES PATENT AND TRADEMARK OFFICE 
CERTIFICATE OF CORRECTION 


PATENTNO. : 4,831,331 


DATED > May 16, 1989 
INVENTOR(S): Bibhas R, De et al 


It is certified that error appears in the above-identified patent and that said Letters Patent is hereby 
corrected as shown below: 


Claim 2, Col. 15, line 12, "(d) adding said first and second" 
should read --(d) adding said first 
and said second--. 


Claim 2, Col. 15, line 15, "a total 180°" should read --a total of 
160°. 


Claim 5, Col. 15, line 27, "1 or 2 Wherein" should read -~1 or 2 
wherein--. 


Claim 6, Col. 15, line 55, "(c) means for providing said" should 
read -~means for moving said--, 


Claim 16, Col. 16, line 63, "((£)" should read -- (f) --. 
’ 


Claim 41, Col. 18, line 29, "signal fashion" should read 
--spiral fashion--. 


Signed and Sealed this 
Twenty-seventh Day of February, 1990 


Attest: 


JEFFREY M. SAMUELS 


Attesting Officer Acting Commissioner of Patents and Trademarks 


UNITED STATES PATENT AND TRADEMARK OFFICE 
CERTIFICATE OF CORRECTION 


PATENTNO. : 4,831,331 


DATED > May 16, 1989 
INVENTOR(S): Bibhas R, De et al 


It is certified that error appears in the above-identified patent and that said Letters Patent is hereby 
corrected as shown below: 


Claim 2, Col. 15, line 12, "(d) adding said first and second" 
should read -=-(d) adding said first 
and said second--. 


Claim 2, Col. 15, line 15, "a total 180°" should read --a total of 
180°, 


Claim 5, Col. 15, line 27, "1 or 2 Wherein" should read --l or 2 
wherein--. 


Claim 6, Col. 15, line 55, "(c) means for providing said" should 
read --means for moving said--. 


Claim 16, Col. 16, line 63, "((£)" should read -- (f) --. 


Claim 41, Col. 18, line 29, "signal fashion" should read 
--spiral fashion--. 


Signed and Sealed this 
Twenty-seventh Day of February, 1990 


Attest: 


JEFFREY M. SAMUELS 


Attesting Officer Acting Commissioner of Patents and Trademarks 


United States Patent 19 


De et al. 


CCAR TEU 


US005146086A. 
[11] Patent Number: 


[45] Date of Patent: 


5,146,086 
Sep. 8, 1992 


[54] METHOD AND APPARATUS FOR IMAGING 
POROUS MEDIA AND METHOD FOR 
FABRICATING NOVEL OPTICAL 
MATERIALS 

Bibhas R. De, Laguna Beach; 

Michael A. Nelson, Glen Avon; 

Kiran K. Pande, Laguna Hills, all of 

Calif. 


Chevron Research and Technology 
Company, San Francisco, Calif. 


[21] Appl. No.: 672,194 


[75] Inventors: 


[73] Assignee: 


[22] Filed: Mar. 19, 1991 
[ST] Mtl sks cin snidanecetne GOIN 21/35 
[52] U.S. Che veces ... 250/253; 250/256 
[58] Field of Search «0.0... 250/253, 255, 256; 
73/38 
[56] References Cited 
U.S. PATENT DOCUMENTS 
4,782,501 11/1988 Dixon et al. oes 250/253 
4,839,516 6/1989 Freeman et al... 290/253 
4,994,671 2/1991 Safinya et al. ...... ws 250/253 
4,996,421 2/1991 Raiet ab. wc cseeseeee 250/253 


OTHER PUBLICATIONS 


“Recommended Practice for Core-Analysis Proce- 
dure”, American Petroleum Institute, API RP 40, First 
Edition, Aug. 1960, Section 3.4 and 3.5 pp. 18-36. 
Research Methods Papers, Journal of Sedimentary Pe- 
trology, vol. 59, No. 4, Jul. 1989, pp. 613-635. 
“Application of Capillary Pressure Measurements to 
the Determination of Connate Water Saturation”, N. R. 
Morrow and J. C. Melrose, Interfacial Phenomena in 
Petroleum Recovery, pp. 257-669. 


BOISE SANDSTONE 


(SLLYVAMOXNDIN ) YSMOd GALLINSNVYL 


L. C. Uren, “Petroleum Production Engineering”, 
Fourth Edition, pp. 660-669. 

B. F. Swanson, “A Simple Correlation Between Perme- 
abilities and Mercury Capillary Pressures”, Journal of 
Petroleum Technolgy, pp. 2498-2504. 


Primary Examiner—Janice A. Howell 

Assistant Examiner—Richard Hanig 

Attorney, Agent, or Firm—Edward J. Keeling; David J. 
Power 


[57] ABSTRACT 


A method and apparatus is described for determining 
petrophysical properties of porous media by the use of 
light transmission through the media. Specifically, by 
using a highly collimated light source, preferably a 
laser, the permeability, grain size, wettability, porosity 
and clay swelling behavior of a sample can be deter- 
mined on the scale of millimeters. A calibration curve 
or function is first established to correlate light trans- 
mission through a slab of sandstone on the order of 5 to 
10 mm under known conditions. This curved is then 
used to determine unknown conditions by measuring 
light transmission and comparing the values to those on 
the curve. Other applications for the phenomenon of 
light transmission through porous media are disclosed 
whereby conditions in a flowing stream of liquid can be 
monitored by use of light transmission. Specifically, 
changes in the refractive index and immiscible contami- 
nation can be detected. Finally, the phenomenon de- 
scribed can be used to model subterranean hydrocarbon 
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METHOD AND APPARATUS FOR IMAGING 
POROUS MEDIA AND METHOD FOR 
FABRICATING NOVEL OPTICAL MATERIALS 


FIELD OF INVENTION 


The present invention relates to the field of labora- 
tory petrophysics, and more specifically, to using light 
to measure rock and fluid properties within a given 
formation. 


BACKGROUND OF THE INVENTION - 


Determining subterranean reservoir properties is an 
important function for a wide array of applications, 
such as the production of minerals, e.g., oil and gas, 
groundwater technology and modeling contaminant 
transport in the subsurface. These properties include the 
permeability, wettability, porosity and grain size of the 
rocks, as well as the large scale flow patterns in the rock 


10 


formations. The present invention lends itself to all of 20 


these areas. 

The permeability of a material is a measure of the 
ability of a porous medium to transmit fluids through its 
pore spaces. For sufficiently slow, linear, incompressi- 


ble steady flow it is described by Darcy’s law: 
guL 
K= (Par 
where 


K=permeability, darcys (d), 

q=flow rate, cc/sec, 

p.=fluid viscosity, centipoise (cp), 

L=sample length, cm, 

AP=fluid pressure differential across the sample, 

atmospheres (atm), 

A=sample cross-sectional area perpendicular to flow 

direction, cm2. 

The most commonly used permeability units are the 
“darcy” (d) and “millidarcy” (md). A rock having one 
darcy permeability conducts 1 cc/sec volume of 1 cp 
viscosity fluid through a cubic sample having 1 cm 
length sides under a pressure differential across the 
sample of I atm. Normally, permeability is determined 
by taking core samples from the reservoir and carrying 
out well defined measurement techniques, well known 
to those skilled in the art, on the samples. Several such 
techniques available for making these measurements are 
described in PETROLEUM PRODUCTION ENGI- 
NEERING - DEVELOPMENT by L. C. Uren, 
Fourth Edition, McGraw-Hill Book Company, Inc., 
1956 and in American Petroleum Institute, API REC- 
OMMENDED PRACTICE FOR CORE-ANAL- 
YSIS PROCEDURE, API RP 40, 1960. 

More particularly, permeability plays a very impor- 
tant role in describing the fluid flow in oil and gas reser- 
voirs. Two primary methods of measurement are prac- 
ticed in the industry: steady state method and the dy- 
namic displacement method. In each a cylindrical core 
is first saturated with water or brine and then oil 
flooded to irreducible water saturation. The core is then 
water flooded or brine flooded and the pressure drop 
across the core is measured along with the oil or water 
or brine production. The average saturations within the 
core are determined from the overall material balance. 
However, the steady state method requires lengthy 
measurement times since it requires stabilization of the 
fluid flow; while the dynamic displacement method 
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does overcome this, it suffers from capillary end effects, 
and is therefore only effective for high flow rates. 

Both of the above described methods of measuring 
permeability are additionally limited in the information 
they provide. Specifically, both methods provide only 
the average or bulk permeability for the core, and can 
only achieve a resolution on the order of about an inch 
or so, where core diameters are in the range of 1 inch to 
4 inches. Additionally, it is well known that permeabil- 
ity can vary considerably in sandstones over a distance 
of only one or two millimeters. The degree of variation 
in permeability has justifiably received much attention 
in recent years. Information about permeability varia- 
tion within a formation can provide valuable informa- 
tion about the geological characteristics of the forma- 
tion. This information is useful in predicting how an oil 
field will develop, the yields that can be expected, and 
how it might compare or contrast with already known 
and developed fields. Running flooding permeability 
tests, as described above, on very small diameter cores 
is not an acceptable solution to the problem, since the 
permeability can still vary considerably along the 
length of the core. Furthermore, running flooding tests 
on small diameter, small length cores is not practical 
because of hydraulic sealing problems and problems 
introduced by boundary effects. 

The only known method that employs light transmis- 
sion for estimating permeability is by examination of 
thin sections of the core. Thin sections of geologic sam- 
ples less than 60 pm thick are customarily examined 
with petrographic image analysis techniques to study 
the micro-features of the sample. From this analysis, 
geological interpretations of the depositional and post 
depositional process which formed the sample can be 
derived. Moreover, physical qualities such as the poros- 
ity and permeability of the sample can be estimated. The 
permeability of the sample, which depends on an inter- 
connected network of individual pores, can be esti- 
mated from the porosity of the sample by a linear re- 
gression formula or by the empirical Kozeny-Carman 
equation. Since a thin section has no meaningful perme- 
ability itself, the above method gives only an indirect 
and qualitative estimate of the permeability of the par- 
ent volume of rock from which the sample was derived. 

A need thus exists for a method to directly measure 
the permeability and permeability variation in a core 
having a macroscopic thickness (i.e. a thickness of many 
layers of grains), with resolution on the order of milli- 
meters. There have been some efforts to obtain permea- 
bility measurements in the “‘less-than-bulk, greater-than- 
microscopic” range, but most proposed methods are 
time consuming and cumbersome. For example, Chan- 
dler, et al. in their article, “A Mechanical Field Perme- 
ameter for Making Rapid, Non-Destructive Permeabil- 
ity Measurements,” J. Sedimentary Petrology, 59, 613 
(1989), discuss an apparatus (an “air minipermeameter”) 
which achieves a resolution on the order of 2 cms. 
While this is an improvement over the bulk method, it 
can still overlook important variations on the millimeter 
level. Swanson, in his article “A Simple Correlation 
Between Permeabilities and Mercury Capillary Pres- 
sures,” J. Petroleum Tech., 33, 2498 (1981), discusses a 
method for determining the permeability of tiny chips 
of rocks using mercury injection porosimetry. While 
this method may provide the desired resolution, it is 
cumbersome and does not provide for a continuous plot 
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or “map” of the permeability and its variation across the 
core. 

The problem of wettability, of great interest in petro- 
Jeum production problems, concerns knowing when 
two fluids (e.g., oil and water) and present in the rock 
pore space, which fluid adheres to, or “wets”, the pore 
walls or the surface of the rock grains, while the non- 
wetting fluid occupies the voids between the grains or 
the pores. The importance and application of this prob- 
lem is discussed in “Application of Capillary Pressure 
Measurements to the Determination of Connate Water 
Saturation,” by N. R. Morrow and by J. C. Melrose in 
Interfacial Phenomena in Petroleum Recovery (N. R. 
Morrow, ed., Marcel Dekker, N.Y., 1991). At present 
there exists no non-invasive, non-destructive method by 
which one can determine whether oil or water wets the 
surface of the rock grains when both fluids are present. 

While the permeability of a rock contains information 
about the grain size, in petrophysical studies it is of 
fundamental importance to make measurements of the 
grain size itself. This parameter can be determined from 
either microscopic measurements, or by measuring with 
sieves of known grid sizes upon disaggregating the rock 
into the constituent grains. Alternately, apparatus are 
commercially available which can determine the grain 
size of a rock sample by use of a laser. One such appara- 
tus is sold under the name MICROTRACR @® by Leeds 
& Northrup Co. Operation of this apparatus is described 
in detail in “Rapid Analysis of Submicron Particles,” by 
P. E. Plantz and H. N. Frock. This apparatus requires a 
sample of the composite medium to be disaggregated, 
the particles are then blown in a cloud through a laser 
beam, with the resulting scattering effect of the laser 
being correlatable to the particle size. This effect is also 
discussed in “Measurement of size and concentration of 
scattering particles by speckle photography” by Gen- 
celi, et al., J. Opt. Soc. Am., V70, 10, 1212, wherein the 
authors admit, “there are moderate restrictions on the 
number of particles per unit volume.” 

The permeability of rocks may be affected by the 
reaction of clays within the rock. Typically, in their 
natural state within a subterranean formation, clays are 
compacted around the grains of the rock. Once the 
natural source fluid is replaced with another fluid, the 
clays may swell, plugging the channels connecting the 
pores, and lowering permeability. It is desirable to 
know not only the total change in permeability, but also 
the rate at which permeability changes. Current means 
in the art for determining the effect that different fluids 
may have on permeability require flooding the core 
with the various fluids and measuring the bulk permea- 
bility. In this manner, however, data can only be col- 
lected at various times, rather than in a continuum. 
Furthermore, it is difficult to determine, with this 
method, precisely how the permeability may be chang- 
ing as a function of time. 

Previous methods of optically characterizing cores 
have concentrated on reflected light. Herbin, et al. dis- 
close In U.S. Pat. No. 4,852,182 a method of longitudi- 
nally splitting a core and imaging the cross section by 
reflected normal light. Along a similar vein, Pruett, et 
al. in U.S. Pat. No. 4,616,134 disclose a method and 
apparatus for optical surface scanning of longitudinal 
cores. However, neither of these methods actually mea- 
sures properties of the core, other than its variation in 
color. 

Recent advances in imaging technology have been 
applied to the area of geologic samples. For example, 
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Smith in U.S. Pat. No. 4,797,906 teaches a method to 
measure porosity in thin sections by use of X-ray reflec- 
tance to stimulate naturally occurring fluorescence in 
the sample. NMR: (Nuclear Magnetic Resonance) has 
been used in U.S. Pat. No. 4,728,892 wherein Vinegar, 
et al. teach a method for determining, among other 
petrophysical properties, permeability by indirect cal- 
culation from measured porosity. Vinegar et al. state 
that a high resolution (1 mm) can be obtained with 
NMR imaging, however, the method is limited to cores 
4.2 cm or less in diameter. In U.S. Pat. No. 4,868,751 
Dogru, et al. claim a method for determining permeabil- 
ity that combines the classic method of confining the 
core to a pressure sleeve with the technology of CAT 
(Computer Aided Tomography) scanning during the 
flooding. Dogru, et al’s method purportedly eliminates 
capillary effects on permeability in the core, but does 
not indicate that the method gives a high resolution 
image of the variation in permeability within the core. 
None of the methods either teach or suggest the use of 
light transmission through the sample to measure per- 
meability. 

The porous media that are the focus of this invention 
are non-absorptive to light. A good example of an ab- 
sorptive material is a thin film of black plastic; while the 
medium is very thin, no light passes because it is ab- 
sorbed by the medium. In non-absorptive materials, 
light might be expected to transmit through the medium 
and emit from the other side, which it does in many 
cases. However, in some instances some or all of the 
light which is transmitted into the non-absorptive me- 
dium may not be transmitted through the medium due 
to diffusion of the light inside the medium. A good 
example of this latter category is a dense pack of irregu- 
lar quartz crystals not unlike sandstone. 

Light transmission through composite media is a | 
subject that has received much attention recently. Cur- 
rent theories and studies of the phenomenon concen- 
trate on the scattering of laser light within the medium 
in order to study the mechanism of propagation. Typi- 
cally, the medium is a smoke or vapor or particles sus- 
pended in a fluid, gel or solid with particle densities so 
low that normally occurring light can easily pass 
through. (e.g., Ishimaru, Kuga, et al., “Scattering and 
Diffusion of a Beam Wave in Randomly Distributed 
Scatters,” J. Opt. Soc. Am. 73, 131-136, 1983). Other 
studies have considered the phenomenon of light trans- 
mission through a porous bed of spherical glass, as for 
example Gate, “The Determination of Light Absorp- 
tion in Diffusing Materials by a Photon Diffusion 
Model”, J. Phys. D.: Appl. Phys., 4, 1049-1056, 1971 
(Great Britain). In these studies, the porous media have 
been transparent to visible light, and so the propagation 
of laser light through them is not surprising. The con- 
cept of laser light transmission through media normally 
opaque to light was considered and discussed in a paper 
by Anderson, “The Question of Classical Localiza- 
tion—A Theory of White Paint?”, Philosophical Maga- 
zine B, 52, 3, 505-509, 1985. In his paper Anderson 
theorized a slab of thickness W of material containing 
random non-absorptive scatters embedded between two 
non-random propagating media. Stating that if the slab 
is thick compared with the conventionally defined 
mean-free-path 1, coherent radiation will be exponen- 
tially attenuated, Anderson went on to state that most of 
the incident radiation will be scattered diffusely and, 
since it has been postulated that there is no absorption, 
the radiation must all come out on one side or the other 
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of the slab W. Thus there is a clear distinction between 
the transmission of diffuse, incoherent radiation of a 
“localized” system and an ‘“‘extended” system, once the 
system is sufficiently large. Unfortunately, “sufficiently 
large” also means “sufficiently opaque” and the experi- 
mental problem may not be all that simple according to 
Anderson, therefore optimum scattering requires struc- 
tures comparable to the wavelength (i.e., particle diam- 
eter is approximately equal to wavelength). In addition, 
although great complexity may be encountered, the 
work may impinge on a number of highly interesting 
and practical systems, such as porous media. This there- 
fore indicates that light transmission through a non- 
absorptive porous medium such as a macroslab of sand- 
stone, as described herein, is not to be expected at all, 
since in sandstone the wavelength of a laser light is 
much less than the particle diameter, and the macroslab, 
due to its thickness, is sufficiently opaque. 

Others have additionally theorized about light trans- 
mission in rock. In a paper titled, “Reflectance and 
Albedo Differences Between Wet and Dry Surfaces”, 
Twomey, et al., Applied Optics, 25, 3, 431-437, 1986, 
the authors discuss why light is preferably absorbed by 
wet rock over dry rock. The theories developed in this 
paper suggest that forward light scattering into porous 
media is dependent on the type of fluid present in the 
voids. However, the focus of the article is absorption by 
the media rather than transmission by the media, and no 
speculation is made as to light continuing to propagate 
through the media, or that the amount of light transmit- 
ted into (i.e., absorbed by) the media might be corre- 
lated to the overall permeability of the media. 

In petroleum production problems and groundwater 
studies, it is necessary to predict fluid flow characteris- 
tics in a stratified earth, where the permeabilities of the 
various strata are known or given. Such characteristics 
are usually determined from computer simulations. In 
conjunction with such simulations, it is useful to have a 
scale-model laboratory “reservoir” where the predic- 
tions of the simulation can be verified. Thus, what is 
needed in the art is a simple, non-destructive, non-inva- 
sive, effective method for determining various proper- 
ties of a porous medium. Such properties include per- 
meability, grain size, wettability, porosity and clay 
swelling of a porous medium, as well as the absence or 
presence of fluid in the porous medium, type of fluid 
present, and scale modelling of fluid flow patterns in a 
stratified porous medium. We have recently discovered 
that laser light will be transmitted through a macroslab 
of sandstone, on the order of 1 mm to 20 mm thick, 
preferably 5-10 mm thick. A phenomenon not expected 
since such a slab is ordinarily opaque to visible light. We 
have further determined that the amount of light trans- 
mitted through the macroslabs is correlated to the per- 
meability of the slab at the point of transmission. Addi- 
tionally, we have discovered that the laser transmission 
through a consolidated sandstone allows for determina- 
tion of average grain size where the sandstone has a 
high concentration of particles per unit volume; thereby 
allowing a nondestructive measurement of said grain 
size. We have still further discovered that the transmis- 
sion of the laser light through a macroslab containing 
water and white oil is dependent upon which of the two 
fluids wets the grain surface, and allows for a method 
for continuously monitoring the change in permeability 
with respect to time in a macroslab core section flooded 
with various fluids. Finally, it has been found that use of 
the transmission phenomenon allows for prediction of 
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fluid flow characteristics when used in conjunction 
with computer simulation, wherein scale model reser- 
voirs are constructed so that light transmission in said 
models can be observed and measured. 


SUMMARY OF THE INVENTION 


The claimed invention provides a method and appa- 
ratus for imaging properties of porous media by light 
transmission through the media. The invention takes 
advantage of a recently discovered phenomenon that 
laser light will penetrate a macroscopically thick slab of 
porous material which is normally opaque to random 
visible light. Further, the amount of light transmitted 
through the medium is correlated to various properties 
of the medium and depends on the types of fluid present 
in the voids of the medium. Thus, by impinging an 
intense, preferably collimated light on one side of a 
porous medium and measuring the amount of light 
transmitted through the medium, one can determine the 
permeability, the grain size, the porosity, and the wetta- 
bility of the medium; as well as the presence, quantities 
and distributions of certain fluids within the medium. 

In one embodiment, the invention is used to measure 
the permeability of porous media as for example sand- 
stones from geologic formations such as petroleum 
reservoirs. A conversion function relating light trans- 
mission to permeability of the rock for a given fluid is 
first established by control laboratory measurements 
and curve fitting. This function can be used for all sam- 
ples from the same formation (i.e., the same rock type). 
A laser beam is then impinged on a core sample a few 
millimeters thick. The sample is scanned by the laser, 
and a sensor opposite the laser records the light trans- 
mitted through the core. The recorded data can then be 
converted to permeability by the conversion function, 
and be displayed graphically by any of several means 
known in the art to give a topological “map” of the 
permeability. In a variation on this embodiment, a core 
sample containing clays is flooded with various liquids; 
the effect the liquid has on the permeability due to clay 
swelling is measured by light transmission through the 
sample. In this way, the change in permeability as a 
function of time can be observed. This same time depen- 
dence of permeability can be observed where an ini- 
tially undisturbed, saturated core section is subjected to 
dynamic flow. The flow will cause fine intergranual 
particles to migrate into the pore throats, blocking the 
flow. 

In a second embodiment the invention is used to de- 
termine the mean grain diameter of porous media, in- 
cluding naturally occurring media such as sandstones. 
Reference curves are generated either from computer 
simulation or measurement of light diffusion in samples 
of known grain size. The experimental results are then 
compared to the reference curves to determine actual 
grain size. 

In a third embodiment the invention is used to deter- 
mine the wettability of porous media such as sand- 
stones. A reference curve is first established correlating 
light transmission to percent oil and water for both 
water wet and oil wet conditions. A sample saturated 
with a known oil/water concentration can then be 
tested and the results compared to the reference curve 
to determine if the grains are oil wet or water wet or are 
in a state of mixed wettability. In a variation of this 
embodiment the heterogeneity of saturation in a par- 
tially desaturated slab of sandstone can be mapped. 
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In a fourth embodiment the invention is used to deter- 
mine the presence or absence of a fluid in a porous 
medium. The fluid can either be static or flowing, and 
the porous medium can be either naturally occurring or 
artificially fabricated. Further, the presence of one of 
two fluids having distinct refractive indices can be de- 
termined, as well as detecting mixing of miscible and 
immiscible fluids with distinct refractive indices. The 
method is similar to that for measuring permeability, but 
the conversion function is established for light transmis- 
sion versus type of liquid, rather than permeability. 

In a fifth embodiment of the invention, porous media 
are used to fabricate scale models of petroleum reser- 
voirs. A light source is then used to simulate fluid flow, 
and light detectors can be used to predict fluid flow 
patterns within the reservoir. The embodiment can 
generally be used to predict fluid flows in any porous 
medium. 

The primary advantage of the claimed method over 
current methods is that the specific petrophysical prop- 
erty of interest can be measured at a high spatial resolu- 
tion and further, that the measurement of permeability 
does not necessitate hydraulic flows in the sample re- 
quiring high pressure plumbing, seals, etc. 

These and other objects, features and advantages of 
the present invention will become apparent from the 
following detailed description, wherein reference is 
made to the Figures in the accompanying drawings. 


BRIEF DESCRIPTION OF THE DRAWINGS 


In order to facilitate the understanding of this inven- 
tion, reference will now be made to the appended draw- 
ings of the preferred embodiments of the present inven- 
tion. 

FIG. 1 is a plot of data showing the relationship be- 
tween laser light transmission through a sandstone and 
permeability of the sandstone. 

FIG. 2 is a drawing of an apparatus for measuring the 
transmission and attenuation of a laser light through a 
dense porous medium. 

FIG. 3 is a drawing of an apparatus for measuring the 
change in permeability of a core sample. 

FIG. 4 is a plot of simulated data predicting the inten- 
sity distribution of the emergent beam from dense po- 
rous media as a function of the distance from the axis of 
the impinging laser beam and the grain size of the me- 
dia. 

FIG. 5 is a sketch of the intensity of laser light trans- 
mitted through a slab of sandstone as a function of the 
distance from the axis of the laser beam. 

FIG. 6 is a drawing of an apparatus for measuring the 
grain size of a porous medium with a laser. 

FIG. 7 is a sketch showing the dependence of light 
transmission through a porous medium saturated with 
two fluids, on the fluid coating the grains of the me- 
dium. 

FIG. 8 is a drawing of an apparatus for monitoring 
fluid passing through a porous medium. 

FIG. 9 is a plot of experimental data showing the 
dependence of the attenuation of a laser through dense 
porous media on the grain size of the media and the 
index of refraction of the fluid saturating the media. 

FIG. 10 is a plot of simulated data predicting the 
dependence of the attenuation of a laser through dense 
porous media on the grain size of the media and the 
index of refraction of the fluid saturating the media. 

FIG. 11 is a drawing showing the relationship be- 
tween light received by a detector in a model reservoir 
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and the proximity of a moving fluid interface to the 
detector. 

FIG. 12 is a drawing of a scale model of a geological 
reservoir using porous media and light to simulate fluid 
flow. 

FIG. 13 is another drawing of a scale model of a 
geological reservoir using porous media and light to 
simulate fluid flow. 


DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 


The claimed invention makes use of a recent discov- 
ery by the inventors, combining two unrelated disci- 
plines to develop the present invention. In the first disci- 
pline, petrophysics, it is known that the size of the 
grains in a porous medium will affect the permeability, 
i.e., in comparing two porous media, one with large 
grains, one with smaller grains, yet having identical 
porosity (i.e., void volume), the medium with the 
smaller grains will have a lower permeability. This is 
attributed to the higher surface area and consequent 
increased viscous drag on the flowing fluid. In the sec- 
ond discipline, optical physics, it is known that, when 
the wavelength of the source light is very small com- 
pared to the size of the particles of the porous medium, 
the diffusive attenuation, or lessening, of light through 
the medium is a function of the particle size. 

This attenuation is due to the greater number of sur- 
faces and the consequent increase in scattering events, 
and thus diffusion. However, as the theory of light 
propagation is generally applied to porous media made 
of glass spheres, and the forward propagation phenome- 
non, known as forward scattering, is usually explained 
in view of this geometry, it is not obvious to what extent 
this phenomenon would apply to sandstone; a medium 
comprising mineral quartz crystals having irregular 
shapes and sizes generally in the range 10-1000 ppm. The 
grains in sandstone, for example, are transparent, some- 
times having a milkiness or a color to them, and, in free 
form, are better known as common sand. As common 
sand is not observed to be transparent there is no reason 
to expect sandstone to be transparent, and in fact, it is 
not transparent to naturally occurring light and most 
artificial light. Even after we theorize that the attenua- 
tion of light transmission through sandstone should 
have a correlation to permeability, in that both are de- 
pendent on grain size, the known opaqueness of sand- 
stone suggests that the phenomenon would not be ob- 
servable, and therefore would not be of any practical 
value. 

It is known from common experience that when 
grains of one transparent medium are surrounded by a 
fluid with an index of refraction approaching that of the 
medium, the medium tends to “disappear” due to the 
lack of refracted light. For example, glass beads in 
water are less “visible” than glass beads in air; the index 
of refraction of water being nearer the index of the glass 
beads than is the index of air. Therefore, one might 
speculate that since a ‘“‘wet” medium is more light-con- 
ductive than a “dry” medium, that light might propa- 
gate better through wet sand than dry, and the phenom- 
enon could be used to advantage. Yet practical knowl- 
edge that wet sand is not transparent in naturally occur- 
ring and most artificial lights teaches away from pursu- 
ing such a theory. 

In experiments attempting to find some correlation 
between permeability of porous media and light trans- 
mission through porous media, we discovered that a 
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sandstone or sand pack, normally opaque, will in fact 
transmit collimated light when the medium is saturated 
with a clear fluid having an index of refraction “n” 
greater than air and most gases (i.e., n=1.0), and fur- 
ther, that the amount of light transmitted through the 
medium is related to the permeability of the medium at 
that point, i.e., the higher the permeability, the greater 
the quantity of light that will be transmitted. However, 
the effect of attenuation of the light signal relative to the 
permeability decreases as the index of refraction of the 
saturating fluid approaches that of the medium itself; at 
that point where the two indices are equal, the medium 
becomes completely transparent and there is no attenua- 
tion of the signal in the medium. The attenuating effect 
will occur for saturating fluids with refractive indices 
greater than and less than the medium itself, but at some 
point if the difference between the two indices becomes 
too large, the attenuation increases to the point where 
little or no light is transmitted. Although some light 
transmission through dry sandstone has been observed 
and is even predictable in some instances, the present 
invention anticipates that the medium will be saturated 
with a liquid when lasers of relatively low power are 
used. For high power lasers, the invention will work 
without the presence of the saturating liquid. 

The laser is preferably used as the light source be- 
cause of its highly collimated nature, wherein coherent 
beam diameters of less than 1 mm are not uncommon, 
which maximize propagation in the forward direction 
for a given amount of incident optical power. Use of a 
laser is not critical; any directed, collimated beam of 
electromagnetic radiation with sufficient intensity can 
also achieve the effect. In fact, as will be seen in one 
embodiment, a less collimated, isotropic light may be 
desirable in some instances where the matrix material of 
the porous medium allows a certain amount of light to 
pass. Use of the collimated light is preferable, however, 
when the matrix material is naturally occurring, such as 
sandstones, as opposed to synthetic, such as glass beads. 
Collimation of an isotropic light source can be achieved 
by any of several means known in the art, such as lenses 
for focusing or parabolic reflectors. Where a laser is 
used, the polarization state of the laser (i.e., either ran- 
dom or linear) is not important for application in this 
invention. However, the same polarizing state should be 
used throughout where data from one application is to 
be correlated to data from a similar application. The 
wavelength of the light used is not restricted to a spe- 
cific value, and can vary from the extreme infrared (IR) 
to the extreme ultra-violet (UV). However, since differ- 
ent wavelengths give different degrees of transmission 
for different applications, some wavelengths may be 
preferable over others. The choice will depend on the 
individual application. The criteria for selection is to 
choose the wavelength that is the most suitable for a 
particular application. “Suitable” includes such factors 
as intensity of light transmitted, dynamic range of trans- 
mission values, filtering effect of natural coloration in 
the porous medium, and possible grain size range of the 
medium. L. F. Gate has discussed the effect of fre- 
quency on diffusive transmission (The Determination of 
Light Absorption in Diffusing Materials by a Photon 
Diffusion Model, J. Phys. D: Appl. Phys., 4, 1049, 1054, 
1971 [Great Britain]) and provides data indicating that 
light scattering, and hence transmission, through porous 
media is a function of the wavelength of the light as well 
as the natural coloration of the medium. Wavelengths 
that are at least one order of magnitude smaller than the 
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mean grain diameter are preferred, and more preferred 
are wavelengths two to three orders of magnitude 
smaller than the mean grain diameter. In light of these 
criteria, the laser should be an ideal instrument for mea- 
suring small scale permeability variations in porous 
media of macroscopic thickness. In the following dis- 
cussion, the word “laser” will be used to describe the 
light source, but this should not be understood as pre- 
cluding the use of other sources of vectored, electro- 
magnetic radiation which qualify under the above dis- 
cussed selection criteria. 

In the first embodiment of the invention, the laser is 
used to measure permeability of a porous medium. One 
requirement is that the individual grains of the medium 
not be opaque to the light being used. This requirement 
allows the invention to include naturally occurring 
media such as those found in geologic formations, such 
as sandstone, and artificial porous media such as sand or 
glass bead packs. In one example of this embodiment, 
the porous medium is Boise sandstone. The first step is 
to determine the correlation between light transmission 
through the medium and the permeability of the me- 
dium for a given fluid. Note that the fluid saturating the 
medium to facilitate the optical measurement is not 
necessarily the same as the hydraulic fluid for which the 
permeability is being determined. One means of deter- 
mining this correlation is by measuring permeability of 
several homogeneous cores by classical means (flood- 
ing), then measuring the corresponding light transmis- 
sion through a macroslab of the corresponding core, 
plotting light transmission versus permeability, and 
determining, by any of several means known in the art, 
a function which fits the plotted data. This function can 
then be used to convert light transmission readings to 
permeability for additional samples of unknown perme- 
ability. Alternately, the function can be plotted to estab- 
lish a calibration curve; this curve should be valid for 
any core taken from the same formation (i.e., the same 
rock type). This assumption is valid as long as the index 
of refraction and clarity of the formation particles re- 
main, fairly constant, and the permeability variation 
results from grain size variation and not, for instance, a 
clogging of the pores by clays, or migration of fine 
particles clogging the pore throats. These latter affects 
will be discussed below. 

Specifically, several small cores (known here as “cali- 
bration cores”), approximately 1” diameter by 1” long, 
from the same formation are obtained. At least six cali- 
bration cores should be used, and preferably 10. How- 
ever, the ideal number is a number of cores which will 
give statistical significance to the data obtained there- 
from. That is, once the data have been curve fitted, a 
Statistical variance can-be determined to see how well 
the curve fits the data. Depending on the accuracy 
desired for the final readings, if the curve does not 
satisfactorily fit the data, more data should be obtained 
to refine the population used to fit the curve, or to 
determine if a different type of curve should be fitted to 
the data. If the primary interest is in permeability varia- 
tion and not absolute permeability measurements, then 
the curve will not need to fit the data as closely. 

The calibration cores should be homogeneous as to 
permeability, with no variation of permeability across 
the sample. One means to determine homogeneity is by 
microscopic examination of a thin section of the calibra- 
tion core, a method well known in the art. Alternately, 
since it is now known that laser light attenuation 
through a macroslab of a porous medium such as sand- 
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stone is dependent on the same variables as permeabil- 
ity, ie., grain geometry, a macroslab, or a section be- 
tween 1 and 20 mm, and preferably between 5 and 10 
mm, of the calibration core can be obtained after the 
permeability is measured by flooding. A laser beam can 
then be projected through this core section at several 
points and the transmitted light measured at the other 
side of the section. If the transmitted light is not seen to 
vary from location to location across the core, then the 
calibration core is homogenous. However, if the trans- 
mitted light does vary, then either the data from that 
core should be rejected, or several light transmission 
readings should be taken at several different points of 
the section and averaged to obtain a light transmission 
number to correspond with the measured bulk permea- 
bility of the original calibration core. If any core section 
exhibits a color variation, then it should preferably be 
excluded from the calibration data. 

A second consideration in selecting the number and 
type of samples to be used in developing the correlation 
function is that the samples should include permeabili- 
ties near the high and low ends anticipated in the forma- 
tion. This will prevent potential erroneous results from 
extrapolation into regions where there is no data to 
provide support. Ideally the thickness of the sections 
used to establish the calibration curve should be the 
same as the sections whose permeability is to be mea- 
sured. However, this does not preclude the possibility 
that results from a test section of a thickness different 
than the calibration section could not be correlated. For 
’ example, by scaling up or down proportionally, correc- 
tions for thickness might be made. In establishing the 
correlation between light transmission and permeability 
for calibration purposes, even in a “homogeneous” sam- 
ple several light transmission readings should be taken 
and averaged. The reason for this averaging is that on 
some scale even “homogenous” materials become het- 
erogeneous. For example, due to the small diameter of 
the laser beam, it is possible that the laser may impinge 
on a surface irregularity detectable on the beam-diame- 
ter scale (<1 mm) which is undetectable on the desired 
mapping scale (>1 mm). This would give erroneous 
results for purposes of calibration. 

Note that the assumption is made that the porous 
media being investigated are nonreactive to the type of 
fluid present in the pores. It is well known that some 
porous media, especially clays, will react with different 
fluids to change permeability. Thus, a core which might 
have a permeability of 500 millidarcies (md) when 
flooded with a 5% brine solution will only have a per- 
meability of 75 md when flooded with fresh water. To 
test for non-reactivity, one core from the formation 
should be flooded with a variety of fluids in core flood- 
ing permeability tests. If the permeability does not 
change with the fluid beyond that mathematically pre- 
dicted from the fluid properties alone, such as viscosity, 
then the porous medium is nonreactive, and the calibra- 
tion curve will not be reactivity-dependent. The fluids 
used in the reactivity test should include oil (preferably 
a white oil so that light absorption is not a factor), brine 
(a brine similar to the brine found in the formation is 
desirable), fresh water, and air. If the core is found to be 
reactive to different types of fluids, then the calibration 
curve will need be developed using only one type of 
fluid, and should be restricted to measurements of cores 
that are saturated with that particular type of fluid. 

Once the calibration curve has been established to.a 
Statistically satisfactory level, and an accompanying 
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correlation function has been determined by curve fit- 
ting means to the data, permeability of a macroslab can 
be measured by use of a laser for any core within the 
same formation. In the example described, five 1" diam- 
eter by 1” long cores with different, yet homogeneous, 
permeabilities were measured for their air permeabili- 
ties. Laser light was then transmitted through a 5 mm 
macroslab of each core and the light transmission mea- 
sured. The data from these measurements gives the 
correlation curve of FIG. 1. While there is not a smooth 
correlation to the data, it should be noted that only five 
points have been plotted, and with the exception of 
point “B” of FIG. 1, a fair correlation does exist given 
the small data set used. 

Generally, the core is not limited in maximum diame- 
ter, but should probably be no smaller than 1 cm in 
diameter. The method is preferably carried out in a dark 
environment to avoid spurious light which may affect 
the results. The section should preferably be between 1 
and 20 mm and more preferably between 2 and 10 mm 
thick. The section should be mounted with one face 
facing the laser and the other face facing the light detec- 
tor. Although in the preferred embodiment the section 
is mounted with its face at 90% to the direction of the 
laser beam and parallel to the light detector, these an- 
gles are not critical. However, the geometry for the test 
apparatus should be the same as the geometry of the 
calibration apparatus. If the geometry is different, how- 
ever, mathematical corrections could be made utilizing 
the variation in geometry and the trigonometric rela- 
tionship therebetween. The preferred embodiment envi- 
sions one laser and one detector in general alignment 
along a line perpendicular to the face of the section 
when the laser and the detector remain static with re- 
spect to each other. In the preferred embodiment, the 
laser and the detector are static while the core section is 
moved through the laser beam until the desired area to 
be scanned has been completely covered. However, it 
would be possible to use a small bundle of multiple 
detectors generally aligned with the laser, as distin- 
guished from a single detector, with the multiple detec- 
tors coupled to an averaging device. Alternately, in the 
preferred embodiment, the core section could be static 
while the laser and detector(s) are dynamic with respect 
to the section and move across the face of the section to 
be surveyed. In a variation on the preferred embodi- 
ment, a moving laser and a static array of detectors 
(with respect to the section) spanning the area of the 
section to be scanned could be used. In still another 
variation a static laser is used with respect to a moving 
detector array and section. 

With reference to FIG. 2, in one example the work- 
piece (1), a water saturated core section 5.5 mm thick 
and 10 cm in diameter, was scanned with a beam (2) 
from a 10 milliwatt, 632.8 nm wavelength, Helium- 
Neon laser (3) with random polarization and 0.68 mm 
beam diameter. The section was placed in a glass dish 
(4) and immersed in a bath of water (5). Typically the 
fluid would be either that found naturally in the section 
(to preserve its natural properties, as discussed above), 
or a fluid which might be injected into the formation, 
such as fresh water or brine in a water flooding opera- 
tion in an oil field. In the example being discussed, the 
laser was located approximately 50 mm from the first 
side of the core. This distance is not critical, however, 
since the laser beam is highly collimated and does not 
diffuse over long distances in air. If dust or mist were 


present in the environment in which the test was being 
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conducted, then the distance from the laser to the sec- 
tion would be more important. In the example the de- 
tector (6) was located within 2 mm of the other side of 
the section. Generally the detector(s) should be located 
as close as possible to the section since the light emitting 
from the detector side of the section is diffuse due to 
diffusion caused within the section. In this example a 
single static laser was used with a single static detector, 
while the section was dynamic. The detector area used 
was 1 cm X1I cm. The section was scanned in an up and 
down pattern with readings of the light transmission at 
the detector recorded at 1 mm intervals, starting with 
one extreme edge of the section. The section was then 
shifted 1 mm to the opposite edge, and readings re- 
corded as the section was again scanned in an up and 
down pattern. This process continues until a grid has 
been formed scanning the complete section in all direc- 
tions. While grid spacing of 1 mm X1 mm was used in 
this example, larger grids may be desirable if less resolu- 
tion is required. 

Smaller or larger detectors may need to be used as a 
compromise between sensitivity and spatial resolution. 
Additionally, detectors with higher or lower sensitivi- 
ties may be required depending on the amount of light 
transmitted through the section and the capabilities of 
the signal conversion and recording devices used in the 
experiment. Since the amount of light transmitted is not 
only dependent on the properties of the section but also 
on the geometry of the experiment and equipment being 
employed, both of which are at the discretion of the 
experimenter, no absolute detector selection criteria can 
be formulated. 

Typically the detector readings are analog voltages. 
These analog readings can be stored by any of several 
means known in the art, or can be converted to digital 
signals by an A/D converter (7) and stored as a re- 
corded value in any digital data collection device, such 
as a computer (8). In this experiment the data was con- 
verted to a digital signal by means of a digital voltmeter 
and then stored in a Hewlett Packard desktop com- 
puter. 

The recorded values obtained can then be displayed 
graphically to give a “map” of the light transmission 
through the section at any given point. Since light trans- 
mission can be correlated to permeability, this “map” 
will give a representation of the permeability variation 
of the core across the section. Alternately, if a mathe- 
matical curve fitting formula has been determined for 
the correlation curve, the recorded values can be ap- 
plied to this formula, and a “map” of the absolute per- 
meability across the section can be plotted. Alternately, 
the data could be displayed in a tabular form or used for 
further calculations. 

In a variation on the first embodiment, the invention 
is used to measure the change in permeability of a core 
section as a function of time. Such a change can occur 
for instance where clay is present in the core section. As 
discussed earlier, often these clays will swell when satu- 
rated with fluids other than the natural reservoir fluid. 
An apparatus for making such measurements is shown 
in FIG. 3. The workpiece core section (1) is placed in a 
clear, watertight container (11) having fluid inlet (12) 
and fluid outlet (13). A photodetector (6) as shown in 
the first embodiment is placed preferably adjacent to 
the side of the container (11) against which the section 
rests. The laser beam (2) from the laser (3) is impinged 
on the side of the section opposite the detector. The 
laser and the photodetector are preferably approxi- 
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mately aligned. The photodetector measurements are 
recorded either on a chart recorder (14) or by a com- 
puter or data logger so that the readings can be corre- 
lated to a time measurement. Other means of accom- 
plishing this objective of recording coincident time 
values and detector measurements which are known in 
the art can be applied. Initially, in recording these time 
values and detector measurements, the container and 
section are flooded with a fluid F1 and the light trans- 
mission through the section is measured in this base 
state. At time t=0 a second fluid F2 is introduced into 
the container by way of the inlet (12). The second fluid 
displaces the first fluid out of the container through the 
outlet (13). The measurements from the photodetector 
are recorded starting at time t=0 and continuing until a 
steady state is reached or the experiment is terminated. 
The recorded values will be correlatable to the permea- 
bility, as taught in the discussion above for the first 
embodiment, and will indicate how permeability 
changes over time as the clay in the section swells or 
contracts. Although the apparatus disclosed only mea- 
sures light transmission at one point on the section, the 
method herein described can be applied at various 
points on the section. Additionally, by first mapping the 
permeability of the section as described in the first em- 
bodiment, variations in the base state permeability can 
be accounted for when monitoring the change over 
time in permeability at various points on the section 
after flooding with the second fluid. 

The time dependent reduction of permeability can 
result from mechanisms other than clay swelling. Spe- 
cifically, the migration of fine intergranular particles 
under the influence of flowing fluid can cause the parti- 
cles to collect in the pore throats, reducing fluid flow 
and thus permeability. Thus the above variation on the 
first embodiment could be used to measure the reduc- 
tion of permeability as a function of time due to the 
migration of fine particles. An apparatus similar to that 
shown in FIG. 3 can be used with the modification that 
the fluid inlet (12) and outlet (13) are oriented on the top 
and bottom of the chamber (11), respectively. 

In a second embodiment the invention described can 
be used to measure grain size and grain size variation of 
porous media. The configuration to determine grain size 
is identical to the configuration described above for 
permeability mapping. However, for this second em- 
bodiment, smaller detectors on the order of 2 mmX2 
mm should be used. This invention makes use of the 
phenomenon that the laser beam, while a pinpoint on 
the incident side of the macroslab when first impinging 
said macroslab, appears as a luminous disk on the emer- 
gent side, with the maximum intensity at the center of 
the disk. The intensity in this disk will decrease as the 
distance from the axis along which the laser is projected 
increases. The distribution of light across the disk will 
have a typically Gaussian shape. In one example, nu- 
merical simulation of the effect was performed by com- 
puter. A three dimensional numerical simulation of ray 
propagation, using a stochastic approach, was used. In 
the simulation, a hypothetical light ray is traced as it 
propagates through the medium. As the ray is scattered, 
its intensity is reduced by the quantity of light reflected, 
and finally, if the ray emerges at the other side of the 
medium, its intensity is recorded. This is done for a very 
large number of rays, on the order of 1 million. In addi- 
tion, the location at which the ray emerges is recorded. 

The results of the simulation are shown in FIG. 4. 
The results have been visually verified and are consis- 
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tent with the type of spreading pattern predicted by 
Genceli, (Measurement of size and concentration of 
scattering particles by speckle photography, J. Opt. 
Soc. Am., V70, 10, 1212) albeit for much lower particle 
densities, yet the same physical phenomenon. In the 
preferred embodiment correlation curves like that of 
FIG. 5 are first established either by numerical simula- 
tion or by measuring the actual emergent beam profiles 
of rock samples having known grain size. Once the 
correlation curves have been established, rock samples 
of unknown grain size can be tested and their emergent 
beam profiles recorded. These stored profiles of sam- 
ples of unknown grain size can then be compared with 
the profiles of samples of known grain size, and the 
unknown grain size determined by suitable interpola- 
tion means known in the art. 

Specifically, an apparatus for determining grain size 
of a porous composite would resemble that of FIG. 2. 
The similarity in apparatus is due to the fact that the 
grain size mapping embodiment is merely a subset of the 
permeability mapping embodiment. However, for grain 
size mapping the single detector (5) used for permeabil- 
ity mapping is replaced by a fixed array of detectors 
covering the entire area of the section to be scanned. 
Detectors having the smallest possible area are pre- 
ferred, with minimal but discernible spacing between 
detectors preferably 0.1 mm or less. This will give high 
resolution, a desirable feature of this embodiment. De- 
tectors having an area of 2 mmX2 mm are commer- 
cially available, as for example from Silicon Detector 
Corp. of Camarillo, Calif. As discussed in the first em- 
bodiment, various combinations of static and dynamic 
arrangements can be configured between the laser 
beam, the core section and the detectors. However, in 
the preferred embodiment the section and detector 
array are static with respect to a dynamic laser beam. 
Using automated means known in the art, a complete 
system can be created for mapping grain size of a core 
section as a function of the surface spatial coordinates of 
the section. Further, these measurements can be part of, 
and simultaneous with, the permeability mapping mea- 
surements described in the first embodiment. With ref- 
erence to FIG. 6, the workpiece (1), a core section, is 
placed on a static array of detectors (20). The core is 
preferably saturated with a liquid such as water or 
white oil, and can be immersed in a bath of the liquid in 
a glass dish as in the example in the first embodiment. 
The glass dish and liquid are omitted from FIG. 6 for 
simplicity. The individual detectors (10a through 10/n) 
of the array (20) are connected by leads [generally (21)] 
to an analog to digital signal converter (22), which 
sends the converted digital signal to a digital computer 
(23). The outputs of the individual detectors can be read 
individually or addressed more efficiently through a 
multiplexer known in the art. The digital computer (23) 
also controls a servo control mechanism (24) which in 
turn controls servo motors (25) and (26) which position 
the laser (3) by any of several means known in the art, 
such as the rack and pinion positioner (27) shown. The 
computer (23) is programmed to move the laser (3) in 
predefined increments in the horizontal and vertical 
directions (corresponding to x and y coordinates) until 
the entire core section has been scanned by the laser 
beam (2). At each position the laser (3) is held steady 
while the (x,y) coordinate position and the correspond- 
ing outputs from the detector array (20) is recorded by 
the computer (23). After the complete core section has 
been scanned the data recorded by the computer is 
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processed to generate an output, which may be a graph- 
ical display (28). In the preferred embodiment the out- 
put is a qualitative topological “map” of the grain size. 
The amplitudes of the elevations in the map correspond 
to the full width at the half-power point of the emergent 
beam as recorded by the detectors, as demonstrated 
graphically by FIG. 8. Areas with finer grains will have 
wider distributions, and hence higher “peaks” in the 
topological map. Alternately, inverses of the emergent 
beam width could be recorded such that the “peaks” 
would correspond to larger grain sizes. As discussed in 
the above paragraph and similarly for permeability 
mapping in the first embodiment, the light transmission 
values themselves could be mapped to show relative 
variation in grain size as a function of location, or, given 
correlation data developed previously, actual grain 
sizes can be plotted as a function of location allowing 
development of a quantitative topological grain size 
map. 

In a third embodiment of the invention, light trans- 
mission through a porous medium is used to determine 
the wettability of the medium. If a porous medium is 
saturated with a combination of two immiscible fiuids, 
such as white oil and water, one fluid will preferentially 
coat the grains while the other fluid will fill the pore 
voids. The fluid coating the grains is known as the 
“wetting” fluid, and the phenomenon is referred to as 
“wettability”. The affinity of the wetting fluid for the 
grains is a function of several variables, including the 
fluid properties and whether one fluid was present be- 
fore the other. The invention makes use of the discov- 
ery that when the wetting fluid has an index of refrac- 
tion closer to that of the porous medium particles than 
does the pore filling liquid, light transmission through 
the medium will approach the maximum. Alternately, 
when the pore filling fluid has an index of refraction 
closer to that of the porous medium particles than does 
the wetting liquid, light transmission through the me- 
dium will approach the minimum. This is due to the fact 
that at optical junctions (the interface between two 
immiscible substances having different refractive indi- 
ces), the nearer the two refractive indices are to one 
another, the more forward scattering there will be. This 
embodiment can be used to determine wettability char- 
acteristics of a porous matrix such as a sandstone under 
normal reservoir conditions and how this matrix will be 
affected by changes in conditions, such as treatment 
with a surfactant. 

Specifically, we start with a macroslab of, for exam- 
ple, dry sandstone of known weight having reasonably 
homogenous permeability. The porosity of the macros- 
lab, which will be referred to as the “section”, is first 
determined by means known in the art. The section is 
then saturated with a first fluid, for example, white oil, 
then weighed. A laser beam is then projected against 
one side of the section and the emergent beam power 
measured and recorded. The first fluid is then incremen- 
tally replaced with a second fluid, which is immiscible 
with the first fluid and having a refractive index differ- 
ent than the first fluid, for example, water. After each 
incremental addition of the second fluid the section is 
weighed and laser light transmission through the sec- 
tion is recorded. The second fluid is added until the 
weight no longer changes. In the water-replacing-oil 
case described here, this stage is known as the state of 
irreducible oil saturation. From the known porosity of 
the section, the known densities of the two fluids, and 
the recorded weights, the percent of each fluid in the 
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section at any given weight can be determined. The 
data is then recorded and would resemble the upper 
curve (31) in FIG. 7. The process is then repeated by 
adding the first fluid to the section initially fully satu- 
rated with the second fluid and recording the weights 
and light transmission values. When this second set of 
data is plotted it will resemble the lower curve (32) in 
FIG. 7, in the case of oil-replacing-water. The complete 
FIG. 7 is now a calibration curve for wettability of this 
porous medium at normal! conditions for the two fluids. 
Conditions can now be varied, for example, the section 
can be saturated with oil and then replaced with water 
which has been treated with a surfactant. 

In one example of the embodiment, a slab of Boise 
sandstone 5 mm thick was first saturated with white oil 
having a refractive index of n=1.45. The refractive 
index of the sandstone matrix medium is that of quartz, 
n=1.55. The oil in the slab was then replaced with 
water with refractive index n=1.33, following the 
method described above, the laser power transmitted 
through the slab was measured and recorded at various 
states of saturation. Once the slab was fully water satu- 
rated (irreducibly oil saturated), the water was then 
flushed from the slab with the white oil, following the 
stepwise method described above. FIG. 7 shows the 
measured points, A, B, C and D, with interpolated 
upper curve (31) indicating the oil-wet grain condition 
and interpolated lower curve (32) indicating the water- 
wet grain condition. As discussed above, such a set of 
curves could now be used to determine what proportion 
of oil wet grains would become water wet if flushed 
with water containing a surfactant. A known volume of 
the surfactant containing water would be added to an 
oil saturated slab for which wettability curves such as 
FIG. 7 has been developed. The percent saturation “‘s” 
is known, the transmitted power “t” is measured, and 
point “x” is determined. This point indicates the portion 
of grains that have become water wet. If “x” falls closer 
to the upper curve, the sample is more oil wet than 
water wet. Consequently, if “‘x” falls closer to the lower 
curve, the sample is more water wet. The oil previously 
covering the grains is now free and recoverable. Infor- 
mation obtained by this method is useful in predicting 
oil recovery increases by use of surfactants in oil fields. 

In the manner that permeability and grain size can be 
mapped as described above in the first two embodi- 
ments, similarly wettability can be mapped. This map- 
ping would be preferably done when one fluid is replac- 
ing another, wherein several maps would be generated. 
In this manner the wetting characteristics of the porous 
material can be observed. By use of a computer, the 
change in wettability as a function of time could be 
calculated, and the change plotted to show a “wettabil- 
ity potential” map. The above method assumes that 
permeability is homogenous throughout the macroslab. 
However, if it is not, a permeability map can first be 
generated for a macroslab by the method of the first 
embodiment, and then corrections can be made for the 
wettability curves using that information. 

Once the permeability for a core section has been 
determined by the method described in the first embodi- 
ment, and the grain size for the same section has been 
determined by the method described in the second em- 
bodiment, the porosity of the section can be calculated 
from the Kozeny-Carman relation, and the results can 
be mapped to produce a topological map of the porosity 
of the section. Thus the invention provides a complete 
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system for measuring and mapping the petrophysical 
properties of porous media such as sandstones. 

In a fourth embodiment of the invention, laser trans- 
mission through porous media can be used to determine 
the presence and type of fluid in the media pores. This 
combination of media and fluid may or may not be 
transparent to non-laser light. However, the embodi- 


-ment assumes laser light will be used due to the colli- 


mated nature of the laser beam which will result in finer 
resolution of the medium characteristics. Further, this 
embodiment assumes that the medium has homogene- 
ous permeability over the section surveyed. Finally, this 
embodiment further assumes the fluid is transparent to 
the light being used. 

As previously discussed, if a porous section is flooded 
with a non-absorbent (to light) fluid of increasing re- 
fractive indices approaching that of the matrix material, 
light transmission through the medium is seen to in- 
crease. Further, as is well known in the art, as the num- 
ber of optical junctions increases, the scattering in- 
creases, and light transmission decreases. An optical 
junction being a junction between two materials having 
different refractive indices. Finally, light transmission 
decreases as the light absorption of a fluid increases. 
That is, light transmission through a porous medium is 
a function of the index of refraction of the fluid, the 
number of optical junctions in the medium, and the light 
absorbance of the fluid. It is important to note what 
light absorbance is being considered, as some fluids may 
be more transparent to one frequency of light than to 
another. Due to this dependence of the light transmis- 
sion on three independent variables, any method for 
determining fluid properties by use of light transmission 
through a porous medium containing the fluid will need 
to hold two of the variables constant. In this embodi- 
ment, it is assumed that the light absorbance of the fluid 
remains constant while either the index of refraction or 
the number of optical junctions vary. Thus, this embodi- 
ment has two variations. 

Since homogenous permeability of the medium is 
required for this embodiment, an artificially fabricated 
medium is preferably used. However, naturally occur- 
ring media can be used, but their homogeneity should 
be first assured by the method described in the first 
embodiment. One means for fabricating the porous 
medium is to create a pack of transparent beads, prefera- 
bly glass or quartz, and preferably irregularly shaped or 
having multiple surfaces to maximize the scattering 
effect. The pack will be referred to hereinafter as a 
“sand pack”. However, as discussed, artificial “sand” 
can be used as well as naturally occurring sand. The 
choice of size of the grains to be used should be made so 
as to optimize the selection criteria; the choice is a com- 
promise between 1) the stronger attenuating effect of 
smaller grain sizes, suggesting the use of larger grains; 
and 2) the larger sensitivity to fluid refractive index 
change for smaller grains, suggesting use of smaller 
grains. The sand pack can be fused together or simply 
contained in a chamber for the purpose of maintaining 
the structural integrity of the pack. The thickness of the 
sand pack will be a function of the particular applica- 
tion, and light transmission and dynamic range of the 
transmission should be the primary criteria for choosing 
the thickness. The attenuating effect of the grains high- 
lights the advantage of the porous medium approach 
over an approach where the refractive index of the fluid 
alone is monitored. In the former approach, there is a 
much higher sensitivity of the transmitted light to 


5,146,086 


19 


changes in refractive index of the fluid than in the latter 
approach. 

In the first variation, the index of refraction of the 
fluid filling the pores of the medium varies. One method 
by which this can occur is if two miscible fluids having 
different refractive indices are mixed, resulting in a fluid 
with an index of refraction different than either of the 
individual fluids. An example of such a mixture is tolu- 
ene and water. This variation could be used to detect 
contamination of one fluid miscible with another having 
a different index of refraction. An apparatus for accom- 
plishing this is shown in FIG. 8. With reference to FIG. 
8, in this example a porous medium (41), fabricated from 
transparent beads (42), is contained in a chamber (43) 
which is used to conduct the fluid (45) to be monitored 
or a slipstream thereof. The chamber has at least two 
opposing points (44a) and (445) which are transparent. 
The remainder of the chamber is preferably opaque to 
minimize light propagation within the walls of the 
chamber to the detector. The beam of collimated elec- 
tromagnetic radiation (2) from the radiation source (3) 
used to monitor the fluid is projected through the trans- 
parent portion of the chamber to a light detector (6) 
adjacent to the chamber. The output (46) from the de- 
tector can be used for any number of functions. For 
example, it can be connected to an analog output indica- 
tor (47), converted to a digital signal and further pro- 
cessed, or coupled to an electrical switch to activate an 
alarm or other device. As the fluid passes through the 
chamber the amount of laser light which is transmitted 
to the detector will vary as index of refraction of the 
fluid varies. One example of an application would be a 
chemical manufacturing process where contamination 
is possible, as for example in alcohol manufacturing 
where water contamination can occur. A slipstream of 
the production stream could be diverted through a 
chamber as described above and checked for variation 
of the index of refraction with the laser as described. If 
the index of refraction varies beyond a predetermined 
acceptable level, the output from the detector could be 
used to signal an alarm or shut down the process. 

Another application of the variation of this embodi- 
ment would be to use the output signal to calculate the 
percent of contaminating fluid present. If the amount of 
light transmitted through the medium when fully satu- 
rated with each liquid is known (by measurement), and 
then the light transmitted in a mixture of the two liquids 
is measured, mixing laws known in the art can be used 
to calculate the percent of each liquid present in the 
mixture. This variation and its applications and exam- 
ples could find uses in chemical, petrochemical, phar- 
maceutical, food and other manufacturing processes. 

A second method by which the index of refraction of 
the fluid saturating the medium can vary is if the me- 
dium sees different fluids with different refractive indi- 
ces at different times. For example, in a phase separation 
process in which two or more fluids with different den- 
sities and different refractive indices are allowed to 
settle out in a vessel and are then pumped out of the 
vessel, the embodiment could be used to detect the 
presence of the different fluids in a common pump out 
line. The output from the detector cell could then be 
used to switch a directional valve to isolate the different 
fluids into different storage vessels. In one example of 
this variation of the embodiment a non-reactive sand 
pack of uniform permeability was saturated in turn with 
toluene (index of refraction n=1.49), with water 
(n= 1.33), and air (n=1.0, i.e., the sand pack was dry). 
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After each saturation, the laser light attenuation 
through the macroslab was measured, with results 
shown in FIG. 9. As FIG. 9 indicates, the general shape 


-of the curve plotting attenuation as a function of grain 


size of the medium remains essentially the same for the 
three saturating fluids (curves 1, 2 and 3). However the 
experimental results indicate that the attenuation curve 
shifts upward to a less attenuated region (i.e., towards 
less negative values) as the index of refraction “n” of the 
saturating fluid increases. In a medium with constant 
grain size, if the attenuation is a known measured value 
then given a set of curves similar to those of FIG. 9, and 
assuming constant light absorption of the fluid, the 
index of refraction of the fluid saturating the medium 
can be determined. If only a few saturating fluids are 
possible, each with a distinctive index of refraction, 
then the saturating fluid can be identified. FIG. 9 also 
shows that smaller grain size of the medium will result 
in greater sensitivity, i.e., greater differences in attenua- 
tion for two fluids having different refractive indices. In 
the example discussed, the attenuation rather than light 
transmission was the recorded value. Attenuation as a 
recorded value is defined as 


Attenuation = 10.0 log (P2/P;) dB 


where P is the power transmitted through a sand pack 
of one thickness and P2 is the power transmitted 
through a sand pack of the same grain size, but a differ- 
ent thickness, where each sand pack is contained be- 
tween two glass plates. Recording the attenuation in 
this manner cancels the effects on the transmission of 
the glass plates, and thus the effect of the medium itself 
can be isolated. This is desirable if experimental results 
are to be compared with theoretical results. 

The experimental results of FIG. 9 are verified by a 
three dimensional numerical simulation of ray propaga- 
tion, as described previously, using a stochastic ap- 
proach. In the simulation, a hypothetical light ray is 
traced as it propagates through the medium. As the ray 
is scattered, its intensity is reduced by the quantity of 
light reflected, and finally, if the ray emerges to impinge 
on a detector at the other side of the medium, it’s inten- 
sity is recorded. The results of one simulation are shown 
in FIG. 10. The simulation generally verifies the experi- 
mental results that show that attenuation of light in- 
creases (i.e., light transmission decreases) as particle size 
decreases and as index of refraction of the saturating 
fluid decreases. 

In the second variation of the fourth embodiment the 
number of interfaces in the medium varies while the 
other factors affecting light transmission remain con- 
stant. This can occur if two immiscible fluids such as air 
and water or oil and water mix in the porous medium, 
or if two phases of two single substances are present, 
such as steam and liquid water. The mixture will create 
new optical junctions where the two immiscible fluids 
meet, creating additional scattering in the medium. The 
increase in the number of junctions will decrease light 
transmission, and thus the presence of, or the variation 
in, the quantity of an immiscible fluid can be detected as 
a result of this decrease in light transmission. An appara- 
tus similar to FIG. 8 can be used for this variation of the 
embodiment as well. One example of an application for 
this variation is a petroleum refining process where 
water is stripped from white oil. If oil becomes en- 
trained in the water, an undesirable effect, the presence 
could be detected by an apparatus similar to FIG. 8 
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which would be installed in an outlet line from a con- 
denser unit condensing the stripped water. The output 
from the detector could then be used to sound an alarm 
or adjust the process. 

In a fifth embodiment of the invention, a fluid satu- 
rated porous media on a larger than core scale but still 
on the scale of a laboratory model can be used to model 
fluid flow behavior in porous structures. This embodi- 
ment is useful for laboratory scale model simulation of 
complex systems and can be used in conjunction with 
computer simulation programs to verify results and 
predictions. The embodiment has the advantage over 
other flow modeling techniques in that hydraulic flows 
are not required, significantly simplifying the modeling 
apparatus. In its broadest form, the system to be mod- 
eled is reproduced on a scale version, with regions or 
zones of differing permeability created from artificial 
sand or glass bead packs with particle sizes preferably in 
the 100-1000 pm range. The mean diameter of the 
grains of the medium should vary no more than does the 
mean diameter of the reservoir being modeled. One or 
more light sources is then located within the model to 
simulate fluid sources, and light detectors are located 
within the model at points at which fluid flow informa- 
tion is desired. Alternately detectors could be located at 
transparent model boundaries. The signals received by 
the detectors will correspond to the fluid flow at the 
detector locations in an equivalent steady state flow 
model, in view of the previous discussion of the effect of 
particle size on permeability and light transmission. 

An example of an application of this embodiment is 
modeling fluid flow in a petroleum reservoir. The 
model can be constructed to either replicate conditions 
in a known reservoir using information obtained from 
cores taken from the reservoir, or the model can simu- 
late a hypothetical reservoir. Because the transmitted 
light from an isotropic source will be attenuated in the 
model reservoir within a very short distance (due to the 
diffusive attenuation as well as the geometric spreading 
of light), it may be desirable to use glass beads since 
they permit more light transmission than quartz grains. 
The reservoir permeability variation is then reproduced 
as closely as possible by using different layers of differ- 
ent bead sizes. The beads may be mechanically bonded 
to one another if desired by sintering, clear epoxy, or 
other means. In an experiment where steady state flow 
characteristics of the reservoir are the desired parame- 
ter to be modeled, the model should be flooded with a 
liquid that will enhance the light transmission. The fluid 
should be clear and preferably non-absorbent to light 
(i.e., “water white’’). The index of refraction of the fluid 
should be selected based on the criteria previously dis- 
cussed in the first embodiment. The simulation can be 
enhanced by modeling the movement of a fluid inter- 
face within the reservoir. For example, if one fluid is 
used in actual field applications as a “drive” fluid, either 
by natural presence (such as ground water or a gas cap) 
to drive oil toward a production point, this event can be 
recreated in the model by adding means for injecting a 
second fluid, having a different refractive index, into the 
model to “drive” the first fluid with which the model is 
initially flooded. As a result, some of the light traveling 
from the light source to the detector will be reflected at 
the fluid interface, and not arrive at the detector. Upon 
arrival of the fluid interface at the detector, the amount 
of light arriving at the detector will change, and the 
occurrence of the interface arrival can be noted. While 
this variation may not provide information regarding 
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flow quantities, it can be useful in predicting dynamic, 
evolving flow patterns, particularly if several detectors 
are located throughout the model. Additionally, data 
collected before the arrival of the fluid interface at the 
detector can provide valuable information. 

As an illustration, FIG. 11 shows an apparatus in 
Diagram A where a permeable zone generally indicated 
by (51) with of length D being isolated between two 
impermeable barriers (52) and (53). A light source (54) 
is isolated from the permeable zone (51) by a glass tube 
(55). At the opposite side of the zone from the light 
source the zone is isolated by a clear wall (56). A light 
detector (6) is pressed flush against the wall! (56) oppo- 
site the permeable zone (51). The permeable zone is first 
saturated with an initial fluid F1 having index of refrac- 
tion n}. A second fluid F2 with index of refraction n2 is 
added to the zone by means of a very small transparent 
injector such as the nozzle (57) which does not interfere 
with the optical measurement. It is assumed that nj is 
not equal to n2, and that Fl and F2 are immiscible. 
Means to allow fluid F1 to escape as fluid F2 is added is 
provided, such as the outlet tubes (58) and (59). As the 
fluid F2 is injected, the interface (I) is created next to 
the interface (1), and moves towards the detector (6). 
The region ahead of the interface (I) is fully saturated 
by fluid F1, and the combination of fluid F1 and the 
matrix material of the porous medium (51) will be re- 
ferred to as Medium 1. The region behind the interface 
(QD) is filled partly with fluid F2 as it is injected and 
partly with fluid F1 (i.e., some residual fluid F1 remains 
behind interface (I) as fluid F2 is injected). This combi- 
nation of fluids F1 and F2 as well as the matrix material 
of the porous medium (51) will be referred to as Me- 
dium 2. The application consists of recording the power 
“P” of the light received at detector (6) as a function of 
time “t”, commencing at t=0 when the injection of 
fluid F2 begins. The recorded result would resemble 
Diagram B of FIG. 11 in the case where n}>n2, and nj 
and nz are less than the refractive index of the matrix 
medium. From Diagram B we can develop Diagram C, 
which gives the location of the interface (I) as a func- 
tion of time, and thus also gives the velocity of the 
interface. An example where this information is useful is 
in the injection of water into a subsurface hydrocarbon 
bearing reservoir. By knowing the fluid interface veloc- 
ity and position, well performance can be predicted. 
Since P is known (by measurement), the distance “‘d” of 
Diagram C is determined from Equation A: 


P= ToAgGQ@)V al Ts (= [at ie el 


where 
‘o=the optical power incident per unit area of the 
interface (I,), which can be determined from the 
properties of the source light (54); 
Ag=area of the face of the detector (6); 
G(d)=the geometric spread factor for light as a func- 
tion of d=ro/d (1, and d as defined in FIG. 11); 
T'g=the power reflection coefficient at the interface 
dn; 

I'7=the power reflection coefficient at the interface 
@; 

T'y=the power reflection coefficient at the interface 
(2); 


1;=the photon mean free path in Medium 1; 
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I,2=the photon mean free path in Medium 2. 

Equation A is solved for d by the following method: 
Since the relationship between P and t is known graphi- 
cally (FIG. 11B) a value for d is selected, P is calculated 
from Equation A, and time t from FIG. 11B that corre- 
sponds to the calculated P is determined, a plot of d and 
tis then made. Successive iterations of this process with 
different values of d, then interpolating a curve between 
the plotted points, will yield the d versus t curve. 

As an illustration of the application of Equation A, 
with reference to Diagram B of FIG. 11, if fluid F1 is 
white oil having index of refraction n; and fluid F2 is 
water having index of refraction n2, at time t=O the 
entire zone (D) of Diagram A is oil-saturated, and the 
detector (6) reads a power level P). In this illustration, 
it is assumed that nj>np, but this is not a requirement, 
and the only restriction is that ny; is not equal to nz, and 
n, and nz are less than the index of refraction (n) of the 
matrix material. At t=0, ['g=T',=I1, where Yr is the 
power reflection coefficient for Medium 1. The value of 
TY) can be determined from a separate reflectance mea- 
surement method known in the art. As soon as the first 
small volume of water is injected by the injector (57) a 
very thin layer of Medium 2 forms around the glass tube 
(55). At that time the power P will suddenly change to 
P2 of Diagram B due to the replacement of I’; with I, 
the power reflection coefficient for Medium 2, and the 
introduction of interface (I), and thus the term I; in 
Equation A. Following this initial event the power P 
then drops more gradually as shown in Diagram B as 
fluid F2 is added and fluid F1 is displaced. When the 
interface (I) gets very near the interface (Ip) (i.e., dD), 
the power will be at P3 of Diagram B. At the point 
where Medium 2 bridges the gap to the interface (12) at 
the glass wall (56), the power will jump to P4 of Dia- 
gram B. This jump is explained by I's, which prior to 
the bridging was equal to T'1, now becomes equal to 2, 
and the elimination of interface (I), and thus the term I’;. 
Following the bridging event there will be no further 
change in P, since at that time an optically static system 
exists. In light of this discussion and diagrams A and B 
of FIG. 11, the solution to Equation A can now be 
solved for d since we know that 
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-continued 


P4 
h=Din irs 


where 

K=T, Ag G(D), a constant for a given apparatus. 

Additional helpful information can be obtained by 
measuring the volume of the effluent of fluid F1 from 
the apparatus as a function of time. 

The above example is only one illustration of how 
fluid flow modeling can be conducted utilizing the pres- 
ent invention. Further refinements of the analysis is 
possible by one skilled in the art. Starting from this 
example, more complex models having different geome- 
tries can be developed. _ 

In modeling the reservoir, the zones should be con- 
toured to replicate undulations or discontinuities in the 
reservoir itself. FIG. 12 shows in cross section one 
example of a model fabricated to simulate fluid flow in 
a petroleum reservoir. A clear glass vessel (61) approxi- 
mately 15” in diameter and having straight sides ap- 
proximately 12” high was used to simulate the reser- 
voir. The bottom of the vessel was internally lined with 
black epoxy (62) to act as a light seal. A tempered clear 
glass tube (63) was sealed to the bottom of the vessel 
with caulk (64). The top portion of the tube protruding 
above the uppermost layer of sand is coated with black 
tape or epoxy (65) to act as a light seal. A tubing (66) of 
3” diameter plastic is run down the inside of the vessel 
(61) to provide a means for emptying fluid from or 
adding fluid to the vessel. A first layer of relatively 
impermeable sand (67) with a mesh size between 140 
and 270 was deposited in the vessel. Further layers of 
differing grain size were deposited to test certain theo- 
ries in experiment, as explained below. The complete 
reservoir model was filled with a fluid (68) (water) 
before the tests were run to improve light transmission. 
A highly compact, high intensity white light should 
preferably be used as the light source. In the example 
being discussed, a Sylvania 750 W tungsten halogen 
lamp (69), suspended from an electrical supply cord 
(70), was used as the source light. Because of the heat 
generated by the intense heat a cooling means such as a 
fan or blower may need to be employed to keep the 
light bulb from overheating and prevent possible dam- 
age to the reservoir model. Alternately, the light can be 
activated for only short periods of time if heat is a prob- 
lem. In the example, short lighting periods were used. 
Although the results were only visually noted in the 
experiments, detectors (71) through (81) could have 
been placed adjacent to the vessel at the different zones 
and measurements taken. A top zone of relatively im- 
permeable fine sand (82) with a mesh size between 60 
and 80 is used as a light seal. 

In a first experiment, a layer of relatively permeable 
glass beads was used to test whether diffusively propa- 
gating light would be confined to high permeability 
channels and not spread out along the vectors 
(100a-100¢), and further, whether it would follow the 
path shown by the vector (101) and not be impeded by 
the contour (83) as it would if it merely followed the 
path shown by vector (102). One side of this layer (84) 
is left flat and the other side (85) is contoured. The bead 
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diameters for this layer ranged between 1.0 and 1.05 
mm. Results demonstrated that light did follow the path 
of vector (101). Although the light at the vessel adja- 
cent to section (85) was not as intense as the light adja- 
cent to section (84), this can be explained by the greater 
path length that light must travel in following vector 
(102) over vector (101) and the resultant increase in 
diffusion. 

In a second experiment utilizing the apparatus of 
FIG. 12, the previous experiment was isolated by a 
region of relatively impermeable glass beads (86a) and 
(86) having bead diameters between 0.10 and 0.11 mm. 
Permeable zones of glass beads (87a), (88a), (88b) and 
(875) were deposited along with an impermeable annu- 
lus barrier (shown in cross section by (89a) and (895)). 
Bead diameters for zones (87a,b) were between 0.45 and 
0.50 mm. Bead diameters for zones (88a,b) were be- 
tween 1.00 and 1.05 mm. Bead diameters for zones 
(89a,b) were between 0.10 and 0.11 mm. The purpose of 
the experiment was to test whether light would follow 
expected fluid flow paths and move along vector (103) 
into zones (87a) and (876) and not flow into zones (89a) 
and (895). Light propagation along vector (104) is gen- 
erally anticipated. The visual results show that light will 
indeed “flow” through the permeable zones in a manner 
simulating the flow of fluids. Specifically, a bright light 
was observed at the vessel wall adjacent to zones (88a) 
and (885), a less bright light was noted adjacent to zones 
(87a) and (876), and virtually no light adjacent to zones 
(89a) and (89b). The difference in light between zones 
(87a,b) and (88a,b) can be explained by the somewhat 
lower permeability of zones (87a) and (87b) compared 
to zones (88a) and (885) and also the longer path length 
the light needs to take in following vector (103) and the 
resulting loss due to diffusion. 

The above results demonstrate the similarity between 
diffusive propagation of light and fluid flow which 
forms the basis of the fifth embodiment of the invention. 
In the dynamic model illustration discussed above 
where one fluid replaces another, this similarity extends 
further to include the similarity between the fluid front 
itself and a “light front” at the interface between the 
two fluids that moves along with the fluid front. 

In fabricating reservoir models, source points simu- 
lating either natural reservoir resources or artificial 
source injection points can be included in the model by 
placing glass tubes at the desired source locations. That 
is, more than one source light may be used. In fact, 
when simulating natural in-situ resources, several 
source points should be used to simulate a broader dis- 
tribution of fluid source. When an artificial source such 
as water or steam injection is to be simulated, the num- 
ber of actual source points should be replicated in the 
model. Production points simulating points at which 
fluid is to be extracted from the reservoir can be simu- 
lated in like manner by placing glass tubes at the desired 
locations within the model and placing detectors inside 
these tubes. Production points are not necessarily points 
at which fluid is to be extracted, and may merely be 
points at which fluid flux is desired to be known. Addi- 
tionally, a source point can act as a production point, 
and vice versa. Within a glass tube, a light source or a 
detector can be moved up or down to simulate injection 
or production at various depths within the reservoir. 
The signal from the detectors can be converted (analog 
to digital) and recorded. The results can be used for 
such purposes as estimating fluid production at a point, 
determining flow patterns in a reservoir, determining 
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optimal location of steam or water injection wells, and 
determining the optimal depth at which to produce 
from a reservoir within a given well. 

An illustration of a reservoir model incorporating the 
above features is shown in FIG. 13. The reservoir 
model, shown generally by (110) is comprised of vari- 
ous layers (11a-11d) of varying permeabilities made up 
of sand or glass beads. Although not shown for simplic- 
ity, the layers would be confined by a vessel enclosing 
at least the four sides and bottom of the reservoir model. 
The reservoir model is fitted with a fill line (66) so that 
liquids may be added to the model and a drain line (112) 
for extracting liquid or draining the model. The layers 
(111a-111d) are contoured to simulate the conditions to 
be modeled. Within the model are several glass tubes 
(1134-113) closed at the bottom and open at the top 
above the uppermost layer (111a). These tubes can 
either be production points or source points, as dis- 
cussed above. In the illustration, tube (113/) acts as a 
source point and has a bright isotropic light (69) sus- 
pended from electrical supply lines (70) acting as a 
source of simulated fluid. Tubes (113d) and (113e) act as 
monitoring points and have light detectors (114) and 
(115) inserted in them. Tubes (113a), (113g) and (113A) 
can be painted flat black to absorb light and simulate 
production points. The shadowed regions behind these 
black tubes simulate real-world reservoir conditions 
where the fluid flow is similarly distributed at produc- 
tion points. The detectors send a signal proportional to 
the light received at the detector to a meter or recorder 
(not shown) via leads (116) and (117), respectively. 
These detectors can be located at different depths 
within the tube to determine the variation in simulated 
source fluid as a function of depth within the model. 
For example, by use of a servo (118) connected to de- 
tector (115) by the detector leads (116), the position of 
the detector (115) can be varied within the tube (113e). 
The position of the detector (115) and the correspond- 
ing light reading could be recorded by a recording 
means, as for example the strip-chart recorder (119), or 
could be sent to a computer (23) for further processing. 
While the isotropic light (69) provides one means of 
simulating the source, an alternate means would be a 
laser. In this variation a laser (3) generates a coherent 
beam of radiation (2) which is projected down the glass 
tube (1135). A mirror (120) located within the tube is 
angled so as to reflect the beam in a vectored direction 
as shown by vector (121). The benefit of using the laser 
as the simulated source is that much greater transmis- 
sion distances can be achieved. It should be observed 
that actual fluid injection points such as described in 
FIG. 11 can be incorporated into the type of model 
shown by FIG. 13 to collect data from moving fluid 
interface modeling. 

While the above embodiments of the invention dis- 
cuss use of an essentially continuous light source, it will 
be observed by one skilled in the art that the measure- 
ments. discussed above can be performed in the time 
domain by using a pulsed light or a pulsed laser and an 
optical oscilloscope. 

From the above description it is evident that the 
present invention provides a method and apparatus for 
measuring properties of porous media by use of light or 
laser transmission through the media. Although only 
specific embodiments of the present invention have 
been described in detail, the invention is not limited 
thereto but is meant to include all embodiments coming 
within the scope of the appended claims. 
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What is claimed is: 

1. A permeability mapping apparatus comprising: 

(a) a matrix medium of closely packed particles; 

(b) a pore filling medium symbiotically interposed 
within said matrix medium; 

(c) a means for generating a coherent collimated 
beam of electromagnetic radiation, said beam of 
electromagnetic radiation vectored to transpierce a 
series of locations defining said matrix medium; 

(d) a means for successively generating an attenuation 
signal for each of the vectored beams of electro- 
magnetic radiation at each location transpierced by 
said beam of electromagnetic radiation, and assign- 
ing a weighted value based on a correlation func- 
tion of known permeability to each of the attenua- 
tion signals so as to compensate each beam attenua- 
tion signal for attenuations at other than one.of said 
matrix locations; 

(e) a means for reconstructing all said weighted val- 
ues over said series of locations defining said matrix 
so as to reconstruct said matrix in terms of each 
attenuation signal at each of said locations; 

(f) a means for displaying each reconstructed value to 
image said matrix medium. 

2. The apparatus of claim 1 wherein the matrix me- 
dium is embodied in a compartment, said compartment 
having a means for allowing both ingress and egress of 
electromagnetic radiation. 

3. The apparatus of claim 1 wherein the pore-filling 
medium has a refractive index greater than the matrix 
medium. 

4. The apparatus of claim 1 wherein the pore-filling 
medium has a refractive index less than the matrix me- 
dium. 

5. The apparatus of claim 1 wherein the matrix me- 
dium and pore-filling medium are optically transparent. 

6. The apparatus of claim 1 wherein the matrix me- 
dium and pore-filling medium are optically translucent. 

7. The apparatus of claim 1 wherein the means for 
generating a coherent collimated beam of electromag- 
netic radiation is a laser source. 

8. The apparatus of claim 1 wherein the means for 
generating a coherent collimated beam of electromag- 
netic radiation is a parabolic reflector and lens combina- 
tion used in conjunction with an isotropic light source. 

9. The apparatus of claim 7 wherein the wavelength 
of the electromagnetic radiation is at least one order of 
magnitude smaller than said matrix medium’s mean 
grain diameter, preferably said wavelength being two to 
three orders of magnitude smaller than said mean grain 
diameter. 

10. The apparatus of claim 9 wherein the matrix me- 
dium is Boise sandstone. 

11. The apparatus of claim 9 wherein the matrix me- 
dium is between 1 mm and 20 mm thick, preferably 
between 2 and 10 mm thick. 

12. The apparatus of claim 11 wherein the means for 
successively generating an attenuation signal is a light 
detector system, said detector system in general align- 
ment with the laser source along a line perpendicular to 
the laser’s collimated beam, said detector system posi- 
tioned on a side opposite of that which is transpierced 
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28 
by the laser’s collimated beam, and generating attenua- 
tion signals proportional to the total attenuation of 
beams passing through said matrix at a plurality of loca- 
tions defining said matrix. — 

13. The apparatus of claim 12 wherein the means for 
reconstructing all weighted values in response to said 
beam attenuation signals is a digital computer compris- 
ing: 

(a) a means to convert said attenuation signals into 

digital signals; 

(b) a means to store said digital signals; 

(c) a means responsive to said stored digital signals 
for converting digital signals into corresponding 
analog signals proportional to said derived 
weighted values. 

14. The apparatus of claim 13 wherein the means for 
displaying each reconstructed value is a cathode ray 
tube means, responsive to said analog signals, for gener- 
ating a pictorial representation of all matrix locations 
transpierced by said beam with a proportional intensity. 

15. The apparatus of claim 1 wherein permeability 
changes over time are determined, further comprising: 

a sealed container having a fluid inlet and fluid outlet, 
said container housing a matrix section; and 

a first fluid, and a second fluid to displace the first 
fluid within the container, said second fluid react- 
ing with clays present in said matrix section to 
cause a change in permeability. 

16. The apparatus of claim 1 wherein grain size and 
grain size variation of the matrix medium are measured, 
and where said means for generating an attenuation 
signal is a light detection system comprising a plurality 
of small light detectors in a fixed array. 

17. The apparatus of claim 16 wherein the light detec- 
tors have an area preferably 2 mm X2 mm. 

18. A method modeling fluid flow within a matrix 
medium comprising the steps of: 

constructing a sealed representation of all zones of 
differing permeability within a matrix medium; 

interposing a light enhancing pore-filling fluid to 
symbiotically interact with said matrix medium; 

flooding said matrix medium with an isotropic light 
source; 

measuring the attenuation of said isotropic light at a 

“ plurality of locations within said matrix, said atten- 
uation being correlatable to steady state fluid flow. 

19. The method of claim 18 wherein dynamic flow 
patterns are modeled further comprising the steps of: 

injecting a second fluid into said matrix to drive the 
initial pore-filling fluid; 

measuring both the change in attenuation of said 
isotropic light at the interface of said second fluid 
and said measuring means, and the time of arrival 
of said second fluid interface at said measuring 
means; 

generating a prediction of fluid flow based on mea- 
sured attenuation as a function of the light source 
and attenuation measuring means properties reflec- 
tion coefficients at the fluid interface, and the pho- 


ton means free path in each zone. 
* * * & * 
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(54) Measuring and treatment tool incorporating broadband stripline aerials and useful in medical technology 


(57) The tool, for determining the nature of fluid in mammal tissue or cancer therapy by hyperthermia transmits, and 
optionally receives electromagnetic radiation over a broad band of 2KHz — 1GHz and can determine resistivity and dielectric 
constant. it incorporates at least one transmitting and receiving antenna. Each antenna comprises a coaxial cable 151 
connected to stripline adapter 153, which is connected to a stripline (155) having a metallic central strip (159). A strip face 
161 is bent at approximately right angles, and has a length that is compatible with the desired frequency coverage. A 
ground plane 165 extends from the stripline adapter to the right angle bend, and a dielectric 167 fills the space between the 
centre strip and the ground plane. An enclosure comprising four metallic walls 181 surrounds the stripline, and is in 
electrical contact with the ground plane and the stripline adapter and a lossless, non-conducting material fills the enclosure. 
The antennae are positioned so that the strip face lies flush with the tool face, to permit electromagnetic energy to be 
transmitted into and out of the material to be analyzed. Amplitude and phase of the received signal are monitored. 
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At teast one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 


V 080 1S¢ ¢ dD 


1/7 


Ht i0| 


FIG_1 


163 
155 
159 161 
~\ KOSS Bi 167 
er 
SSS5SS5N5555 SZ 


(Wate 


LLL 


Le 


FIG_2 


FIG_2A 


FIG_S 


ot 7 


BROADBAND ANTENNA BROADBAND ANTENNA 
IN AIR IN'E' 


I NYNLSY 
t a 
ADaAOW 
P) SSOT NYNL3SY 
— 
: 
ow 
XD 
> 
m 


A -20 -2 
-25 -25 
a -30 -30 
=o5 738 
sss3ss 833338 
06 8888 8 08 8 o 0 98 
FREQUENCY (MHz) FREQUENCY (MHz) 
4 BROADBAND ANTENNA 
2 IN 'E' BRINE EXPERIMENTAL 
z (NO METAL PLATE) SET-UP 
= 
a . TRANSMITTING RECEIVING 
Y 219 ANTENNA ANTENNA 
O 
2-20 
5 35 
~ -50 TOOL 
~~60 PAD 
-70 
~-80 Ni ff, 


i) mn © O 
S$ 86888 
00606066 


OOZI 


(REFLECTOR) 


FREQUENCY (MHz) Vi Lisi 
FIG_4C FIG_4D 


G™9ls 
(ZHW) AONANDSYS 


00g Osi O91 Ovi oz oOo1 os O09 Ov O02 C302 


06- 


O8- 


6/4 


OL- 
09- 
SSO1 NOISSINSNVYL 
OSs- wep 


Ov- 


ZHH S=!i SSO1 NYUNLSY 7 O¢- 


WEP £ YOJ AONSNOIYS YAMO) 
| Oc- 
SSO1 NYNL3ay 
Oi- 


O- 


(SLTOAITTIN) S9VLIOA 


(SLIOAITIUIW) SOVINOA 


TRANSMITTING THROUGH 
'e' BRINE -2" 
- FROM REFLECTOR 


30 
25 
20 
15 
1) 
5 
e) 


- eits ES GER aie a 


0.0 10 20 30 40 507 
TIME (MICROSECONDS ) 


FIG_6A 


TRANSMITTING THROUGH 
'e' BRINE - 4" 
FROM REFLECTOR 


SLIOAITIIN) SOVLIOA 


SLIOAITTIAN) SOVLIOA 


( 


0.0 10 20 30 40 50 
TIME (MICROSECONDS ) 


FIG_6C 


Paar | 


TRANSMITTING THROUGH 
'E' BRINE -3" 
FROM REFLECTOR 


5 
00 10 20 30 40 50 
TIME (MICROSECONDS ) 


FIG_6B 


TRANSMITTING THROUGH 


‘e' BRINE 

35 
30 
as NO REFLECTOR 
20 

15 

10 

5 

oO 
-5 


00 10 20 30 4050 
TIME ( MICROSECONDS ) 


FIG_6D 


01 
02 
03 
04 
0S 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 


2251030 
ais 


METHOD AND APPARATUS FOR BROADBAND ELECTROMAGNETIC 
ENERGY COUPLING 


The present invention relates generally to the 
electromagnetic coupling and analysis. More specifically, 
this invention provides an antenna which can combine the 
functions of various resistivity and dielectric constant 
devices into a single tool, capable of operating over a wide 
range of frequencies. It is particularly useful in the 
field of medical technology. 


In the field of medical technology, it is well known that 
electromagnetic energy is useful in various types of 
diagnoses and treatments. For example, recent statistics 
show that pulmonary and cardiopulmonary diseases are 
responsible for more than three million hospital admissions 
and 30,000 deaths every year in the United States. 
Pulmonary abnormalities are virtually always associated with 
changes in lung water content or distribution. Most workers 
agree that there is no single, nondestructive method 
available to detect accurately early changes in lung 

water content. 


For a clinically useful technique, it is desirable to detect 
early changes in both the extravascular lung water, which 
strongly reflects most pulmonary abnormalities, and the 
intravascular compartment. Recently, the use of the 
electromagnetic methods to detect changes in lung water 


content have shown promising initial results, particularly 


for detecting small variations in water content. 
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Particularly at microwave frequencies, changes in the 
dielectric properties of tissue are closely related to the 
amount of water present. Electromagnetic techniques, 
therefore, basically utilize changes in the permittivity and 
conductivity of the lung region to detect changes in lung 
water content. This method has the advantage of using 
highly penetrating electromagnetic signals rather than 
ultrasonic signals which are both highly attenuated and 
dispersed in the lung. Furthermore, electromagnetic 
techniques have the potential for continuous monitoring of 
patients in intensive care units, such as those with heart 
failure or extensive burns. 


U.S. Patent No. 4,240,445 issued to Iskander et al. and is 
incorporated herein by reference for all purposes. Iskander 
teaches a method of coupling electromagnetic energy into a 
material such as tissue, to measure water content. 

Measuring lung water content is an especially useful 
application. However, Iskander’s device is so large that 
only a few antennas can be place on the chest, and the 
antenna cannot be described as a point source. Also, the 
electric field vanishes at some distance from the antenna, 
as the electric fields in the two parallel slots are 
oppositely directed. Furthermore, a resistor is included in 
the antenna, which dissipates much of the electromagnetic 
energy in the antenna itself and introduces a limitation in 
the power handling capability of the antenna. Additional 
prior work includes: M. F. Iskander and C. H. Durney 
(1980): “Electromagnetic Techniques for Medical Diagnosis: 
A Review", Proceedings of IEEE, vol. 68, no. 1. and 

M. F. Iskander et al (1982): “Two-dimensional Technique to 
Calculate the EM Power Deposition Pattern in the Human 
Body", Journal of Microwave Power, vol. 17, no. 3. There is 
thus a need for a device that is compact enough to permit 
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placing of many antennas forming an array on a chest to 
obtain a well-defined image of the chest cavity, a device 
that has an antenna. that can be mathematically described as 
a point source, and one which does not suffer from 
cancellation of the electric field at a certain distance. 


A dielectric transmitting and measuring device can also be 
used to heat an interior portion of a mammalian body to 
destroy or reduce the size of tumors. Tumor reduction 
therapy, or cancer therapy by hyperthermia, combined with 
radiation or drugs is known in the art to either stop or 
slow down the growth of cancer cells, or cause the death of 
the cancer cells. (See, for example, Streffer, C., "Cancer 
Therapy by Hyperthermia and Radiation", Urban and 
Schwarzenberg, Munich, F.R.G., 1980 and Dethylefsen, L.A. 
(Editor), "The Third International Symposium: Cancer 
Therapy by Hyperthermia, Drugs and Radiation, Colorado State 
University, Ft. Collins, U.S.A., 1980.) 


One such device is disclosed by J. Scheiblich et al. 
"Radiofrequency-Induced Hyperthermia in the Prostate", 
Journal of Microwave Power, vol. 17, no. 3, 1982, Ottawa, 
Canada. Scheiblich et al’s device works only at a 
single frequency. 


A propagating electromagnetic wave has two fundamental 
characteristics, amplitude and phase. By comparing the 
amplitude and phase of an electromagnetic wave as it passes 
receivers, propagation characteristics of the probed medium 
may be studied. Measurement of these two characteristics 
may be used to determine the dielectric constant and 

the conductivity of the media through which the wave 


is propagated. 
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However, no one tool in the prior art is capable of probing 
or coupling energy into a material over a broad band of 
frequencies. It is therefore advantageous to extend the 


frequency range. 


The largest hurdle to developing such a broadband dielectric 
tool has been the lack of a suitable broadband antenna that 
can couple electromagnetic energy to and from a material, 
and that is compact enough to fit within the confines of 

a tool. 


The prior work is limited in the attempts at electromagnetic 
coupling, analysis, and treatment, in that no suitable 
single antenna element has been designed which can couple 
electromagnetic energy into a material, such as mammal 
tissue, over a broad range of frequencies, that is also 
sufficiently compact and is capable of handling high power 
levels. There is therefore a need for a device and a method 
for use in such broadband applications. 


The present invention is surprisingly successful in 
providing a method and apparatus for combining the functions 
of various conductivity and dielectric constant devices and 
electromagnetic energy coupling devices into a single 
device, capable of operating over a wide range of 
frequencies. It is especially useful in medical 

technology applications. 


A measuring or electromagnetic coupling tool, having a tool 
face, also has a novel transmitting antenna and a novel 
receiving antenna. Electromagnetic energy is transmitted to 
a transmitting antenna. A stripline adapter permits the 
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energy to flow to a stripline having a metallic central 
strip. A strip face of the central strip is bent at 
approximately right angles, and has a height that is 
compatible with desired frequency coverage. 


A ground plane extends from the stripline adapter to the 
right angle bend, so that a distal end of the central strip 
extends away from it, and a void is created between the 
center strip and the ground plane. 


A dielectric is positioned to nearly fill the void. The 
dielectric is comprised of a material having a very high 
dielectric constant and a very low energy loss. The 
transmitting antenna is positioned so that the ground plane 
is fixedly connected to the measuring tool, and the strip 
face lies flush with the tool face, so that electromagnetic 
energy can be transmitted into the material to be analyzed. 


An enclosure surrounding the stripline is comprised of four 
metallic walls which are positioned in electrical contact 
with the ground plane and the stripline adapter, so that the 
strip face is nearly centered in the opening created by the 
walls and the ground plane. 


A loss-less, non-conducting material fills in any remaining 
open space in the enclosure, so that the non-conducting 
material forms an additional wall that is really flat with 
the strip face. 


A receiving antenna is comprised in essentially the same 

manner as the transmitting antenna, and is positioned in the 
tool so that it can receive the electromagnetic energy which 
has traveled through the material being probed. A means for 
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monitoring the received energy detects the phase and 
the amplitude. 


In another embodiment of this invention, broadband 
measurements are taken to determine the quantity of a fluid 
in a material, such as water in a lung. 


It is one object of this invention that electromagnetic 
energy is transmitted and received over a wide frequency 
range, specifically from a few KHz to a few GHz. A commonly 
used frequency range is from 2 KHz to 4 GHz. 


The tool may further comprise a pad, which substantially 
conforms to the surface of the mammal tissue, and holds the 
antennas. At least one transmitting antenna is necessary. 
No receiving antenna is necessary, although a plurality of 
each is often desirable. 


The above and other embodiments, objects, advantages, and 
features of the invention will become more readily 
apparent from the following detailed description of the 
invention, which is provided in connection with the 
accompanying drawings. 


Figure 1 is a schematic, sectional view of the inventive 
device positioned adjacent to mammal tissue. 


Figure 2 shows a top, front, and side view of the novel 
transmitting antenna. 


Figure 2A is the same view as Figure 2, further illustrating 
the enclosing metallic walls. 
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Figure 3 shows an antenna mounted on a tool face. 


Figure 4 shows three graphs of transmission and return loss 


as a function of frequency. 


Figure 5 ig a graph of transmission and return loss as a 
function of frequency, for low frequencies. 


Figure 6 shows four graphs of time-domain transmission 
measurements at various distances from a metal reflector 


plate in a brine. 


In accordance with the present invention, a new improved 
method and apparatus for coupling electromagnetic energy 
into a material for determining the nature of various 
materials and the fluids contained therein and to induce 
hyperthermia, using a broadband measuring apparatus, has 
been developed. 


Referring to the drawings, a first embodiment of the 
inventive broadband tool 101 is shown in Figure 1, 
positioned around a portion of a mammal body such as a chest 
cavity 103. A means such as a belt mount 109 positions tool 
face 111 near the mammal skin 104, such that transmitting 
antennas such as Tl and T2 and receiving antennas such as Rl 
and R2 are positioned touching the skin surface of 104. The 
tool face 111 is defined as the surface of the belt 

mount 109 containing the aperture plane of the antennas, and 
is preferably a continuous metallic surface. The belt 
mount 109 may be made of any suitable flexible material that 
can be strapped around the portion of interest of the mammal 
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body. A conducting compound such as a conducting grease may 
be applied at the interface 113 between the tool face 111 
and the skin surface 104 to improve coupling between the 
antennas and the chest cavity 103. 


The region of the mammal body to be investigated may not be 
electrically homogeneous. In the chest cavity 103 for 
example, there are organs such as the heart 105, the lung 
region 106, the vertebra 107, and there may also be a 
tumor 108. It is often desirable to analyze or treat 
selected portions of such a cavity 103. 


An analysis of the chest cavity 103, for example, can be 
done by a dielectric imaging of the cavity. This is done by 
transmitting electromagnetic energy at a suitable frequency 
across the chest cavity 103 from a transmitting antenna such 
as Tl, and receiving this energy at a receiving antenna such 
as Rn. In this way the phase and the amplitude of the 
propagated electromagnetic wave for the path T1Rn (shown in 
dashed line) is determined. Since there can be a 
multiplicity of transmitting antennas Tn and a multiplicity 
of receiving antennas Rn, a multiplicity of such paths 
crisscrossing the entire chest cavity can be studied. From 
this information, using well known techniques, a dielectric 
image of the chest cavity can be generated. Such an image 
displays the various organs in the cavity, and when suitably 
made, can reveal the presence of tumor 108. The dielectric 
properties, and thus a dielectric image, can be determined 
as a function of position within the material being probed. 
Since dielectric image is very sensitive to the presence of 
water, it can also give an assessment of the lung water 
content; Cf. "Microwave Methods of Measuring Changes in Lung 
Water", by M. F. Iskander and C. H. Durney, Journal of 
Microwave Power, vol. 18(3), 1983, p. 265. 
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Note that although the antennas have been labeled as either 
transmitting or receiving antennas, any given antenna can 
serve either function. 


The broadband capability of the antennas is an advantage in 
the above applications for the following reasons: 
structures (e.g., heart, tumor) of different sizes require 
different frequencies for their best definition in the 
image; highly lossy regions such as fluids may require 
employment of relatively low frequencies so that the 
electromagnetic losses are acceptable; in time-domain 
application, simultaneous information at a multiplicity 

of frequencies can be developed. 


In the treatment mode, it is desirable to reduce or 
eliminate the tumor 108 by hyperthermia, i.e., by 
selectively heating only the tumor region 108 to a high 
temperature. Thus, by selecting a suitable group of 
antennas to transmit, one can selectively deposit 
electromagnetic energy in the region of the tumor 108; Cf. 
"Two-dimensional Technique to Calculate the EM Power 
Deposition Pattern in the Human Body", by M. PF. Iskander, 
P. F. Turner, J. B. DuBow and J. Kao, Journal of Microwave 
Power, vol. 17(3), 1982, p. 175. 


The broadband capability of the antennas is an advantage in 
the above application because for a given situation, one can 
select the frequency that simultaneously produces the 
optimum deposition of power and localization of the heating 
using known techniques. 


An example of the inventive transmitting antenna 150 is 
shown in Figure 2. A coaxial connecting means, such as 
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coaxial connector 151 is electrically connected to a 
stripline adapter 153 which is capable of transmitting 
electromagnetic energy from the coaxial connector 151 to a 
stripline section with metallic central strip 155. An 
especially useful stripline adapter is a model No. 
3070-1404-10 designed by Omni-Spectra, or other types of 
microwave stripline adapters. Other types of transmission 
means may be utilized to transmit electromagnetic energy to 
the antenna. For example, a strip transmission line may be 
electrically connected to the stripline section having the 
metallic central strip 155. As a commercial 
coaxial-to-stripline transition means has been utilized, the 
dimensions included herein reflect this means. One 
knowledgeable in the art would realize that the 

dimensions may be altered to change frequency coverage 

and to fine-tune performance. 


Metallic center strip 155 has a front end 157, a flat strip 
body 159, a flat strip face 161, and a distal end 163. The 
front end 157 is electrically connected to the center 
conductor 169 of the stripline adapter 153. Solder is a 
particularly useful connecting means. Flat strip body 159 
may also be tapered to come to a point at front end 157 to 
provide a smooth electrical transition between the center 
conductor 169 and the center strip 155. The strip face 161 
is bent at approximately right angles to strip body 159, and 
has a height that is measured from the right angle bend to 
distal end 163. The height is compatible with the desired 
frequency coverage. The longer the height, the more lower 
frequency coverage is allowed. A \" height permits a 
frequency range of approximately 2 KHz > 1 GHz. A Sma 
height extends the upper frequency limit to approximately 

2 GHz. An upward frequency limit of 4 GHz is attainable as 


01 
02 
03 
04 
0s 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 


=Tt= 


well. The metallic center strip 155 can be made of any 
metal. Copper, brass, or aluminum are especially useful. 


A ground plane 165 extends from stripline adapter 153 to the 
right angle bend in the center strip 155, so that the distal 
end 163 extends away from the ground plane 165 and so that a 
void exists between the center strip 155 and the ground 
plane 165. Ground plane 165 is comprised of a metal. 
Commercial grade stainless steel is particularly useful. It 
is desirable to keep the ground plane and center strip as 
short at possible, to permit the apparatus to remain as 
compact as possible and to allow the use of as many antennas 


as possible. 


The void between the ground plane 165 and the center 
strip 155 is largely filled with a dielectric 167. The 
Gielectric 167 should have a very high dielectric constant 
and a very low loss. By loss, we mean the dissipation of 
energy. The dielectric 167 can be a ceramic dielectric, and 
comprised of material such as Barium Titanate or Lead 
Zirconate Titanate. A crystalline dielectric may also be 
used, although more expensive. The thickness of the 
dielectric 167 is determined by the stripline adapter 153 
used. The dielectric 167 acts to make the capacitance of 
the center strip 155 very large. 


The construction of the antenna is completed by enclosing 
the center strip 155 by metallic walls 181, 182, 183, and 
184, which contact the ground plane 165 and the adapter 153 
electrically, as shown in Figure 2A. The walls add rigidity 
and prevent leakage of the electromagnetic radiation. The 
strip face 161 is approximately centered in the rectangular 
opening created by the edges of the walls and the edge of 
the ground plane 165. Thus, the distance between an edge of 
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the strip face 161 and the adjacent edge of a wall is 
substantially the thickness of the dielectric 167. The 
entire void space in the antenna enclosed by the walls, 
including the set back 168 at the dielectric edge, is filled 
with a loss-less, non-conducting material such as a mixture 
of epoxy and alumina which sets hard, seals the antenna, and 
makes it more rugged. 


The ground plane 165 and the walls 181, 182, 183, and 184 
are fixedly connected to an electromagnetic coupling or 
analyzing tool as seen in Figure 3. The strip face 161 is 
positioned to lie flush with the tool face 171 (which is the 
same as the tool face 111 of Figure 1), so that the 
transmitting antenna 150 can transmit electromagnetic energy 
into a material such as mammal tissue. A conductive 
substance, known in the art, is usually placed on the 
outside of the mammal tissue, to permit a sufficient flow of 
electromagnetic energy into the tissue. void space 173 is 
filled with a loss-less, non-conducting material such as an 
epoxy-aluminum compound. The ground plane 165 and the walls 
181, 182, and 183 connect to the tool face. 


A receiving electromagnetic antenna is comprised in 
essentially the same manner as the transmitting antenna, and 
is positioned in the tool in the same manner as the 
transmitting antenna, so that the receiving antenna can 
receive the electromagnetic energy which has traveled 
through the material that is analyzed. 


The present invention is especially useful in the field of 
microwave diagnostics of fluid content and fluid quantity. 
For example, the apparatus can couple electromagnetic energy 
into mammal tissue. The electromagnetic energy can be 
monitored to provide an indication of the amount and 
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distribution of a fluid, such as water, inside the mammal 
tissue. One particularly useful application is to measure 
the water content in a lung. The present apparatus is very 
compact, and therefore requires a much smaller skin contact 
area. Also, many antennas can be placed on a chest cavity, 
to obtain a well defined image of the chest cavity. The 
inventive antennas can be mathematically described as a 
point source, thus making analysis of the data easier. A 
conductive substance should be placed on the outside of the 
chest cavity, to permit a sufficient flow of electromagnetic 
energy into the chest cavity. 


The prior art (Iskander et al.) has the drawback that the 
electric field vanishes at some distance from the tool face, 
since the fields in the two parallel slots are oppositely 
directed. No such cancellation occurs with the present 
invention. Furthermore, the incorporation of a resistor in 
Iskander et al’s antenna introduces a power limitation. 


In another embodiment, the present invention can be used in 
the field of microwave hyperthermia. The apparatus can 
couple electromagnetic energy into the interior portion of a 
mammal, so that the electromagnetic energy is focused to 
heat and thereby reduce the size of or destroy a tumor. 
Tumor reduction therapy or cancer therapy, by hyperthermia, 
combined with radiation or drugs is known in the art to 
either stop or slow down the growth of cancer cells, or 
cause the death of the cancer cells. 


The present invention has the advantage over the prior art 
that many frequencies can be selected. Because there is no 
limitation to the power handling capability in the inventive 
antenna, the present invention is particularly suited for 
depositing microwave power into a localized area inside a 
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mammal, such as a human. Either a single antenna or an 
array of antennas could be used. 


In yet another embodiment, the apparatus can be implanted 
inside the body of a mammal, and used as a radio frequency 
antenna. Either a single antenna or an array of antennas 
could be used. As the inventive antenna can be made very 
small (as small as approximately 10 mm long and 
approximately 5 mm high), it is particularly suitable to 
this application. As the antenna gets smaller, the 
frequency coverage shifts to higher frequencies. The 
apparatus can be constructed with a commercial smicro-coaxial 
connector. However, smaller devices can be constructed 
through the use of a customized coaxial connector. 


The apparatus can operate in the frequency domain, using a 
Single frequency, multiple frequencies (such as 
simultaneous, selectable, or time-multiplexed for example), 
or swept frequency techniques. Or, the apparatus can 
operate in the time domain, using pulses of a wide variety 
of shapes, widths, rise and fall times, etc. When the 
pulses are transformed to the frequency domain, either 
electronically using a spectrum analyzer, or numerically 
using mathematical transforms, the same information is 
obtained as would be given by a frequency domain tool. 


A prototype tool was constructed, with the inventive 
antennas. The tool consists of one transmitting and one 
receiving antenna, the distance between them being variable. 


An acceptable dielectric antenna must meet the following 
criteria: 
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(i) It must be able to couple sufficient energy into 
and from the material at its operating frequency 
to allow probing of the material; 


(ii) This probing energy must penetrate into the 
material, rather than clinging to the surface of 
the tool (i.e., it must travel as a freely 
propagating wave rather than a surface wave guided 
along the tool face). 


In the present instance, the above two conditions must hold 


over the entire range of the frequency of operation. 


The first of the above criteria is tested by measuring the 
return loss for the transmitting antenna, and the 
transmission loss from the transmitting to the receiving 
antenna - both as a function of frequency. These 
measurements are shown in Figure 4 where the tool is placed 
in air and against brine of conductivity 0.5 mho/m (to 
represent a biological medium). The return loss curve in 
brine shows that sufficient energy is entering the brine 
over the frequency range of the measuring device 
(Hewlett-Packard HP8505A Network Analyzer; 500 KHz - 

1300 MHz) to permit probing. The transmission loss shows 
that sufficient energy is being received at the receiving 
antenna to permit measurements. 


Measurements were made by using another measuring device 
(HP3577A Network Analyzer; 5Hz - 200 MHz) to test the low 
frequency limitation of the antenna. The results are shown 
in Figure 5, showing that the low frequency limitation is 
about 5 KHz. The improved return loss performance in the 
200 MHz region (at Figure 4) results from a drying (curing) 
of the epoxy alumina filling between measurements. 
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Figure 6 shows time-domain transmission measurements at 
various distances (d) to a metal reflector plate in the 
brine. The change in amplitude of the received pulse as a 
function of the distance of the metallic reflector shows 
that the energy has penetrated into the brine out to the 
location of the plate. 


While a preferred embodiment of the invention has been 
described and illustrated, it should be apparent that many 
modifications can be made thereto without departing from the 
Spirit or scope of the invention. Accordingly, the 
invention is not limited by the foregoing description, but 
is only limited by the scope of the claims appended hereto. 
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CLAIMS: 


Apparatus for coupling electromagnetic energy into 
materials comprising a measuring tool having a tool 
face, said measuring tool further comprising an 
electromagnetic transmitting antenna, said transmitting 
antenna further comprising: 


(a) 


(b) 


(c) 


(d) 


a coaxial cable connecting means and means to 
transmit electromagnetic energy therethrough; 


a stripline adapter capable of transmitting 
electromagnetic energy from said coaxial cable 
connecting means to a stripline having a metallic 
central strip, said center strip having a front 
end, a flat strip body, a flat strip face, and a 
distal end, said front end electrically connected 
to a center conductor of said stripline adapter, 
said strip face bent at approximately right angles 
to said strip body and having a height measured 
from said right angle bend to said distal end that 
is compatible with a desired frequency coverage; 


a ground plane which extends from said stripline 
adapter to said right angle bend, so that said 
distal end extends away from said ground plane and 
so that a void exists between said center strip 
and said ground plane; 


a dielectric largely filling said void, said 
dielectric comprised of a material having a very 
high dielectric constant and a very low energy 
loss, so that said transmitting antenna is 
positioned so that said ground plane is fixedly 


(e) 


(f£) 


(g) 


(h) 
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connected to said measuring tool and said strip 
face is positioned to lie flush with said tool 
face so that said transmitting antenna can 
transmit electromagnetic energy into said 
material; 


an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in 
electrical contact with said ground plane and said 
Stripline adapter, so that said strip face is 
nearly centered in the opening created by said 
walls and said ground plane; 


a loss-less, non-conducting material which fills 
in any remaining open space in said enclosure so 
that said non-conducting material forms an 
additional wall that is nearly flat with said 
strip face; 


said receiving electromagnetic antenna comprised 
in essentially the same manner as said 
transmitting antenna, said receiving antenna 
positioned in said measuring tool in the same 
manner as said transmitting antenna, so that said 
receiving antenna can receive said electromagnetic 
energy which has traveled through said material; 
and 


means for monitoring the amplitude and the phase 
of said received electromagnetic energy. 


Apparatus as recited in Claim 1 further comprising a 
means for positioning said tool face near said 
material. 


10. 


279.6 


Apparatus as recited in Claim 1 wherein said 
electromagnetic energy is focused to heat and thereby 
reduce the size of a tumor in a mammal. 


Apparatus as recited in Claim 1 wherein said 
electromagnetic energy is focused to heat and thereby 


destroy a tumor in a mammal. 


Apparatus as recited in Claim 1 wherein said antennas 
are positioned on a belt-mounted device. 


Apparatus as recited in Claim 1, further comprising a 
plurality of receiving antennas. 


Apparatus as recited in Claim 6 further comprising a 
plurality of transmitting antennas. 


Apparatus as recited in Claim 7 wherein said materials 
are mammal tissue and water. 


Apparatus as recited in Claim 1, wherein broadband 
measurements are taken to determine said dielectric 
properties as a function of position within 

said material. 


An apparatus for coupling electromagnetic energy to 
determine the quantity of a fluid in a material, said 
apparatus having a tool face and further comprising a 
first electromagnetic transmitting antenna, said first 
transmitting antenna further comprising: 


(a) a coaxial cable connecting means and means to 
transmit electromagnetic energy therethrough; 
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(d) 


(e) 
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a stripline adapter capable of transmitting 
electromagnetic energy from said coaxial cable 
connecting means to a stripline having a metallic 
central strip, said center strip having a front 
end, a flat strip body, a flat strip face, and a 
distal end, said front end electrically connected 
to a center conductor of said Stripline adapter, 
said strip face bent at approximately right angles 
to said strip body and having a height measured 
from said right angle bend to said distal end that 
is compatible with a desired frequency coverage; 


@ ground plane which extends from said stripline 
adapter to said right angle bend, so that said 
distal end extends away from said ground plane and 
so that a void exists between said center strip 
and said ground plane; 


a dielectric filling most of said void, said 
dielectric composed of a material having a very 
high dielectric constant and a very low energy 
loss, so that said first transmitting antenna is 
positioned so that said ground plane is fixedly 
connected to said logging tool and said strip face 
is positioned to lie flush with said tool face so 
that said first transmitting antenna can transmit 
electromagnetic energy into said material; 


an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in 
electrical contact with said ground plane and said 
stripline adapter, so that said strip face is 


01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 


11. 


12. 


13. 


a9 Ve 


nearly centered in the opening created by said 
walls and said ground plane; 


(f£) a loss-less, non-conducting material which fills 
in any remaining open space in said enclosure so 
that said non-conducting material forms an 
additional wall that is nearly flat with said 
strip face; 


(g) said receiving electromagnetic antenna comprised 
in essentially the same manner as said 
transmitting antenna, said receiving antenna 
positioned in said apparatus in the same manner as 
said transmitting antenna, so that said receiving 
antenna can receive said electromagnetic energy 
which has traveled through said material; and 


(h) means for monitoring the amplitude and the phase 
of said electromagnetic energy, so that the 
quantity of said fluid can be determined. 


Apparatus as recited in Claim 1 or 10 wherein said 
transmitting antenna transmits and said receiving 
antenna receives electromagnetic energy over a 
frequency range of 2 KHz to 4 GHz. 


Apparatus is recited in Claim 1 or 10 wherein said 
transmitting antenna can alternately function as a 
receiving antenna and said receiving antenna can 
alternately function as a transmitting antenna. 


Apparatus as recited in Claim 12 further comprising a 
belt-mount, said belt-mount substantially conforming to 
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the outside of a mammal tissue and holding said 


transmitting and receiving antennas. 


Apparatus as 
plurality of 


Apparatus as 
plurality of 


Apparatus as 
are water. 


Apparatus as 
ground plane 


Apparatus as 


recited in Claim 13 further comprising a 
receiving antennas. 


recited in Claim 14 further comprising a 
transmitting antennas. 


recited in Claim 10 wherein said fluids 
recited in Claim 1 or 10 wherein said 
is no greater than 10 mm in length. 


recited in Claim 1 or 10 wherein said 


strip face has a height that is no greater than 5 mm. 


Apparatus as 


recited in Claim 1 or 10 wherein said 


electromagnetic energy is monitored to provide an 
indication of the amount and isecibubion of a fluid 


inside mammal tissue. 


Apparatus as 


recited in Claim 1 or 10 wherein no 


receiving antenna is incorporated. 


Apparatus as 
apparatus is 


recited in Claim 1 or 10 wherein said 
implanted inside a mammal, as a radio 


frequency antenna. 


Apparatus as 


recited in Claim 21 wherein said apparatus 


does not incorporate a receiving antenna. 
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Apparatus as recited in Claim 1 or 10 wherein a strip 
transmission line is electrically connected to said 
stripline, so that electromagnetic energy can be 
transmitted thereto. 


Apparatus as recited in Claim 10 wherein said nature of 
said fluid is determined as a function of position in 
said material. 


Method for coupling electromagnetic energy into 
materials comprising the steps of: 


forming a measuring tool having a tool face, an 
electromagnetic transmitting antenna and a receiving 
antenna, said transmitting antenna further comprising: 


(a) a coaxial cable connecting means and means to 
transmit electromagnetic energy therethrough; 


(b) a stripline adapter capable of transmitting 
electromagnetic energy from said coaxial cable 
connecting means to a stripline having a metallic 
central strip, said center strip having a front 
end, a flat strip body, a flat strip face, and a 
distal end, said front end electrically connected 
to a center conductor of said stripline adapter, 
said strip face bent at approximately right angles 
to said strip body and having a height measured 
from said right angle bend to said distal end that 
is compatible with a desired frequency coverage; 


(c) a ground plane which extends from said stripline 
adapter to said right angle bend, so that said 
distal end extends away from said ground plane and 
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so that a void exists between said center Strip 
and said ground plane; 


a dielectric largely filling said void, said 
dielectric comprised of a material having a very 
high dielectric constant and a very low energy 
loss, so that said transmitting antenna is 
positioned so that said ground Plane is fixedly 
connected to said measuring tool and said strip 
face is positioned to lie flush with saia tool 
face so that said transmitting antenna can 
transmit electromagnetic energy into said 
material; 


an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in 
electrical contact with said ground plane and said 
Stripline adapter, so that said Strip face is 
nearly centered in the opening created by said 
walls and said ground plane; 


a loss-less, non-conducting material which fills 
in any remaining open space in said enclosure so 
that said non-conducting material forms an 
additional wall that is nearly flat with said 
strip face; 


said receiving antenna comprised in essentially 
the same manner as said transmitting antenna, and 
positioned in said measuring tool in the same 
manner as said transmitting antenna, so that said 
receiving antenna receives said electromagnetic 
energy which has traveled through said material; 
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interconnecting said measuring tool with a means for 
monitoring said electromagnetic energy whereby said 
dielectric properties can be measured; and 


interconnecting said measuring tool with a source of 
electromagnetic energy. 


Method as recited in Claim 25 further comprising a 
means for positioning said tool face near said 
material. 


Method as recited in Claim 25 wherein said 
electromagnetic energy is focused to heat and thereby 
reduce the size of a tumor ina mammal, 


Method as recited in Claim 25 wherein said 
electromagnetic energy is focused to heat and thereby 
destroy a tumor in a mammal. 


Method as recited in Claim 25 wherein said antennas are 


positioned on a belt-mount device. 


Method as recited in Claim 25 further comprising a 
Plurality of receiving antennas. 


Method as recited in Claim 30 further comprising a 
Plurality of transmitting antennas. 


Method as recited in Claim 31 wherein said materials 
having dissimilar dielectric properties are 
mammal tissue and water. 


Method as recited in Claim 25 wherein said broadband 
measurements are taken to determine said dielectric 
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properties as a function of position within 
said material. 


Method for coupling electromagnetic energy to determine 
the quality of a fluid in a material, comprising the 
steps of: 


forming an apparatus having a tool face, an 
electromagnetic transmitting antenna, and a receiving 
antenna, said transmitting antenna further comprising: 


(a) a coaxial cable connecting means and means to 
transmit electromagnetic energy therethrough; 


(b) a stripline adapter capable of transmitting 
electromagnetic energy from said coaxial cable 
connecting means to a stripline having a metallic 
central strip, said center strip having a front 
end, a flat strip body, a flat Strip face, and a 
distal end, said front end electrically connected 
to a center conductor of said Stripline adapter, 
said strip face bent at approximately right angles 
to said strip body and having a height measured 
from said right angle bend to said distal end that 
is compatible with a desired frequency coverage; 


(c) a ground plane which extends fron said stripline 
adapter to said right angle bend, so that said 
distal end extends away from said ground plane and 
so that a void exists between said center strip 
and said ground plane; 


(d) a dielectric largely filling said void, said 
dielectric comprised of a material having a very 
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high dielectric constant and a very low energy 
loss, so that said transmitting antenna is 
positioned so that said ground plane is fixedly 
connected to said measuring tool and said strip 
face is positioned to lie flush with said tool 
face so that said transmitting antenna can 
transmit electromagnetic energy into said 


material; 


an enclosure surrounding said stripline comprising 
four metallic walls, said walls positioned in 
electrical contact with said ground plane and said 
stripline adapter, so that said strip face is 
nearly centered in the opening created by said 
walls and said ground plane; 


a loss-less, non-conducting material which fills 
in any remaining open space in said enclosure so 
that said non-conducting material forms an 
additional wall that is nearly flat with said 
Strip face; 


said receiving antenna comprised in essentially 
the same manner as said transmitting antenna, and 
positioned in said apparatus in the same manner as 
said transmitting antenna, so that said receiving 
antenna receives said electromagnetic energy which 
has traveled through said material; 


interconnecting said measuring tool with a means for 
monitoring said electromagnetic energy whereby said 
nature of said fluid can be determined; and 
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interconnecting Said apparatus with 
electromagnetic energy. 


Method as recited in Claim 31 or 40 


@ source of 


wherein said 


transmitting antenna transmits and said receiving 
antenna receives electromagnetic energy over a 


frequency range of 2 KHz to 4 GHz. 


Method as recited in Claim 25 or 34 


wherein said 


transmitting antenna can alternately function asa 
receiving antenna and said receiving antenna can 
alternately function as a transmitting antenna. 


Method as recited in Claim 36 further comprising a 
belt-mount, said belt-mount substantially conforming to 
the outside of a mammal tissue and holding said 


transmitting and teceiving antennas. 


Method as recited in Clain 37 further comprising a 


Plurality of receiving antennas. 


Method as recited in Claim 38 furthe 
Plurality of transmitting antennas. 


Method as recited in Claim 39 wherei 
antennas are positioned on Said tool 
antennas are positioned on said belt 


Method as recited in Claim 34 wherei 
water, 


Method as recited in Clain 25 or 34 
Plane is no greater than 10 mm in le 


t comprising a 


n some of said 
face and some 
~mount, 


n said fluids are 


wherein said ground 
ngth. 
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Method as recited in Claim 25 or 34 wherein said strip 
face has a height that is no greater than S mm. 


Method as recited in Claim 25 or 34 wherein said 
electromagnetic energy is monitored to provide an 
indication of the amount and distribution of a fluid 
inside mammal tissue. 


Method as recited in Clain 25 or 34 wherein no 
receiving antenna is incorporated. 


Method as recited in Claim 25 or 34 wherein said 
apparatus is implanted inside a mammal, as a radio 
frequency antenna. 


Method as recited in Claim 46 wherein said apparatus 
does not incorporate a receiving antenna. 


Method as recited in Claim 25 or 34 wherein a strip 
transmission line is electrically connected to said 
Stripline, so that electromagnetic energy can be 
transmitted thereto. 


Method as recited in Claim 40 wherein said quantity of 
said fluid is determined as a function of position 
within said material. 
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